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{5,1) 

2V3  (10p5  +3p  -  12p-’)cos(0) 

17 

(5.1) 

“  2V3  (lOp5+3p-12p-'')sin(0) 

18 

(5,3) 

2  VT  (5  p^  -  4  p-^)  cos(3  0) 

19 

{5,  3) 

2V3  (5p5-4p^)sin(30) 

20 

(5,  5) 

2  V3  p^  cos(5  0) 

21 

(5,5) 

2  V3  p**  sin(5  0) 

22 

(6,  0) 

VT  (20p^  +  12p2-30p^-l) 

23 

{6,  2) 

VTT  (15p^’  +6p2  -2Op^)sin(20) 

24 

{6,  2} 

VTT  (15  p^’  +  6  p2  -  20 p‘^)  cos(2 0) 

25 

(6,  4) 

VTT  (6p‘'’-5p'‘)sin(4  0) 

26 

(6,4) 

ViT  (6  p^’  -  5  p'*)  cos(4  0) 

27 

{6,  6) 

VTT  p^’  sin(6  0) 

28 

(6,  6) 

VTT  p^’  cos(6  0) 

29 

{7,  1) 

4  (35  p’  +  30  p-^  -  4  p  -  60  p'’ )  sin(0) 

30 

(7,  1) 

4  (35  p^  +  30 p^  - 4 p  - 60 p^)  cos(0) 

31 

{7,  3) 

4  (2 1  p’  +  1 0  p"*  -  30  p^ )  sin(3  0) 

32 

(7,3) 

4  (21  p’  +  10 p‘^  -  30 p^)  cos(3 0) 

33 

(7,5) 

4(7p’  -6  p^)sin(5  0) 

34 

(7,5) 

4  (7  p’  -  6  p^ )  cos(5  0) 

35 

(7,7) 

4  p’  sin(7  0) 

36 

(7,7) 

4  p’’  cos(7  0) 
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UK  DUAL-USE  ADAPTIVE  OPTICS  PROGRAMME 
Alan  Greenaway 

DERA  Malvern,  St  Andrews  Road,  Malvern  WR14  3PS,  UK 


Abstract 

A  dual-use  prograriime  is  being  established  in  the  UK  in  order  to  exploit  current  research 
on  Adaptive  Optics  for  large  UK  telescopes  and  to  exploit  a  new  programme  on  military 
uses  of  Adaptive  optics.  The  dual-use  programme  will  concentrate  on  the  use  of  Adaptive 
Optics  in  laser  materials  processing.  This  will  address  only  a  restricted  range  of  civilian 
uses  of  Adaptive  Optics  technology.  The  presentation  wiU  describe  the  scope  of  the  dual- 
use  programme,  the  range  of  applications  that  the  programme  is  intended  to  address  and 
the  approaches  that  will  be  adopted  in  order  to  satisfy  the  requirements  of  the  selected 
application  areas. 


Introduction 

In  terms  of  serious  investment  of  manpower  and  finance,  the  UK  is  a  late-entrant  in 
Adaptive  Optics.  The  UK  astronomy  community,  through  the  Particle  Physics  and 
Astronomy  Research  Council  (PPARC)  led  the  way  in  1994  with  a  significant  programme 
intended  to  equip  two  large-aperture  telescopes  with  Adaptive  Optics  correction  and  to 
use  laser  wavefiront-sensing  beacons  to  increase  sky  coverage  beyond  that  which  would  be 
available  if  only  natural  guide-stars  were  used.  A  UK  programme  on  Adaptive  Optics  for 
military  applications  has  just  started  in  April  1997. 

Athhough  Adaptive  Optics  has  been  used  in  military  applications  for  a  considerable  period 
and  in  restricted  commercial  applications  (e.g.  auto-focus)  for  a  similar  length  of  time, 
concerted  attempts  to  apply  Adaptive  Optics  with  more  than  one  degree  of  fi-eedom  to 
industrial  processes  are  still  relatively  few.  Entering  a  high-technology  research  area  well 
behind  other  participants  offers  both  advantages  and  a  serious  challenge.  Because  much 
of  the  fimdamental  technological  developments  (such  as  a  range  of  wavefiront  modulator 
systems,  a  range  of  wavefront  sensor  systems  and  closure  of  feedback  control-loops)  have 
already  been  demonstrated,  a  programme  on  industrial  applications  can  be  kept  within 
modest  costs  and  does  not  require  ‘proof-of-principle’  stages.  However,  the  challenge  is 
to  ensure  that  the  programme  can  deliver  new  and  useful  results  in  a  short  timescale  and 
without  the  extensive  heuristic  experience  that  other  groups  have  gained. 

The  proposed  UK  Dual-use  programme  is  intended  to  bridge  a  gap  between  the  funded 
UK  programmes  in  astronomy  and  military  applications  and  cost-critical  industrial 
applications.  The  areas  of  application  selected  have  been  chosen  to  match,  as  closely  as 
possible,  areas  of  existing  mutual  interest  to  UK  industry,  military  and  academic  bodies. 
The  programme  requires  coordination  of  funding  from  various  and  the  gestation  period  for 
establishment  of  the  funding  has  been  protracted  and  is  still  incomplete. 


Strategy  and  Plans  of  UK  Dual-use  Programme 

Neither  the  existing  PPARC  programme,  nor  the  newly-started  UK  military  programme  on 
Adaptive  Optics,  will  address  directly  the  important  issues  concerned  specifically  with  the 
control  of  lasers  in  materials  processing.  The  industrial  application  of  for  cutting,  welding, 
drilling,  patterning  and  annealing  are  of  significant  interest  to  UK  industry  and  provide 
benefit  for  the  UK  military  and  UK  defence  industry.  The  UK  academic  community  also 
has  a  significant  programme  in  this  area,  funded  through  the  Engmeering  and  Physical 
Science  Research  Council  (EPSRC).  For  these  reasons  the  application  of  Adaptive  Optics 
to  improvement  of  laser-based  materials  processing  was  selected  as  a  suitable  vehicle  for  a 
dual-use  programme,  funded  by  the  Nfinistiy  of  Defence  (MoD),  the  Department  of  Trade 
and  Industry  (DTI),  Industry  and  the  EPSRC. 

For  materials  processing  applications  the  use  of  lasers  as  a  tool  offers  advantages,  over 
mechanical  or  chemical  working  of  the  material,  which  include  speed  of  working,  the 
ability  to  cut/pattem  very  fine  detail  over  large  workpieces,  the  absence  of  tool  wear  and 
vibration  and  the  minimisation  of  induced  thermal  distortion  in  the  workpiece. 

The  interaction  of  the  laser  beam  with  the  material  being  worked  upon  is  a  complex 
process,  involving  not  merely  the  removal  of  material  from  the  area  that  is  being  actively 
worked  but  also  the  control  of  the  ejected  material.  Adaptive  optics  offers  an 
unprecedented  opportunity  for  accurate  and  dynamic  control  of  the  manner  in  which  the 
material  is  processed,  through  high-speed  and  flexible  control  of  the  laser-beam  profile  at 
the  point  of  working. 

Discussions  with  UK  Industry  have  shown  an  encouraging  level  of  interest  in  participation 
in  a  pre-competitive,  dual-use  Adaptive  Optics  programme  to  address  issues  specifically 
concerned  with  laser  materials  processing  Requirements  that  have  emerged  as  offering 
significant  potential  benefit  fi'om  application  of  adaptive  optics  to  industrial  materials 
processing  include: 

•  Accurate  control  of  beam  focus  at  the  point  of  working.  This  is  particularly  important 
when  precise  patterning  and  drilling  of  surfaces  with  significant  surface  curvature  or 
relief  is  required  or  when  precise  control  of  hole  shape  and  depth  is  required. 

•  Accurate  control  of  beam  shape  at  the  point  of  working.  Such  control  can  facilitate  the 
use  of  beam-shapes  tailored  to  the  operation  being  undertaken,  for  example,  focal  spots 
elongated  in  the  direction  of  working  can  be  beneficial  in  cutting  and  welding 
operations,  production  of  a  halo  around  the  working  focus  may  help  to  control 
deposition  of  ejected  material  and  thus  improve  the  finish  achieved,  accurate  control  of 
the  energy  distribution  between  multiple  foci  can  assist  in  accuracy  and  uniformity 
when  parallel  working  is  used  to  increase  throughput. 

•  Intra-cavity  control  of  the  laser  wavefi'ont,  to  increase  the  quality  of  the  beam  exiting 
the  laser  cavity. 

A  prime  concern  that  impinges  on  the  range  of  industrial  applications  to  which  adaptive 
control  of  laser  materials  processing  might  be  put,  is  the  unit  cost  associated  with  adaptive 


control  systems.  A  consensus  opinion  appears  to  be  that  if  the  overall  cost  of  an  adaptive 
system  for  laser-beam  control  is  £100k  or  more  very  few  such  units  would  be  used  and 
these  only  for  very  specialist  or  military  applications.  If  unit  costs  could  be  reduced  to 
£50k  it  appears  that  a  modest  number  of  units  are  likely  to  be  employed,  but  to  ensure 
widespread  application  one  should  aim  for  a  unit  cost  that  does  not  exceed  £10-20k.  At  a 
£10k  unit  cost,  opinion  was  that  adaptive  optics  units  would  be  fitted  to  almost  every  laser 
system  for  industrial  materials  processing  applications. 

Since  the  objective  of  the  proposed  dual-use  programme  is  to  address  successfully  a  wide 
range  of  industrial  applications,  the  reduction  of  unit  costs  to  below  £10k  per  unit,  for 
production  of  1000  or  more  units  of  low-order  extra-cavity  adaptive  optics  systems 
capable  of  handling  industrially-useful  power  levels,  has  been  selected  as  a  prime  objective 
of  the  programme.  Use  of  bi-morph  mirrors  or  use  of  large-scale  spatial  light  modidators 
based  on  liquid  crystals  would  appear  to  offer  the  potential  for  low-cost  waveflont 
modulators.  The  achievement  of  cheap  and  effective  wavefront  sensing  through  a  plasma 
and  cloud  of  ejected  material,  and  coupling  such  sensing  through  cost-effective  electronics 
to  achieve  the  requisite  control,  appears  to  offer  a  significant  challenge. 

Industrial  interest  in  the  programme  spans  the  full  range  of  laser  types  and,  in  order  to 
keep  the  programme  as  generic  as  possible,  it  is  intended  that  the  main  emphasis  will  be  on 
extra-cavity  adaptive  optics  for  beam-shaping  and  control.  For  these  extra-cavity 
applications  it  is  intended  that  the  scientific  output  of  the  programme  will,  once  intellectual 
property  rights  (IPR)  have  been  protected,  be  freely  disseminated  through  the  scientific 
press  and  a  series  of  open  seminars.  It  is  expected  that  use  of  such  IPR  will  be  available 
under  non-exclusive  licence  arrangements.  However,  provision  will  be  made  for 
participants  in  the  programme  to  delay  the  granting  of  licences,  or  to  delay  dissemination 
of  the  results,  if  they  consider  that  this  is  necessary  to  protect  a  significant  commercial 
advantage. 

The  use  of  intra-cavity  adaptive  optics  is  not  so  easily  managed  in  a  collaborative  and  pre- 
competitive  programme,  because  of  the  requirements  of  commercial  confidentiality  in  laser 
manufacture.  Thus  it  is  intended  that  work  applying  technological  developments  within 
the  programme  to  intra-cavity  adaptive  systems  can  be  addressed  by  the  programme 
participants,  but  that  such  studies  would  be  on  a  bi-lateral  or  small-group  basis,  with 
further  funding  being  required  where  significant  effort  was  to  be  put  into  such  studies. 

The  wavefront  sensing  required  to  ensure  accurate  control  of  the  laser  beam  profile  at 
point  of  working  is  a  signiJBcant  topic  to  be  addressed  within  the  programme.  Assessment 
of  focus  may  be  made  between  laser  pulses  using  various  optical  and  mechanical 
techniques.  Current  approaches  tend  to  rely  on  indirect  methods,  including  mechanical 
sensing  of  the  position  of  the  surface  to  be  worked  or  on  maximisation  of  the  scattering 
from  the  plasma  and  cloud  of  ejected  material  produced  when  the  material  is  worked. 

Both  of  these  approaches  have  the  merit  of  simplicity,  but  are  unlikely  to  lead  to  the 
accurate  beam  control  that  is  required  from  this  programme.  More  sophisticated  methods 
using  spectral  diversity  [1]  have  been  shown  to  be  capable  of  sensing  the  penetration  of 


the  keyhole  region  during  laser  welding,  even  with  cw  laser  radiation.  However,  none  of 
these  methods  provide  direct  information  concerning  the  condition  of  the  laser  beam  at  the 
point  of  impact  with  the  workpiece. 

If  the  scattering  from  the  plasma  generated  by  laser  interaction  with  the  material  being 
worked  is  too  strong,  insufBcient  energy  will  reach  the  surface  and  the  plasma  generation 
will  become  self-limiting.  It  thus  seems  probable  that  eflBcient  laser  working  of  the  surface 
will  always  be  done  under  conditions  where  a  large  fraction  of  the  incident  energy  is 
unscattered  before  it  reaches  the  surface  being  worked.  Under  such  conditions  it  should 
be  possible  to  image  directly  the  beam  incident  on  the  work-area  and,  thus,  to  provide 
accurate  feedback  on  the  beam  condition. 

Such  direct  imaging  of  the  laser  beam  at  the  point  of  impact  with  the  workpiece  will  be 
severely  limited  in  contrast  by  emissions  from  the  plasma  and  by  scattering  and  emissions 
from  the  ejected  material  and  the  workpiece.  The  challenge,  therefore,  is  to  exploit 
mechanisms  that  will  enhance  the  contrast  in  the  image.  Amongst  the  mechanisms  that 
will  be  examined,  individually  and  in  combination,  to  facilitate  such  imaging  are: 

•  Polarisation-sensitive  imaging,  using  the  polarisation  to  suppress  emissions  from  the 
workpiece,  plasma  and  ejected  matter  and,  perhaps,  scattering  from  the  ejected  material 

•  Narrow-band  imaging,  to  suppress  emissions  from  the  workpiece,  plasma  and  ejected 
material 

•  Multi-spectral  imaging,  to  provide  a  mechanism  for  discrimination  between  emissions 
and  different  scattering  mechanisms 

•  Coherent  (e.g.  heterodyne)  detection,  to  discriminate  between  emissions  and  different 
scattering  mechanisms. 

Exploitation  of  mechanisms  to  enhance  contrast  is  particularly  important  to  permit  the 
measurement  of  beam-shape  in  the  wings  of  the  beam  profile.  Although,  these  wings  are 
unlikely  to  be  used  for  cutting  or  patterning  they  may  be  used  to  control  deposition  of 
ejected  material.  In  each  case  the  ‘imaging’  or  wavefront  sensing  must  be  achieved  using 
a  low-cost  system,  or  a  system  that  can,  in  the  near  future,  be  re-engineered  to  provide  a 
low-cost  solution. 


Exploitation  and  benefits  expected  from  UK  Dual-use  Programme 

As  far  as  possible  consistent  with  protection  of  IPR  results  from  the  programme  will  be 
published  at  an  early  point.  However,  if  the  programme  is  successful  in  delivering  results 
that  a  programme  participant  considers  to  offer  a  substantial  commercial  advantage, 
provision  will  be  made  to  restrict  the  licence  access  of  non-participants. 

It  is  expected  that  each  of  the  industrial  participants  (currently  5,)  who  between  them  have 
offered  involvement  in  the  programme  at  a  level  of  £470k  over  the  3 -year  programme 
duration,  will  participate  in  the  design  and  subsequent  testing  of  a  prototype,  extra-cavity, 
adaptive-optics  laser-beam  control  system.  The  system  design  will,  as  far  as  practical,  be  a 


common  module  for  the  whole  applications  range  to  be  tested.  For  this  reason  it  is 
unlikely  that  the  system  will  be  hilly  optimised  for  any  one  of  the  applications  to  be  tested, 
but  an  important  aspect  of  the  programme  will  be  to  make  an  assessment  of  the  scope  for 
cost-effective  optimisation  of  the  common  module  for  each  of  the  applications.  Using  this 
approach  is  expected  to  lead  to  the  maximum  commonality  between  applications  and, 
thereby,  the  most  cost-effective,  if  not  fully  optimum,  solution  to  each  requirement. 

Important  potential  is  considered  to  exist  in  the  manufacture  of  large-scale  hybrid  circuits, 
where  the  application  of  adaptive  control  of  laser  heat  treatment  offers  an  unprecedented 
ability  to  initiate  complex  thermal  cycles  individually  tailored  to  the  requirements  of 
individual  components  in  the  circuit. 

The  full  programme  has  been  costed  at  £1.4m  over  its  3-year  duration.  Applications  that 
will  be  addressed  include  drilling,  patterning,  welding  and  annealing  . 

In  terms  of  potential  return  on  the  investment,  it  is  noted  that  patteming/perforation  of 
aerofoils  has  been  estimated  to  have  the  potential  to  effect  a  saving  of  8-10%  in  aviation 
fuel  costs[2].  If  increased  accuracy  in  drilling/patterning  can  contribute  by  1%  toward  the 
level  of  fuel  saving  actually  achieved,  the  total  cost  of  the  programme  would  be  recovered 
by  the  MoD  from  savings  in  aviation  fuel  in  under  one  year.  Additionally,  it  has  been 
suggested  [3]  that  perforated  aerofoils  can  be  used  to  increase  the  stability  of  an  aircraft 
when  flying  slowly  in  delicate  operations  such  as  landing. 

Other  applications  that  will  be  considered  during  the  programme  will  include  ophthalmic 
uses  covering  retinal  imaging  and  therapy,  where  the  increased  resolution  can  provide 
increased  diagnostic  capability  and  reduced  CO-lateral  damage  in  retinal  operations,  and 
potentially,  in  corneal  operations. 


Collaboration 

The  programme  participants  would  welcome  the  opportunity  to  discuss  with  other 
interested  groups  further  collaborations  on  an  equitable  basis. 
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The  Main  Problems  to  Be  Solved  by 
Methods  of  Adaptive  Optics  in  Laser 

industry 


1.  Stabilisation  and  optimisation  of  different  laser 

radiation  parameters. 

2.  Formation  and  maintenance  of  given  intensity 
distribution  of  laser  beam  on  given  surface. 


Applications: 

Laser  microtechnology. 

Laser  etching. 

Laser  heating  technology. 
Laser  ophtalmosurgery. 
Laser  dermatology. 


Different  methods  are  used  to  obtain  the 

desired  result: 

1.  Complex  optical  systems. 

2.  Methods  of  kinoform  flat  optics. 

3.  Methods  of  nonlinear  adaptive  optics. 

4.  Extracavity  adaptive  optical  systems. 

5.  Intracavity  adaptive  optical  systems. 


Design  of  adaptive  optical  system. 


block  of  control 


Aberrations  to  be  reproduced  or  corrected. 


Type  of  laser 

Type  ot 
aberrations 

Amplitude 

Frequency 

range 

Solid  state  YAG  laser 
(output  power  500  W) 

Range  scale: 

astigmatizm 

coma 

spherical  ab. 

P-V  -  8a 
P-V  -  2a 
P-V  -  3/1 

up  to  30Hz 

CO2  technological 
laser 

(output  power  5  kW) 

Range  scale: 
detocus 
astigmatizm 
coma 

P-V  -  10/x 
P-V  -  A/f 
P-V  -  3/1 

up  to  1  Hz 

Semi-passive  bimorph  mirror. 

Main  features:  aperture  -  40  mm,  stroke  -  7/i, 

number  of  actuaturs  -  8,  13,  18, 
mirror  quality  -  0.1 /x  (P-V), 
bandwidth  -  2000  Hz, 
material  -  glass,  quartz,  copper. 


Scheme  of  electrodes  position. 


Photo  of  the  bimorph  mirrors. 


The  interferometric  patterns  of  the  mirror: 
defocns,  astigmatism,  one  electrode  response  profile, 
profile  for  -fllOV  at  odd  electrodes  and  -120V  at  even  ones. 


RMS  of  a  Several  Aberrations 
Approximation  with  17-electrode 

Bimorph  Mirror. 

Defocus  -  0,3% 

Astigmatism  -  0,7% 

Coma  -  8% 

Spherical  aberration  -  3% 


InflueYice  function  of 
different  electrodes. 


El.  N  le. 

Appl .  wol .-150U 
P-V  -  1.41H.  i 
Step  -  0.2H- 


El.  N  9. 

0ppl.  WO1.-150M 
P-U  -  1.58m. 

Step  -  B.ZM. 


El.  N  6. 

Appl.  «h>1.-150V 
P-V  -  0.57m. 
Step  -  0.05m. 


El.  N  1. 

Appl.  vol.— 150V 
P-V  -  1.89m. 
Step  —  0.25j4. 


El.  N  2. 

Appl.  wol. -150V 
P-V  -  1.98m. 
Step  -  B.^M. 


El.  N  3. 

Appl.  wol. -150V 
P-V  -  1.76m. 
Step  -  B.25m. 


Water-cooled  Bimorph  Mirror. 

Main  features:  aperture  -  100  mm,  stroke  -  10/^, 

number  of  actuators  -  18, 
mirror  quality  -  O.lfi  (P-V), 
bandwidth  -  1500  Hz, 
material  ^  copper,  molybdenum. 

Design  of  water-cooled  mirror. 


Photo  of  the  experimental  prototype. 


Formation  of  Spatial  Characteristics 
of  the  Excimer  Laser  Beam. 

Experimental  setup. 


Main  output  parameters: 

duration  of  a  pulse  -  30  ns; 
wavelength  -  308  nm; 
energy  of  pulse  -  30  mj. 

Photo  of  spatial  intensity  distribution. 


Control  of  Copper-Vapor  Laser 
Radiation  by  Adaptive  Mirror. 


Experimental  setup. 


Divergency  decreased  from  0.25  mrad  (without  mirror 
control)  to  0.08  mrod.  Calculated  divergency  -  0.04  mrad. 

Photo  of  spatial  intensity  distribution 

in  the  far  field. 


e.iM 

>  I  ...j 
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Control  of  CO2  Laser  Radiation  by 
Water-cooled  Bimorph  Mirror. 

Exprimental  setup  with  commercial 

laser  TL-5, 


0  0.26  0.52  0.78  t.ms 

Influence  of  the  surface  shape  of  adaptive 
mirror  on  the  output  mode  structure. 


I jarb .un . 


5 


Schematic  setup  of  the  CW  CO2  laser  (A,  =  10.6  pm) 


Construction  of  bimorph  deformable 

mirror. 


1-  reflecting  coating 

2-  glass  plate 

3-  piezoceramic  actuator  disks 

4-  grounded  electrode 

5-  common  electrode  e1 

6-  segmented  ring  electrodes  e2,  e3 

7-  baseplate 


Scheme  of  segmented  electrodes 


The  reflection  from  both  sides  of  mirror  (1)  created  a  Fizeau  interferometric  fringe  pattern 
registrated  by  CCD  camera  (10)  and  then  developed  with  a  special  program  in  IBM  computer 
(1 1).  The  output  laser  beam  profile  was  registrated  by  photodiode  (5). 

Intracavity  light  power  is  about  160  W. 


Laser  beam  profile  vs.  cavity  mirror  deformation 
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Temporal  behavior  of  the  intracavity  laser 
mirror  under  light  heating 


4 - ► 

2.5  mm 

These  interferogramms  presents  the  time 
change  of  intracavity  mirror  surface  under 
light  heating.  Time  step  between  pictures 
At=0.5  sec.  First  picture  corresponds  to  the 
moment  of  the  starting  laser  generation. 
Transmission  coefficient  of  mirror  under  test 
was  P=  0.2%. 


This  curve  presents  the  time 
dependance  of  laser  cavity 
mirror  deformation  Iq.  The 
whole  process  of  reaching 
stationary  regime  takes  about  2 
seconds.  This  curve  is  well 
approximated  by  following 
expression: 


Pictures  1  -  5  showes  the  change  of  the 
curveture  radius,  while  on  pictures  6-9  one 
can  see  only  a  slight  fluctuation  of  the  fringe 
pattern  on  the  edge  of  the  interferogramm. 


Io=Ios(l  “  exp(-t/T)) 
where  t=0.7  sec. 


Conclusion 


O 

CO  ’S 

g3  S 

ti  tC  ” 

P,  -K  :§ 

§)•§  .S 
•S  2  S 

*  'S)^ 

^  s  t: 

+-•  (U  2 
rP  rS 
O  ttf 

4J  b  ‘-I 

^  g 


(N 


CO 

(D 


CO 


intracavity  beam  power  160  W.  This  time  depends  strongly 
on  the  quality  of  the  dielectric  mirror  cover. 


Compensation  of  thermal  lens  aberrations. 

Experimental  setup. 


Amplitude  of  thermal  lens  aberrations. 

A1,  A2,  -  astigmatizm,  Kl,  K2  -  coma,  SP  -  spherical  aberration 
Pump  current  -  26  A.  Results  are  given  in  microns. 


A1 

A2 

Kl 

K2 

SP 

Before  correction 

-0.22 

1.90 

-0.27 

-0.24 

-0.06 

After  correction 

-0.02 

0.41 

-0.17 

-0.05 

0.03 

IntepfepogpaH  of  active  elenent 


without  thermal  aberrations  of  correction  of  this 

lens  thental  lens  aberrations  by 

adaptive  mirror 


Control  and  Formation  of  the 
Solid  State  Laser  Radiation. 

Experimental  setup. 


Main  output  parameters  of  laser  radiation: 
output  power  in  multymode  regime  -  60  W; 

in  single  mode  regime  -  7  W, 
pump  current  -  32  A,  5  kW  lamp. 

'Decreasing  the  divergency  of  a  laser  beam. 


A  -  loop  opened 

B  -  loop  closed 
(divergency  reduced 
three  times) 


Far  field  distribution. 


Different  far  field  intensity 
distribution  of  a  solid  state  laser 

radiation. 


Single  node  regine. 


Hultynode  regine. 
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Laser  beam  intensity  profile  control  by  adaptive  optics 


Koshichi  Nemoto,  Takashi  Fujii,  Naohiko  Goto,  Takuya  Nayuki 
Komae  Research  Laboratory,  Central  Research  Institute  of  Electric  Power  Industry,  11-1  Iwado-Kita  2-chome,  Komae-shi, 
Tokyo  201  JAPAN,Fax:intl.+81  3  3488  6697,  e-mail:  nemoto@criepi.denken.or.jp 

and  Yoshi-kazu  Kanai 

Leonix  corporation,  4-15-3  Higashi-Koigakubo,  Kokubunji,  Tokyo  185  JAPAN 

Laser  irradiation  efficiency  has  a  strong  affect  on  the  cost  of  laser  isotope  separation.  We  are  developing  a  laser 
beam  intensity  profile  control  device  based  on  coordinate  transformation  in  order  to  improve  the  laser  irradiation  efficiency.  In 
this  paper,  we  attempt  to  obtain  the  optimum  surface  profile  of  the  deformable  mirror  automatically. 

The  Fresnel  equation,  which  describes  the  beam  propagation  characteristics,  shows  that  the  complex  amplitude  at 
the  location  of  the  image  can  be  controlled  by  modifying  the  phase  distribution.  Iherefore,  if  we  can  generate  the  phase 
distribution  necessary  for  intensity  profile  transformation  using  a  deformable  mirror,  we  can  transform  an  arbitrary  intensity 
profile  to  other  profiles.  This  method  will  transform  beams  to  various  profiles  as  long  as  the  deformable  mirror  can  produce  the 
necessary  surface  profiles  and  the  compatibility  condition  is  satisfied. 

The  feedback  control  loop  is  shown  in  Fig.  1 .  The  beam  reflected  by  a  deformable  mirror  is  monitored  on  the  screen 
by  a  CCD  camera.  The  image  data  is  transfeixed  from  the  image  processor  to  a  computer  through  aGPIB  bus,  and  evaluated  by 
the  computer.  The  computer  calculates  the  siq^plied  voltage  of  nine  actuators  based  on  the  evaluation  values  of  the  profiles.  The 
deformable  mirror  is  designed  to  transform  a  Gaussian  profile  to  a  rectangular  one,  and  has  nine  moving  points  and  four  fixed 
points  1).  When  we  control  the  laser  beam  intensity  profile,  the  relation  between  the  surface  profile  of  the  deformable  mirror 
and  the  intensity  profile  of  the  reshaped  laser  beam  is  more  complicated  than  that  of  phase  compensation  because  the  complex 
amplitude  of  the  reshaped  beam  is  the  Fourier  or  Fresnel  transformation  of  the  beam  reflected  on  the  deformable  mirror. 
Moreover,  the  characteristics  of  actuators  are  not  linear  and  depends  on  the  stress.  Thus,  the  transfer  function  of  this  system  is 
complicated.  We  adopted  Genetic  Algorithm  (GA)  for  controlling  a  deformable  mirror  because  we  can  treat  this  system  as  a 
black  box  using  GA.  GA  is  one  of  the  optimum  value  calculation  algorithms  for  obtaining  the  maximum  or  minimum  value  of  a 
function  f(x)  in  restricted  region.  The  only  feedback  value  is  the  fitness  value,  which  describes  the  degree  of  similarity  of  the 
beam  profile  to  a  rectangular  contour. 

Figure  2  (a)  and  (b)  show  the  initial  Gaussian  beam  profile  and  the  optimum  beam  profile  reshaped  to  a  rectangular 
contour,  respectively.  The  optimum  intensity  profile  was  obtained  within  the  limits  of  the  deformable  mirror  system.  In  every 
trial,  each  maximum  fitness  value  converges  to  the  same  value  after  more  than  30  iterations.  Moreover,  each  average  fitness 
value  also  converges  to  the  same  level.  However,  since  the  mutation  rate  vras  set  to  large  value,  the  average  value  did  not 
converged  to  the  same  value  as  the  maximum  value.  The  spatial  coherency  of  the  beam  is  shown  to  be  conserved  after  intensity 
profile  transformation,  since  we  can  observe  the  interference  fringes  between  the  reshaped  beam  and  the  reference  beam. 

Reference  1)  K.  Nemoto  et.al.  Optics  Lett.,  21(3),p. 168(1995) 


Deformable  Mirror 


(a)  (b) 

Fig.2  Initial  (a)  and  obtained  beam  profile  (b) 


Fig.l  Schematic  diagram  of  the  control  system 


Development  of  an  Adaptive  Optics  System  for  SILVA 


P.-Y.  Madec,  V.  Michaud,  D.  Rabaud,  F.  Charbonnier,  R.  Deron,  L.  Rousset-Rouviere 
Office  National  d'Etudes  et  de  Recherches  Aerospatiales 
B.P.  72, 92322  Chatillon  -  France 

C.  Gallou,  E.  Delavaut,  P.  Dupuy,  G.  Delarue 
Commissariat  a  I'Energie  Atomique 
Centre  d'Etudes  de  Saclay  -  91 191  Gif  sur  Yvette  -  France 


1.  Introduction 

To  propose  an  alternative  to  currently  used  uranium 
enrichment  processes,  the  Combustible  Cycle  Department 
of  the  CEA  (the  French  Atomic  Energy  Commission)  is 
developing  a  process  called  SILVA,  a  French  acronym  for 
Atomic  Vapour  Isotope  Separation  by  Laser.  To  optimise 
the  performance  of  the  SILVA  prototype,  ONERA  and 
CEA  are  developing  and  implementing  an  Adaptive 
Optics  system  designed  to  reshape  the  wavefront  of  the 
dye  laser  beams  used  to  ionize  the  uranium  vapour. 

The  Adaptive  Optics  system  consists  of  a  deformable 
mirror  with  48  actuators  and  a  wavefront  sensor  (WFS) 
including  21x21  analysis  points.  From  the  measurements 
supplied  by  the  analyzer,  a  real-time  computer  determines 
the  action  to  be  applied  to  the  deformable  mirror.  The 
control  computer  system  architecture  is  based  on  a 
modular  and  upgradable  design.  The  Adaptive  Optics 
system  developed  can  thus  easily  be  reconfigured  and 
especially  adapted  to  the  SILVA  improvements  expected 
during  the  optimization  phases. 

This  paper  first  recalls  the  principle  of  the  SILVA 
process,  then  describes  the  Adaptive  Optics  system 
developed  and  the  characteristics  of  its  main  components. 
The  system  performance  simulated  by  dedicated  software 
developed  by  ONERA  and  CEA  is  indicated.  Finally,  first 
experimental  characterization  results  are  presented. 

2.  SILVA  Principle 

The  SILVA  process  is  now  under  study  at  CEA.  It 
could  ultimately  replace  gas  diffusion  for  uranium 
enrichment  with  isotope  235  [2.1,  2.2,  2.3]. 

The  SILVA  process  includes  two  subsystems  : 

•  The  separator  which  generates  the  atomic  vapour 

•  The  laser  system  which  supplies  the  photons  used  to 
selectively  ionize  the  uranium  235. 

Selective  ionization  of  uranium  235,  carried  out  in 
three  steps,  requires  using  a  high-power  laser  beam 
obtained  by  multiplexing  several  visible  wavelengths. 
The  laser  system  outputs  this  beam  which  is  reflected  in 
the  separator  by  a  set  of  mirrors  to  illuminate  the  vapour. 

For  optimal  ionisation  of  all  the  uranium  235  atoms, 
the  beam  illumination  profile  must  remain  sufficiently 
uniform  throughout  its  propagation.  The  beam  wavefront, 
which  plays  a  major  role  in  the  variation  of  the 
illumination  profile,  must  therefore  be  controlled.  This  is 
done  using  an  Adaptive  Optics  system  with  a  deformable 
mirror,  a  wavefront  sensor  and  real-time  computer. 

Adaptive  optics  systems  were  designed  and  have 
been  in  use  for  several  years  at  ONERA  for  military  and 


civilian  applications.  A  dedicated  system  has  been  adapted 
to  correction  of  a  dye  laser  beam  wavefront.  This  laser, 
consisting  of  an  oscillator  and  dye  amplifiers,  shows 
particular  spatial  characteristics.  The  rectangular  geometry 
cross-section  laser  beam  exhibits  high  spatial  frequency 
intensity  non-uniformities,  in  particular  a  large  over-shoot 
in  intensity  near  the  edges,  which  may  introduce 
wavefront  measurement  errors.  This  problem  is  mentioned 
in  Section  4.2. 

In  the  future,  several  from  those  systems  could  be 
required  to  shape  the  SILVA  process  beam,  at  least  one 
per  color. 

3.  General  Description  of  an  Adaptive  Optics  System 

An  Adaptive  Optics  system  is  a  servo  closed-loop 
opto-mechanical  set-up  which  compensates  for  the  phase 
errors  of  an  optical  wave  in  real  time  [3.1]. 

This  is  made  possible  by  the  use  of  a  deformable 
mirror  located  in  the  optical  path.  The  mirror  consists  of  a 
thin  glass  plate  and  actuators  located  behind  it.  The 
driving  of  these  actuators  can  control  the  profile  of  the 
mirror's  reflective  surface  to  compensate  for  the  phase 
errors  in  the  reflected  optical  wave. 

The  second  major  component  of  an  Adaptive  Optics 
system  is  the  wavefront  sensor.  This  component  measures 
the  phase  errors  to  be  compensated  for  by  the  deformable 
mirror.  There  are  several  types  of  wavefront  sensors.  The 
one  chosen  for  the  SILVA  system  is  based  on  the 
Hartmann-Shack  principle  [3.2]. 

Finally,  to  control  the  deformable  mirror  by  a 
feedback  loop,  it  is  necessary  to  have  a  computer 
architecture  to  process  the  wavefront  sensor  signals  in  real 
time  to  determine  the  control  voltages  to  be  applied  to  the 
deformable  mirror. 

The  specifications  of  these  components  are 
determined  by  the  spatial  and  temporal  characteristics  of 
the  errors  to  be  corrected  and  the  correction  accuracy  to  be 
achieved  [3.2,  3.3].  In  particular,  the  number  of  degrees 
of  freedom  of  the  Adaptive  Optics  system,  generally 
determined  by  the  number  of  actuators  of  the  deformable 
mirror,  is  given  by  the  spatial  spectrum  of  the  phase  errors 
to  be  compensated  for.  Similarly,  the  desired  correction 
bandwidth  is  given  by  the  temporal  spectrum  of  the  errors 
to  be  corrected. 

Preliminary  measurements  made  on  the  dye  laser 
beam  of  the  prototype  SILVA  system  showed  that  the 
main  phase  errors  to  be  compensated  for  were  defocusing 
and  astigmatism,  as  well  as  aberrations  at  a  higher  spatial 
frequency  whose  amplitude  was  sufficient  to  significantly 
degrade  the  quality  of  the  dye  laser  beams.  An  optical 


quality  better  than  "kJlO  rms  is  required  for  the  corrected 
wavefront.  This  specification,  to  which  is  added  the 
requirement  for  adapting  to  the  rectangular  geometry  of 
the  optical  pupil  of  SILVA,  led  to  use  a  deformable  mirror 
with  8x6  actuators  in  a  square  mesh. 

Since  so  far  we  do  not  have  a  sufficiently  fast 
wavefront  sensor,  the  characteristic  variation  times  of  the 
errors  to  be  corrected  are  not  known.  At  the  request  of 
CEA,  the  Adaptive  Optics  system  developed  should 
privilege  accurate  spatial  analysis  of  the  errors.  The 
correction  bandwidth  obtained  is  not  a  crucial 
requirement.  However,  the  real-time  computer  system 
designed  must  be  sufficiently  powerful  and  reconfigurable 
to  be  able  to  evolve  toward  large  bandwidth  Adaptive 
Optics. 

4.  Description  of  the  Components 
4.1.  Deformable  Mirror 

The  deformable  mirror  obtained  for  SILVA  was 
manufactured  by  CILAS  [4.1]  (Figure  4.1).  It  includes  8  x 
6  actuators  located  in  a  square  mesh  with  a  spacing  of  9.5 
mm.  Each  motor  has  a  mechanical  excursion  of  ±S  pm 
obtained  with  a  control  voltage  of  ±400  V.  The  first 
mechanical  resonance  frequency  is  higher  than  14  kHz. 

The  mirror  surface  roughness  is  very  weak.  The 
errors  at  rest  are  mainly  due  to  defocusing  and 
astigmatism.  In  closed-loop  operation,  the  mirror  flatness 
is  better  than  X/10  peak-to-peak  (X/45  rms). 


Figure  4. 1  :  Deformable  mirror  with  48  actuators  made 
by  CILAS;  the  mechanical  cover  was  removed  to  show 
the  piezoelectric  actuators. 

This  type  of  mirror  with  a  backside  water-cooling 
system  of  its  reflective  surface,  can  operate  in  presence  of 
laser  beam  power  up  to  5  kW.  The  water  flow  is 
generated  by  nozzles  visible  in  the  photograph  of  Figure 

4.1. 

4.2.  Wavefront  sensor 

The  Hartmann-Shack  analyser  principle  is  based  on 
geometrical  optics  [4.2].  The  wavefront  is  sampled  by  a 
matrix  of  microlenses,  each  defining  an  analysis  subpupil. 
The  phase  errors  then  occur  as  a  simple  tilt  of  the 
wavefront.  The  tilt  causes  shifting  of  the  image  formed  in 
the  focus  of  the  subpupil  with  respect  to  a  reference 
position  given  by  the  position  of  each  image  for  a 
perfectly  flat  wavefront.  All  these  shiftings  are  measured 


to  determine  the  local  wavefront  tilt  in  each  subpupil. 
These  measurements  can  be  used  either  to  reconstruct  the 
wavefront  by  the  help  of  specific  algorithms  or,  to  control 
a  deformable  mirror,  which  is  also  a  form  of 
reconstruction. 

The  local  shifting  of  the  images  is  measured  by  a 
CCD  camera  located  in  the  focus  plane  of  the  microlenses 
whose  video  output  signal  is  digitized  and  processed. 
Processing  consists  simply  of  calculating  the  position  of 
the  center-of-gravity  of  the  image  formed  in  the  focus 
plane  of  each  subpupil. 

The  wavefront  sensor  developed  for  SILVA  must 
satisfy  two  requirements  :  first  to  be  able  to  control  a 
deformable  mirror  to  provide  an  Adaptive  Optics  system 
and  secondly,  in  the  absence  of  correction,  characterise 
precisely  the  wavefront  errors  of  the  dye  laser  system. 
This  second  requirement  makes  it  necessary  to  define  an 
analyser  whose  spatial  resolution  is  as  high  as  possible 
and  which  therefore  has  a  large  number  of  subpupils.  To 
optimise  the  control  of  the  deformable  mirror,  the 
geometry  of  the  latter  must  be  adapted  to  that  of  the 
analyser.  It  is  preferable  to  have  a  whole  number  of 
subpupils  between  each  actuator.  To  meet  these  two 
requirements,  one  sensor  with  21  x  15  subpupils  was 
chosen. 

The  choice  of  the  analyser  CCD  sensor  was  guided 
by  the  number  of  pixels  required  in  each  subpupil.  This 
number  dimensions  the  field  of  analysis  of  each  subpupil, 
i.e.  the  sensor  measurement  dynamic  range.  A  large 
number  of  pixels  allows  analysis  of  large  phase  errors. 
However,  too  many  pixels  is  penalising  for  the  response 
time  of  the  Adaptive  Optics  system.  Indeed,  the  more 
pixels  in  the  sensor,  the  higher  the  feedback  loop  control 
frequency  and  therefore  the  smaller  the  correction 
bandwidth.  A  256  x  256  pixels  CCD  sensor  made  by 
DALSA  company  capable  of  a  rate  of  150  images  per 
second  was  chosen.  TTiere  are  therefore  12  x  12  pixels  for 
each  subpupil. 

Finally,  after  settling  the  physical  specifications  of 
the  analyser,  we  needed  to  evaluate  its  measurement 

accuracy  [4.3].  Where  is  the  analyser  measurement 
noise  variance  expressed  in  rad^  we  have: 

2  _  5?  1  /Nt\^  Ns^ 

^  “  2  <N>  \ND/  3  <n>2  Nd^ 

where: 

•  <N>  is  the  average  number  of  photons  detected  per 
subpupil  and  per  frame,  expressed  in  electrons  per 
pixel  per  frame 

•  Ce  is  the  standard  deviation  of  the  CCD  electronic 
read  noise  expressed  in  electrons  per  pixel  per  frame 

•  Nx  is  the  full  width  at  half  maximum  of  the  image 
profile  in  the  focus  plane  of  the  microlenses 
expressed  as  a  number  of  pixels 

•  Np  is  the  full  width  at  half  maximum  of  the 
diffraction  limited  profile  in  the  focus  plane  of  the 
microlenses  expressed  as  a  number  of  pixels 

•  Ns^  is  the  number  of  pixels  used  to  compute  the 
center-of-gravity . 


Since  the  dynamic  range  is  defined  as  the  ratio 
between  the  maximum  number  of  electrons  that  can  be 
cumulated  in  each  pixel  and  the  standard  deviation  of  the 
read  noise,  the  DALSA  sensor  used  has  a  dynamic  range 
of  500.  Furthermore,  at  half  the  CCD  saturation,  <N>  is 
equal  to  approximately  35,000  electrons  per  pixel  per 
frame.  For  SILVA,  the  errors  to  be  measured  are  mainly 
defocusing  and  astigmatism.  For  instance,  the  full  width 
at  half  maximum  of  the  image  in  the  focus  of  each 
microlens  is  limited  by  diffraction.  We  therefore  have 
Nj/Nd  =  1.  For  Shannon  sampling  of  each  image,  we 
determine  the  focal  length  of  the  microlenses  such  that  Np 
=  2.  Finally,  the  number  of  pixels  used  to  calculate  the 
center-of-gravity  is  set  to  around  10. 

Based  on  these  characteristics,  the  measurement 
accuracy  is  Xyi60  rms,  where  X  is  the  mean  wavelength  of 
the  wavefront  sensor,  i.e.  approximately  0.6  pm.  It  should 
be  noted  that  it  is  very  difficult  to  achieve  such  an 
accuracy  experimentally,  mainly  due  to  local  turbulences 
in  the  test  room  and  vibrations. 

Finally,  although  it  can  be  formally  demonstrated 
that  the  calculated  position  of  the  center-of-gravity  of  the 
image  in  each  subpupil  is  equal  to  the  spatial  mean  of  the 
first  derivative  of  the  corresponding  wavefront,  it  should 
be  noted  that  this  result  is  true  only  if  the  intensity  is 
uniform  in  the  subpupil.  However,  one  of  the  features  of 
the  SILVA  dye  laser  beams  is  to  exhibit  non-uniform 
intensity  at  spatial  scales  smaller  than  a  subpupil. 

The  wavefront  sensor  measurement  accuracy  was 
simulated  at  ONERA.  Typical  phase  errors  of  the  SILVA 
system  were  considered.  The  intensity  distribution  used 
was  based  on  experimental  measurements  made  on  the 
CEA  bench  facilities.  The  conclusion  of  the  simulations  is 
that  the  intensity  variations  encountered  are  not  a  limiting 
factor  for  characterisation  of  the  phase  errors  expected  on 
the  SILVA  dye  laser  beams.  Criteria  were  defined  to 
evaluate  the  phase  measurement  accuracy  that  can  be 
achieved  for  a  given  intensity  profile.  These  criteria  are 
based  on  the  spatial  spectrum  of  the  intensity  distribution 
in  the  pupil. 

4.3.  Real-Time  Computer 

The  computer  was  derived  from  developments 
conducted  at  ONERA  in  collaboration  with  the  French 
company  SHAKTI  [4.4]  for  space  observation 
applications.  Its  architecture  is  capable  of  controlling  an 
Adaptive  Optics  system  with  a  large  bandwidth.  Using 
two  CPU  cards  with  a  total  of  12  Texas  Instrument 
TMS320C40  Digital  Signal  Processors,  it  is  capable  of 
controlling  100  actuators  at  frequencies  up  to  2500  Hz. 

However,  for  applications  requiring  less  computer 
power,  it  is  possible  to  use  a  powerless  dedicated  version 
of  these  cards  which  can  be  inserted  in  an  ordinary  PC. 

The  computer  includes  all  the  electronics  used  to 
compute  the  deformable  mirror  control  voltages  from  the 
analogue  video  signal  output  by  the  camera.  Its  main 
components  include: 

•  A  remote  digitization  unit 

<  A  digital  video  signal  equalizer  unit 

•  A  real-time  CPU  unit 


•  A  video  display  unit. 

Developed  entirely  by  SHAKTI,  it  is  capable  of 
controlling  the  48  actuators  of  the  deformable  mirror  using 
the  video  signal  from  the  DALSA  CCD  camera. 

The  maximum  control  frequency  achievable  with  this 
architecture  is  30  Hz  which  allows  a  correction  bandwidth 
of  a  few  Hz.  This  computer  could  achieve  a  control 
frequency  of  500  Hz  without  any  modification,  simply  by 
using  a  64  X  64  pixel  CCD  sensor. 

The  architecture  developed  includes  two  racks.  The 
first  contains  the  digitization  unit  and  the  second,  to  the 
VME  standard,  contains  the  equalizer,  CPU  and  graphic 
display  cards.  Both  racks  are  controlled  via  a  workstation 
providing  the  user  interface.  The  workstation  is  used  for 
initialising  the  cards  and  having  them  perform  the 
functions  required  for  management  of  the  system. 

Remote  Digitization  Unit 

This  unit  consists  of  a  rack  in  single  Eurocard 
format.  It  is  capable  of  digitizing  the  video  signal  from 
any  camera  on  12  bits  at  rates  up  to  10  Mpixels/s.  The 
digitization  noise  is  less  than  one  bit.  The  unit  performs 
many  operations,  including  dynamic  sensor  windowing, 
offset  adjustment  and  digitizing  by  double  correlated 
sampling.  The  digitization  clocks  are  set  from  the 
workstation  by  software. 

After  being  digitized,  the  signal  is  sent  to  the  VME 
computer  containing  the  equalization  card  and  CPU  card. 
The  digital  video  signal  can  therefore  be  transmitted  over 
long  distances  without  any  loss  of  quality,  even  in  severe 
electromagnetic  environments. 

Video  Signal  Equalizer  Card 

So  as  not  to  introduce  any  bias  on  calculation  of  the 
intensity  center-of-gravity  in  each  subpupil,  it  is  necessary 
for  all  the  pixels  to  have  the  same  flux  response.  If  I(i,j)  is 
the  raw  signal  supplied  by  the  pixel  on  the  i**’  line  and  j**’ 
column,  I(i,j)  depends  on  the  received  flux  F(i,j)  as 
follows: 

I(i,j)  =  G(i,j)F(i,j)-HO(i,j) 

where  G(i,j)  is  the  pixel  sensitivity  and  0(i  j)  its  offset. 

The  equalizer  card  is  used  to  obtain  a  uniform 
response  on  all  the  pixels  by  subtracting  an  offset  plane 
representing  all  the  0(i,j)  from  the  image  consisting  of  all 
the  I(i,j)  then  multiplying  the  result  by  a  gain  plane 
consisting  of  all  the  1/G(i,j).  All  these  operations  are  done 
in  real  time  The  offset  and  gain  planes  can  be  recorded 
experimentally.  They  are  loaded  into  a  memory  space  on 
the  card  which  performs  equalization  in  real  time  at  rates 
up  to  10  Mpixels/s. 

After  being  equalized,  the  video  signal  is  thresholded 
to  eliminate  any  unilluminated  pixels  whose  only 
contribution  to  calculation  of  the  center-of-gravity  is  the 
read  noise.  Therefore,  only  around  10  pixels  are 
preserved  for  calculation. 

The  equalized  and  thresholded  video  signal  is  then 
sent  to  the  CPU  card  via  a  private  video  bus  located  on  the 
card  front  edge.  The  data  rate  on  the  bus  can  be  up  to  10 
Mwords/s.  This  private  bus  frees  the  VME  bus  to  which 


the  card  is  connected  and  which  is  used  only  to  set  the 
parameters,  such  as  the  gain  and  offset  planes  and  acquire 
the  experimental  data  in  real  time  (analyser  measurements, 
deformable  mirror  commands). 

Display  Card 

This  card  is  used  for  video  display  of  the  digitized 
and  equalized  video  signal  supplied  by  the  equalizer  card. 
So  as  not  to  slow  down  the  private  video  bus,  this  card  is 
located  in  wiretap  mode:  it  merely  extracts  data  from  the 
bus  as  needed  without  disturbing  transmission.  It  is  one  of 
the  basic  components  for  adjusting  the  Adaptive  Optics 
system. 

CPU  Card 

This  card  is  at  the  heart  of  the  computer  architecture. 
It  provides  the  computation  power  required  for  real-time 
control  of  the  loop. 


Figure  4.2.  :  SHAKTI  motherboard  equipped  with  a 
CPU  module  with  3  DSPs  and  a  64-output  DAC  module 

It  consists  of  a  motherboard  on  which  various 
modules  can  be  mounted  (Fig.  4.2).  The  motherboard 
dialogues  with  the  VME  bus  and  retrieves  the  digitized 
data  transferred  via  the  private  bus.  It  includes 
programmable  logic  components  so  that  it  can  accept  any 
type  of  data  format  from  the  video  bus  on  its  input  and 
send  the  computed  data  in  any  format  on  its  output.  The 
power  of  this  card  resides  in  the  possibility  of  completely 
reconfiguring  it  by  dynamically  downloading  new 
programmes  into  its  programmable  logic.  This  means  that 
the  same  computer  architecture  can  be  easily  used  for 
different  applications.  This  capability  has  been  widely 
used  in  recent  years  at  ONERA. 

Different  modules  can  be  mounted  on  the 
motherboard.  For  SILVA,  two  CPU  modules  are  used, 
each  including  three  Texas  Instrument  TMS320C40  DSPs. 
The  first  module  computes  the  centres  of  gravity.  The  first 
DSP  is  used  for  this  computation  as  such,  the  second 
stores  the  computed  data  in  a  buffer  and  the  third  stores 
the  equalized  images  from  the  wavefront  sensor.  These 
buffers  are  accessible  by  both  the  C40  DSPs  and  the  VME 
bus  which  allows  the  user  to  retrieve  the  data  in  real  time 
for  analysis.  This  operation  is  possible  even  when  the 
system  is  operating  in  closed-loop  mode.  The  second 
module  computes  the  deformable  mirror  control  voltages 
from  the  results  supplied  by  the  first  module.  The 
computation  consists  of  multiplying  the  precalibrated 
control  matrix  discussed  below  by  a  vector  containing  the 


analyser  measurements.  The  multiplication  is  performed 
by  the  first  DSP,  and  the  second  DSP  stores  the  computed 
control  voltages  in  a  buffer  accessible  to  user. 

The  card  performance  is  limited  by  the  center-of- 
gravity  computation.  This  operation,  the  most  costly  in 
CPU  time  and  data  transfer  time,  can  only  be  performed 
by  a  single  DSP.  Each  of  the  processors  is  potentially 
capable  of  processing  5  Mpixels/s.  However,  for  SILVA, 
the  maximum  rate  achieved  is  limited  to  2  Mpixels/s, 
which  corresponds  in  our  case  to  a  frame  frequency  of  30 
Hz.  This  limitation  is  due  to  the  real-time  software  used 
to  compute  the  centres  of  gravity.  If,  during  testing  on  the 
Adaptive  Optics  system,  it  appears  that  the  specified 
bandwidth  is  too  small,  the  software  will  be  first  optimised 
to  achieve  the  maximum  performance  capabilities  of  the 
component.  If  the  bandwidth  is  still  too  small,  the  use  of  a 
sensor  with  fewer  pixels  will  have  to  be  considered.  For 
instance,  a  CCD  camera  with  64  x  64  pixels  would  allow  a 
control  frequency  of  500  Hz  to  be  achieved  without 
optimization. 

The  third  module  mounted  on  the  motherboard  is  a 
DAC  with  64  converters.  Each  operates  on  12  bits  with  a 
conversion  speed  of  2  V/ps  and  one-half  bit  of  noise.  The 
output  conversion  dynamic  range  is  continuously 
adjustable  up  to  -i- 10  V. 

4.4.  Control  Software 

The  software  architecture  used  to  control  the  SILVA 
Adaptive  Optics  system  was  entirely  developed  by 
ONERA.  This  architecture  includes  function  libraries 
available  in  servers  accessible  on  the  TCP/IP  network  and 
a  graphic  user  interface.  The  software  can  be  executed  on 
a  workstation  or  a  PC  tyf)e  computer. 

Programming  of  the  function  libraries  is  structured 
on  three  different  levels: 

•  The  lowest  level  includes  the  functions  close  to  the 
machine.  The  programmes  depend  on  the  computer 
hardware  used.  Programmes  are  either  in  C  or  C40 
assembler  and  can  only  be  written  by  real-time 
computing  specialists. 

•  The  medium  level  includes  the  basic  functionalities 
required  for  an  Adaptive  Optics  system,  such  as 
measurement  acquisition  and  mirror  control.  The 
programmes  are  sequences  of  low  level  functions. 
They  are  in  C  or  C++  and  can  be  written  by  software 
developers. 

•  The  highest  level  includes  tools  for  using  the  Adaptive 
Optics  system,  such  as  control  matrix  computation  or 
offset  plane  acquisition.  The  programmes  in  this  case 
are  sequences  on  medium  level  functions.  Written  in 
TCL,  an  interpreted  language  in  the  public  domain, 
they  use  function  servers  created  from  medium  level 
libraries.  They  can  be  developed  by  engineers 
specialised  in  Adaptive  Optics. 

The  graphic  user  interface  is  used  to  control  the 
complete  system  from  a  graphic  menu  using  the  mouse. 
Written  in  TCL-TK,  which  is  the  graphic  extension  of 
TCL,  it  is  a  sequence  of  high  level  functions  which  allow 
the  user  to  select  the  physical  parameters  to  be  set  for  the 
system. 

This  programming  philosophy  has  many  advantages. 
Programming  on  several  levels  means  that  only  part  of  the 
software  needs  to  be  rewritten  in  case  of  upgrade  of  the 


equipment  or  modirication  of  the  Adaptive  Optics  system 
specifications  or  even  a  change  of  application. 

It  is  obvious  that  this  saves  considerable 
development  time  and  thereby  reduces  the  cost.  It  also 
allows  the  different  levels  to  be  developed  in  parallel  and 
requires  only  limited  number  of  specialists  in  real-time 
computing.  Finally,  the  use  of  an  interpreted  language 
and  function  servers  makes  it  easy  for  the  user  to  modify 
the  interface  or  test  new  capabilities  by  constructing  his 
own  management  application,  leaving  an  enormous 
amount  of  space  for  imagination  and  creativity. 

5.  Control  Algorithms 

5. 1 .  Control  Matrix  Computation 

All  the  physical  phenomena  relating  the  wavefront 
sensor  to  the  deformable  mirror  are  assumed  to  be  linear. 
We  ca  take  the  example  of  an  aberrant  wave  for  which  the 
analyser  outputs  a  measurement  vector.  The  aim  is  to  find 
a  control  matrix  which,  multiplied  by  the  measurement 
vector,  yields  a  mirror  control  voltage  vector  that 
minimises  the  residual  variance  of  the  wave  phase  after 
correction. 

This  is  done  by  defining  a  matrix  characterising  the 
influence  of  the  mirror  on  the  wavefront  sensor.  This 
matrix,  called  interaction  matrix  and  denoted  D,  is 
obtained  by  applying  a  unit  command  to  each  motor  in 
turn  and  recording  each  measurement  vector  associated 
with  each  motor.  All  the  measurement  vectors  are 
arranged  in  columns  to  form  matrix  D. 

Then  for  any  control  vector  1V>,  it  is  possible  to 
calculate  the  measurement  vector  1e>  by: 

le>  =  D  1V> 

Calculating  IV>  from  1e>  amounts  to  invert  the 
above  equation. 

As  a  general  rule,  interaction  matrix  D  is  rectangular 
(there  are  more  measurements  than  there  are  actuators  to 
be  controlled).  A  criterion  must  therefore  be  chosen  to 
calculate  the  control  matrix.  Since  the  resulting  matrix 
must  minimise  the  residual  errors  after  correction,  the 
least  squares  minimisation  criterion  is  used,  yielding  the 
following  result: 

(‘D  D)  IV>  =  to  lE> 
where  is  the  transposed  matrix  of  D. 

The  control  matrix  is  therefore  computed  by 
inverting  the  matrix  product  (^D  D),  a  symmetric  square 
matrix.  The  singularity  of  this  matrix  is  analysed  by 
determining  its  eigenvectors  and  eigenvalues.  This  yields: 

(‘D  D)  =  P  A  tp 

where: 

•  P  :  column  matrix  of  eigenvectors 

•  A  :  diagonal  matrix  of  eigenvalues. 

An  eigenvector  represents  a  spatial  deformation 
profile  of  the  deformable  mirror.  The  corresponding 
eigenvalue  represents  its  sensitivity  to  the  wavefront 
analysis.  If  one  of  the  eigenvalues  is  zero,  matrix  (^D  D)  is 
not  invertible.  This  mode  with  a  zero  eigenvalue  can  then 
be  filtered  out  to  calculate  the  control  matrix.  The  result 
becomes: 


D*=  PA-1  ‘P‘D 

where 

•  D*  :  optimal  control  matrix 

•  A'l  :  inverse  of  the  eigenvalue  matrix  in  which  the 
inverse  of  the  eigenvalues  corresponding  to  the  modes  to 
be  filtered  out  is  replaced  by  0. 

This  makes  it  possible  to  compute  the  control  matrix 
after  filtering  out  the  modes  that  are  relatively  insensitive 
to  the  wavefront  analysis  to  prevent  them  from  causing 
too  much  noise  on  the  command. 

The  tools  used  to  record  the  interaction  matrix, 
determine  the  eigenvalues  and  eigenvectors,  display  them, 
select  the  modes  to  be  filtered  and  finally  compute  the 
control  matrix  are  available  in  the  computer  function 
libraries. 

This  last  matrix  is  used  by  the  CPU  card  to  generate 
the  deformable  mirror  control  voltages  in  real  time. 

5.2.  Temporal  Correction 

An  Adaptive  Optics  system  is  a  servo  closed-loop 
system.  It  is  governed  by  control  system  equations.  It  is  a 
multivariable  servo  system  which  can  be  separated  into  as 
many  single-variable  feedback  systems  operating  in 
parallel  as  there  are  actuators  on  the  deformable  mirror.  In 
this  case,  if  the  control  variables  are  the  eigenvectors 
described  above,  then  they  become  eigenmodes  of  the 
system  and  they  are  independent  because  of  the 
orthogonality  property  of  these  modes. 

To  ensure  feedback  loop  stability  and  accuracy,  a 
recursive  algorithm  is  executed  by  the  real-time  computer 
CPU  card  to  apply  various  temporal  correctors.  For 
instance,  the  user  can  easily  defined  a  PID  type  controller 
or  a  Smith  controller  or  any  other  type. 

In  the  case  of  SILVA,  since  the  temporal  constraints 
are  not  very  severe,  a  simple  integral  controller  is  used  to 
cancel  the  static  errors  of  the  dye  laser  optical  beam  and 
ensure  a  bandwidth  of  a  few  Hz. 

6.  Performance 

6. 1 .  Simulations 

6.1.1.  Bandwidth 

ONERA  has  been  conducting  theoretical  research  on 
the  control  of  Adaptive  Optics  for  several  years  [6.1].  The 
temporal  behaviour  of  each  in  the  system  components  is 
well  known.  The  open-loop  transfer  function  of  the 
SILVA  Adaptive  Optics  can  therefore  be  completely 
simulated.  In  the  case  of  use  of  simple  integral  action,  it 
can  be  written: 

1  -  e'Ts 

G(s)  =  Ki^e-'^^ 

where: 

s  is  Laplace  variable  (s  =  jto  ;  ©  =  2TCf :  =  -1) 

Kj  is  integral  controller  gain 
T  is  control  period 
T  is  pure  control  delay. 


Frequency  [Hz] 


Figure  6.1.  :  Bode  diagram  of  the  simulated  open-loop 
transfer  function 

For  a  30  Hz  sampling  frequency  and  a  measured  pure 
delay  of  28  ms  due  only  to  reading  of  the  CCD,  the  gain 
Kj  is  set  to  a  value  of  0.65  to  obtain  a  phase  margin  of  45 
degrees.  Therefore,  the  closed-loop  gain  does  not  exceed 
2.3  dB.  Figure  6.1  shows  the  gain  and  phase  of  the 
simulated  open-loop  transfer  function.  The  0  dB  cut-off 
frequency  is  2.7  Hz. 

Experimental  measurements  were  made  on  the 
SILVA  Adaptive  Optics.  The  results  are  given  in  Figure 
6.2,  showing  the  gain  and  phase  of  the  open-loop  transfer 
function  for  an  integral  action  gain  of  0.65.  This  confirms 
the  simulation  results,  with  a  0  dB  cut-off  frequency  of  2.7 
Hz  and  a  phase  margin  of  45  degrees. 
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Figure  6.2.  ;  Bode  diagram  of  the  experimental  open- 
loop  transfer  function 


6.1.2.  Correction  Accuracy 

ONERA  has  a  range  of  numerical  simulation 
software  to  model  the  behaviour  of  each  of  the 
components  of  an  Adaptive  Optics  system  and  the  phase 
errors  to  be  corrected.  The  deformation  profile  of  each  of 
the  actuators  is  taken  into  account  as  well  as  operation  of 
the  Hartmann-Shack  analyser  and  its  measurement 
accuracy  according  to  the  number  of  subpupils,  the  flux 
and  the  read  noise.  The  control  matrix  is  computed  from 
simulation  of  the  optical  interaction  matrix.  The  modal 


filtering  methods  described  above  are  taken  into  account. 
Finally,  the  loop  transfer  function  is  represented. 


Figure  6.3-a  :  Defocusing  before  correction 


Figure  6.3-b  :  Defocusing  after  correction 

From  this  set  of  software,  the  SILVA  Adaptive 
Optics  system  correction  accuracy  was  modelled  for  the 
main  errors  present  on  the  dye  laser  beam.  Two 
parameters  were  calculated: 

•  The  ratio  between  the  peak-to-peak  error  before 
correction  and  the  peak-to-peak  error  after  correction, 
denoted  Gcc 

•  The  ratio  between  the  peak-to-peak  error  before 
correction  and  the  standard  deviation  after  correction, 
denoted  Gcc/rms 

For  defocusing,  we  found  Gcc  =  9.3  and  Gcc/rms  = 
66.  For  astigmatism,  we  obtained  Gcc  =  8  and  Gcc/rms  = 
60.  The  correction  accuracy  achieved  is  mainly  limited  by 
the  number  of  actuators  of  the  deformable  mirror  and  the 
deformation  profile  of  each  actuator.  Figures  6.3-a  and 
6.3-b  show  the  3D  defocusing  before  and  after  correction. 
The  two  graphs  are  to  the  same  scale.  For  instance,  for 
defocusing  exhibiting  a  peak-to-peak  error  of  X,  a  residual 
error  of  A/66  is  obtained  after  correction. 


Q 


6.1.3.  Propagation  Simulation 

Figures  6.4-a  and  6.4-b  show  the  illumination  and 
typical  wavefront  profile  of  a  laser  beam  from  the  SILVA 
system.  The  measurements  were  made  using  a  CCD 
camera  for  illumination  and  a  Hartmann-Shack  wavefront 
sensor  for  the  optical  phase.  Figure  6.4-a  shows  the  edge 
effects  visible  in  one  dimension  of  the  illumination 
profile,  characteristic  of  beam  guidance  by  amplifier 
ducts. 
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Figure  6.4-a  :  Measured  dye  laser  beam  illumination 
profiles 


Figure  6.4-b  :  Measured  dye  laser  beam  wavefront 
profiles 

Figures  6.5-a,  6.5-b  and  6.5-c  illustrate  the 
improvement  expected  in  the  beam  shape  during 
propagation.  These  profiles  result  from  numerical 
simulation  of  propagation  by  a  double  Fourier  transform 
computation  code.  The  initial  intensity  profile  is  the 
same,  and  only  the  wavefront  considered  has  changed: 

•  Figure  6.5-a  :  The  two  main  curvatures  (mean 
curvature  and  astigmatism)  of  the  initial  wavefront  are 
corrected.  This  case  is  representative  of  the  results 
that  can  be  achieved  in  the  absence  of  Adaptive 
Optics.  Correction  is  obtained  using  spherical  and 
cylindrical  lenses. 

•  Figure  6.5-b  :  The  first  15  Legendre  modes  of  the 
initial  wavefront  are  corrected  using  Adaptive  Optics. 
It  is  recalled  that  since  the  pupil  geometry  is 
rectangular,  the  optical  phase  can  be  decomposed  on 
the  orthonormal  base  of  the  Legendre  polynomials. 

•  Figure  6.5-c  :  The  wavefront  is  perfectly  flat. 

The  light  intensity  distribution  of  the  beam  on  either 
side  of  the  pupil  image  is  computed  to  analyse  its 
uniformity  during  propagation  in  the  uranium  vapour 
separator.  The  propagation  distance  considered  in  the 
simulation  is  equal  to  (a^/4X,),  where  a  is  the  beam  radius 
and  A,  is  the  wavelength.  It  is  checked  that  better 
uniformity  of  the  beam  illumination  distribution  during 
propagation  is  preserved  after  total  or  partial  correction  of 
the  wavefront. 


Figure  6.5-a  :  Simulated  illumination  profile  without 
Adaptive  Optics:  defocusing  and  astigmatism  corrected. 


Figure  6.5-b  :  Simulated  illumination  profile  with 
Adaptive  Optics:  first  15  Legendre  polynomials  corrected. 


Figure  6.5-c  :  Simulated  illumination  profile  for  a  flat 
wavefront. 

6.2.  Experimental  Characterisations 
6.2. 1 .  Correction  Accuracy 

The  purpose  is  to  measure  the  maximum  correction 
accuracy  of  the  Adaptive  Optics  system  developed.  For 
this  purpose,  the  experiment  consists  of  correcting  the 
intrinsic  flatness  errors  of  the  deformable  mirror  and 
evaluating  the  resulting  surface  quality.  The  experimental 
setup  used  is  described  in  Figure  6.6. 

We  first  used  a  reference  flat  mirror  whose  flatness 
evaluated  by  the  Zygo  was  equal  to  X./16  peak-to-peak. 
The  measurement  supplied  by  the  wavefront  sensor  were 
recorded  to  be  used  as  reference  during  the  tests.  The 
reference  flat  mirror  was  then  replaced  by  the  CILAS 
deformable  mirror  of  which  only  an  area  of  24  x  36  mm 
was  used  (5x4  actuators).  The  feedback  loop  was  then 
closed  and  the  surface  condition  of  the  deformable  mirror 
after  correction  was  evaluated  by  the  Zygo. 


Figure  6.6  :  Experimental  setup  for  evaluating  the 
correction  accuracy  of  the  Adaptive  Optics  system 

Figure  6.7-a  shows  the  surface  condition  of  the 
deformable  mirror  at  rest.  It  exhibits  a  peak-to-peak  error 
of  XU  (for  an  rms  error  of  Ay  10).  Figure  6.7-b  shows  the 
surface  condition  of  the  deformable  mirror  after 
correction.  It  exhibits  a  peak-to-peak  error  of  A./1 1  (rms 
error  of  A/46).  The  correction  accuracy  achieved  is 
excellent  and  consistent  with  the  stated  objectives. 


Figure  6.7-a  :  Surface  condition  of  the  deformable 
mirror  at  rest 


Figure  6.7-b  :  Surface  condition  of  the  deformable 
mirror  after  correction 

6.2.2.  Reshaping  of  the  Laser  Beam 

The  purpose  of  this  experiment,  described  in  Figure 
6.8,  is  to  correct  the  phase  of  an  He-Ne  laser  beam  (beam 
1)  to  make  it  identical  to  another,  so-called  reference 
beam  (beam  2). 


Deformable 


CHANNEL  A  CHANNEL  B 

Figure  6.8  :  Experimental  laser  beam  shaping  setup 


The  phases  of  the  two  beams  are  compared  using  the 
corresponding  far-field  images  recorded  by  the  CCD 
camera  of  the  setup  described.  With  beam  1  masked,  the 
wavefront  sensor  measures  the  phase  errors  of  beam  2 
which  are  then  used  as  reference  in  the  rest  of  the 
experiment.  The  image  supplied  by  the  CCD  camera  is 
recorded.  The  reference  beam  is  then  masked  to  use  only 
beam  1.  The  feedback  loop  is  closed.  The  wavefront 
sensor  (WFS)  measures  the  phase  error  between  the  two 
beams  and  the  real-time  computer  evaluates  the  voltages 
to  be  applied  to  the  deformable  mirror  to  give  beam  1  the 
same  phase  profile  as  beam  2.  The  image  supplied  by  the 
CCD  camera  is  recorded  again. 

Figure  6.9  shows  the  far-field  images  obtained  with 
each  of  the  two  beams.  The  recorded  image  of  the 
reference  beam  is  on  the  left  and  the  recorded  image  of 
beam  1  reshaped  by  the  Adaptive  Optics  system  is  on  the 
right.  Comparison  of  these  two  images  shows  the 
operational  quality  of  the  system  developed. 


Figure  6.9  :  Far-field  image:  beam  2  on  the  left,  beam 
1  on  the  right 

7.  Conclusion 

An  adaptive  optics  system  was  developed  for  the 
SILVA  prototype.  Its  development  benefited  from 
research  conducted  by  ONERA  for  several  years  on  large 
bandwidth  Adaptive  Optics  dedicated  to  space 
observation. 

The  system  consists  of  a  deformable  mirror  with  48 
actuators  and  a  wavefront  sensor  with  21x21  subpupils. 
It  corrects  the  wavefront  errors  present  on  the  dye  laser 
systems.  The  control  frequency  of  the  system  is  limited  to 
30  Hz.  The  hardware  and  software  computer  architecture 
designed  is  flexible  and  reconfigurable,  which  is  very 
desirable  for  a  prototype  during  the  testing  and 
optimisation  phase.  In  particular,  it  will  allow  a  control 
frequency  of  500  Hz  to  be  achieved  simply  by  using  a 
64x64-pixel  CCD  camera. 

This  system  was  integrated  and  tested  in  the 
laboratory  at  ONERA.  Its  performance  is  consistent  with 
the  simulations  made.  It  is  now  installed  at  CEA  where  it 
is  operating  on  one  of  the  dye  laser  beams  of  the  SILVA 
prototype. 
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ABSTRACT 

In  inertial  confinement  fusion,  different  shapes  and  intensity  distributions  of  focal  spots  are  desired  for  different 
physical  experiments.  In  this  paper,  the  two  methods  of  phase  retrieval  (Geometric  Transformation  and  Iterative 
Fourier  Transformation)  in  controlling  laser's  focal  intensity  profile  are  been  introduced.  Deformable  mirrors  with 
different  actuator  numbers  and  arrangements  (triangle  and  rectangle)  are  simulated  to  compare  fitting  capabilities  for 
the  given  wave-front  profile.  The  possibilities  of  using  Adaptive  Optical  technology  with  these  two  methods  to  control 
laser's  focal  intensity  profile  are  been  compared  and  analyzed  by  computer  simulations. 

Keywords:  beam  shaping.  Adaptive  Optics,  deformable  mirror,  computer  simulation 

1.  INTRODUCTION 

There  are  some  applications  in  laser  technology  to  form  some  kinds  of  focal  intensity  profiles  or  intensity  distributions. 
Especially  in  inertial  confinement  fusion  (ICF)  system,  a  high  power  laser  pulse  is  focused  on  a  target  to  induce  nuclear  fusion. 
Different  shapes  and  intensity  distributions  of  focus  spots  are  desired  for  different  physical  experiments.  Some  demands  are 
critical  such  as  uniformity  of  light  intensity  distribution  and  sharp  edge  of  the  focus  spot,  very  low  background  intensity  level 
outside  the  spot.  In  recent  years,  there  are  some  methods  for  laser  beam  shaping.  Reflective'*^  and  refractive'^^  surfaces  have 
been  used  to  redistribute  the  energy  of  a  laser  beam  to  a  uniform  beam  with  circular  shape.  These  methods  generally  require 
complicated  aspherical  surfaces  that  are  difficult  to  fabricate  and  are  useful  only  in  a  limited  set  of  situations.  Some  binary  and 
diffractive  optical  components''*"’"®'  are  also  used  to  form  the  beam  shaping  system.  These  components  have  the  adequate 
freedoms  and  precision  to  form  various  types  of  phase  distribution.  The  drawback  of  them  is  that  they  can  not  tolerate  for  daily 
changing  of  original  phase  of  laser  beam.  The  original  phase  is  often  changing  that  may  influence  the  desired  intensity 
distributions  at  the  focus  spot.  By  now,  another  effective  method  for  getting  desired  focal  intensity  profiles  is  to  obtain  the  phase 
modulation  by  inserting  deformable  mirror  (DM)  in  incident  beam. 

DM  is  a  reflective  component  that  it  has  the  better  performances  of  transmissive  efficiency,  destroying  threshold  and  the 
adaptability  for  different  environments  than  the  other  components.  Because  the  DM  can  adjust  its  surface  shape  at  any  time  to 
adapt  the  wave-front  changed,  it  tolerates  to  different  demands  of  experiments.  In  this  paper,  two  methods  of  phase  retrieval 
{Geometric  transformation  and  Iterative  Fourier  Transformation)  of  using  DM  to  retrieve  beam  phase  at  pupil  plane  for 
controlling  focal  intensity  profile  will  be  described.  The  fitting  capabilities  of  DM  with  different  actuator  arrangements  are  also 
analyzed  by  computer  simulations.  The  possibilities  of  using  Adaptive  Optical  technology  to  control  laser's  focal  intensity  profile 
are  discussed. 


2.  THEORY 

The  relation  between  complex  amplitude  a(x,  y)  at  mirror  and  complex  amplitude  A(u,v)  at  focal  spot  can  be  represented  by 
far-field  diffractive  equation'*®' : 


A(u,v): 


iV 


(1) 


Here,  (x,y)  and  (u,v)  are  the  coordinates  at  the  mirror  and  the  focal  plane,  X  is  wave  length,  f  is  the  focus  length  of  lens,  X  is 
pupil  area  and  k=27c/X. 


The  main  ideal  to  reshape  the  intensity  distribution  of  focal  spot  is  to  insert  a  phase  function  <l>{x,y)  into  a(x,  y),  and  the 
expected  intensity  distribution  at  focus  spot  is  gotten.  The  principle  is  expressed  as  fig.  1. 


Incident  beam  enters  the  optical  system  and  the  phase  distribution  of 
wave-front  is  detected  S-H  sensor.  At  the  same  time,  the  focal  intensity 
profile  is  detected  by  an  intensity  detecting  sensor  (CCD  camera). 
Comparing  the  data  by  CCD  with  that  by  S-H  sensor,  we  can  control  the 
DM  to  insert  the  desired  phase  distribution  into  the  beam.  After  inserting 
the  additional  phase,  the  wave-fi'ont  at  focus  spot  becomes: 


A(u,v) 


i/lz 


E 

(j){x,y)  is  the  phase  produced  by  DM,  <l>^{x,y)  is  the  original  phase  of 
incident  beam,  and  (l>^,{x,y)  =  <l>i{x,y)  +  <l^x,y). 


The  intensity  distribution  on  focal  plane  is: 

l{u,v)  =  \A{u,vf 


(3) 


From  descriptions  above,  the  focal  intensity  profile  is  determined  by  phase  distribution  resulted  from  the  surface 

shape  of  DM.  (I){x,y)  can  be  solved  by  following  two  methods  through  the  relation  between  intensity  distribution  at  pupil  plane 
and  that  at  focal  plane. 


2.1.  Geometric  transformation  method  (GTM) 


Approximating  the  far-field  integral  in  eq.(2)  by  stationary  phase  method  originated  by  Olof  Bryngdahl^”',  it  can  be  shown 
that  a  point  (x,y)  in  the  pupil  plane  is  mapped  to  a  point  (u,v)  in  the  focal  plane  given  by  the  relation 


'  f  Ik- d^^{,x,y^l  ax 

UJ 

IvJ 

flk-apAx,y)l^_ 

If  the  desired  mapping  is  defined  by 


'x^ 

’v(x,y) 

A. 

I'j 

y(x,y)_ 

The  necessary  phase  {x,y)  to  produce  this  mapping  is  the  solution  of  the  partial  differential  equations  as  eq.(5): 


(4) 


& 


f 


£ 

/ 


■v{x,y) 


(5) 


(p{x,y)  can  be  obtained  by  p^ix,y')- p^(x,y)  if  Pi{x,y)  is  measured  by  S-H  sensor  and  the  relationship  between  point  (u,v) 
and  (x,y)  is  known  .  The  DM  surface  shape  S(x,y)  is  given  as 

S(x,y)=P{x,y)/{k/2)  =  [p^(x,y)-P'ix,y)]/(k/2)  (6) 

In  our  application,  incident  beam  truncated  by  rectangular  window  is  expected  to  be  reshaped  to  the  uniform  focus  spot. 
Assume  the  size  of  truncate  window  at  pupil  plane  is  (±x„,±>’,)  and  the  size  of  focus  spot  is  (±Mj,±Vj).  The  relationship 
between  point  (u,v)  and  (x,y)  can  be  solved  based  on  energy  conservation  principle  given  by  Parseval's  theorem,  the  ener^^  of 
pupil  plane  is  equal  to  that  of  focal  spot: 


nyo 

I(x,y)dxdy  =  (Mo,v,)  =  •2v„ 

-ya 


Here  a  is  uniform  intensity  distribution  at  focal  plane.  The  relation  between  {x,y)  and  (K,v)can  be  obtained  similarly  as  point 
(^o>>'o)  8®*  following  equation: 


1  f'o  r-^« 

a=- -  \l{x,y)dxdy 

4U(|V,  J-xoJ-ya 


uv  =  u„v„ 


r-'o  fyo 

Iix,y)dxdy 

J-Jifi  J-vn 


(7) 


u .  V  are  the  function  of  x  and  y.  In  this  transformation  system  by  using  this  method,  the  compatibility  condition^^'  must  be 
satisfied  by  the  input  and  output  beams. 


Mx,y)  ^  di'(x,y) 

A 

The  eq.(8)  expresses  is  that  ^^(x,y)  must  simultaneously  satisfy  both  equations  in  eq.(5).  From  above  list,  the  desired  focal 
intensity  profiles  are  obtained  if  the  DM  can  produce  the  necessary  surface  profiles  determined  ly  eq.  (6).  (7)  and  (8). 


2.2.  Iterative  Fourier  Transformation  method  (IFTM) 

The  relation  between  pupil  plane  and  focal  plane  can  be  expressed  by  Fourier  Transformation  (FT)  because  of  the  existence  of 
focusing  lens.  The  complex  amplitude  at  pupil  plane  is  transformed  to  that  at  focal  plane  by  FT.  The  desired  intensity 
distribution  will  be  obtained  by  inserting  a  phase  function  that  can  be  gotten  by  IFTM  into  original  phase.  IFTM  is  also  called 
Gerchberg-Saxton  algorithm^'^^^  which  was  originally  invented  in  connection  with  the  problem  of  phase  retrieval  from  two 
intensity  measurements  (pupil  and  focal  plane),  and  had  been  modified  by  J.R.Fienup*®^  and  Guozhen  Yang  etc.^^”*'  in  some 
applications  successfully.  The  main  principle  of  IFTM  is  as  following  and  shown  as  fig.2. 

©Constructing  original  complex  amplitude  a^{x,y)  by  intensity  I(x,y)  and  phase  function  (j>^{x,y)  which  can  be  formed  by 
random  number; 

Oo  ix,y)  =  ■Jl(x,y)  ■  exp[iV„  (x,y)] 

This  algorithm  in  k-th  iterative  consists  of  the  following  four  steps  @  ~  ® : 


Fig.2 ,  Principle  of  Iterative  Fourier  Transformation  method 


©Fourier  transform  the  complex  amplitude  of  pupil  plane  to  that  of  focal  plane: 

(w,v)  =  1^  (i/,v)|exp[;'F*  (m,v)]  =  FT{|o4  (x,>')|exp[/^i  (jc,>')]} 

©Replace  the  modules  of  the  computed  results  |>l^(w,v)|  with  the  expected  Fourier  modules  ^I^(u,v)  to  form  a  new 
amplitude  at  focal  plane: 


(w.v)  =  expr/^*  (m.v)] 

©Inverse  Fourier  transform  the  complex  amplitude  at  focal  plane  to  that  at  pupil  plane: 


a'k  (x,y)  =  \al  (x,>/)|  (x,y)]  =  Fr^[Al  (k,v)]  =  Fr^yj„(u,v)  exp[iT*  (k,v)]] 

©Replace  the  modules  of  the  computed  results  with  the  measured  ones  y]l{x,y)  to  form  a  new  amplitude  of  the 

pupil  plane: 

ak+i(x,y)  =  yll(x,y)-e\p[i<^[  (x,^)]  =  Jl(x,y)  -expl^^^ix^y)] 

©Return  to  step  @  until  mean-square  error  of  |ai  to  |n*+,(x,>/)| 


satisfies  the  setting  constraint.  The  expected  phase  at  pupil  plane  and  intensity  at  focus  spot  will  be  output. 


(9) 


This  algorithm  is  used  for  any  problems  in  which  particular  constraints  are  known  in  both  pupil  and  focal  plane.  The 
iterations  continue  until  the  computed  Fourier  transform  satisfies  the  Fourier-domain  constraints  or  its  inverse  transform  satisfies 
the  pupil  plane  constraints,  then  one  has  found  a  solution  that  a  Fourier  transform  pair  satisfies  all  the  constraints  in  both 
domains.  We  use  three  parameters  to  represent  the  intensity  quality  of  focal  spot  quantitatively: 


©The  mean  square  error  of  simulated  amplitude  (MSEuv)  to  desired  amplitude  at  focal  spot. 


Z  z  ,V„  )|  -  (w„ , v„ )  ] 


MSEuv  = 


M~l  N-l 


©Light  energy  efficiency  (LEE)  defined  as  the  ratio  of  energy  concentrated  in  effective  dimension  at  focal  spot  to  all  energy 
of  incident  beam. 


LEE  = 

/ 

j‘jl(x,y)-dx-dy 

/ 

_-ao->0 

X  100% 


minimum  to  peak  value. 


3.  COMPUTER  SIMULATION 

There  is  a  requirement  in  a  physics  experiment  that  the  intensity  of  incidence  beam  with  six-order  super-Gauss  rectangular 
profile  which  is  truncated  by  rectangular  window  be  formed  rectangular  uniform  intensity  at  the  focus  spot  by  the  optical  system 
with  parameters  as  table  1. 


Table  1 .  System  parameters 


1 

wave  length 

1.054pm 

2 

laser  aperture 

220mm  x  220mm  (rectangle) 

3 

size  of  focus  spot 

300pm  X  300pm  (rectangle) 

4 

original  energy  distribution 

I  {x,y)  =  exp[-2(x  /  x„ )‘  ]  •  exp[-2(>'  )‘  ] 

5 

focus  lens 

fiT)=7 

The  intensity  distribution  at  pupil  plane  can  be  expressed: 

I(x,y)  =  (10) 

Here,  x„  and  Vq  are  the  Gaussian  half  widths  and  also  the  half  widths  of  truncate  window.  Assume  .;C||=yj=150mm,  the  beam 
intensity  distribution  is  as  fig.3: 


fctoisiy 


•10S  «  •  M 


(a)2-D  graph 


;<(Qxn) 

(b)l-D  graph 


Fig.3  >  Beam's  origin  intensity  distribution  at  mirror 
3.1.  Geometric  transformation  method  (GTM) 

Assume  incident  beam  is  planar  wave-front.  Joining  eq.(5),  (7).  (8)  and  (10) ,  the  desired  phase  function  will  be  solved  as: 

<Kx,y)  =  J  \'piyi)dyi 

p(yo)->o 


and  p(/)  =  J 

The  computed  desired  phase  distribution  is  as  fig.4. 


phase 


(a)2-D  graph 

Fig.4 ,  Computed  value  of  phase  distribution 

After  removing  defocus,  the  actual  needed  phase  distribution  is  as  fig.  5. 


(a)2-D  graph 

Fig.5,  Phase  distribution  after  removing  defocus 

The  diffractive  intensity  distribution  is:  I(u,v)  =  — ^-y[j?^(w,v)+ J?^(u,v)] 

(Az) 


and 


=  JJ|a(^,>')|-cos 


<l>ix,y)--j{vx  +  vy) 


uxdy 


sin 


^x,y)—j:{ux  +  vy)  ^dy 


The  computed  intensity  distribution  at  focus  spot  is  as  fig.  6. 


Iitensky 


(a)2-D  graph 


From  fig.  6,  there  are  oscillations  on  the  top  of  intensity  profile.  These  are  caused  by  finite  truncate  window  at  pupil  plane.  By 
our  simulations,  the  smaller  the  truncation  window  is,  the  more  oscillations  the  intensity  distribution  is.  When  the  truncate 
window  is  enlarged  twice  (600mm  x  600mm),  the  oscillations  are  fewer  and  the  center  is  flatter  as  fig.  7. 


(a)2-D  graph  (b)l-D  graph 

Fig.7,  Intensity  distribution  at  focal  spot  with  truncate  window  600mm  x  600mm 

We  can  get  the  value  of  MSEuv,  LEE,  Grad  as  table  2  expressed  from  fig.6.  The  shortage  of  this  method  is  the  large 
amount  of  computation  if  the  higher  computing  precision  to  be  satisfied. 

3.2.  Iterative  Fourier  Transformation  method  (IFTM) 


The  convergence  of  this  method  had  been  proved  by  J.R.Fienup  that  the  algorithm  decreases  the  mean-square  error  MSExy  (or 
stays  the  same)  in  eq.  (9)  at  each  iteration.  In  practice,  this  algorithm  usually  decreases  the  MSExy  rapidly  for  the  first  few 
iterations  but  much  more  slowly  for  later  iterations.  The  speed  of  convergence  also  depends  on  the  type  of  constraints  imposed^'®^ . 
It  has  been  found  that  the  computed  convergent  phase  is  strongly  depending  on  the  different  original  phase  distributions.  Fig.8 
shows  the  one-dimensional  intensity  distribution  after  1200  iterations  with  two  original  phase  distributions  which  are  random  in 
ranges  of  (-7r,7i)  and  (-7i/10,7i/10)  respectively. 


lateusity 


u(mm) 


Intensity 


u(mm) 


(a)Initial  phase  area  (-n,n) 


(b)Initial  phase  area  (-7r/10,7t/10) 


Fig.8,  Intensity  distribution  with  1200  iterations  (ID) 


After  analyzing  for  the  computing  results,  we  find  that  the  more  approximate  the  original  phase  is  to  the  theory  value,  the 
more  exact  the  resulting  intensity  distribution  is  to  the  desired  ones  and  the  iterative  times  are  relatively  fewer.  For  this  reason, 
the  computed  phase  distribution  by  GTAf  may  be  used  as  the  original  phase  in  this  algorithm.  Fig.9  shows  the  computed  results 
by  this  method  with  120  iterations. 


(a)3-D  graph  (b)2-D  graph 

Fig.9%  Intensity  distribution  with  120  iterative 


From  fig.9,  MSEuv^  LEE  and  Grad  by /FTA/are  gotten  and  compared  with  those  by  GTM  from  fig.6  as  table  2. 


Table  2  ^  The  performance  parameters  of  two  methods 


method 

MSEuv 

LEE 

Grad  (pm) 

GTM 

0.030162 

93.82% 

52.8 

IFTM 

0.018838 

99.19% 

7.42 

By  using  IFTM,  the  edge  of  the  intensity  distribution  is  much  sharper,  the  center  of  the  intensity  is  flatter,  and  the  lateral  lobes 
result  from  diffraction  is  smaller  than  those  by  using  GTM.  This  method  has  performances  of  fast  computing  speed  by  using^oy/ 
Fourier  transformation,  favorable  stability  and  convergence. 

4.  FITTING  CAPABILITY  OF  DEFORMABLE  MIRROR 

DM  is  the  key  component  of  this  beam  shaping  system.  It  changes  the  beam  phase  when  the  actuators  of  DM  move  to  change 
the  mirror  surface.  Its  surface  can  be  adjusted  arbitrarily  and  maintained  for  a  long  time.  The  fitting  capabilities  of  DM  to 
expected  surface  shape  with  different  acttiator  arrangements  and  actuator  numbers  have  been  investigated.  The  four  probable 
arrangements  as  Fig.  1 0  are  simulated  when  the  incident  angle  is  45°. 


Fig,  10  >  Several  actuator  arrangements  with  different  number 

Deformations  of  DM  caused  by  all  actuators  form  the  complicated  surface  profile.  It  is  very  important  that  the  obtained  profile 
approximates  to  the  desired  mirror  surface  as  well  as  possible.  The  fitting  capability  of  DM  with  above  four  actuator 
arrangements  is  investigated  as  following. 

Assume  the  influence  function  /(/-)  can  be  represented  as  Gaussian  shape:  f  (r)  =  exp[-(|F |  /  rOj )“ ] 

Here,  a  is  Gaussian  index,  co^  is  half  width,  r  is  the  distance  apart  from  driving  actuator.  The  values  of  and  a  depend  on  the 
actuator  space  ^  mirror  stiffness  ^  connecting  structure  between  actuators  and  mirror  etc.  The  following  parameters  are  used  in 
the  simulations:  a=2.2  and  cof-9.  165mm.  The  resulted  wave-front  can  be  represented  as: 


/  I  —  I 

=  -^J  ■ 

'-‘  (  L  ^  J 

Here,  F  is  an  arbitrary  point  on  the  mirror,  N  is  actuator  number,  r.  is  the  coordinate  of  the  i-th  actuator,  and  Cj  is  the  driving 
signal  of  the  i-th  actuator,  represented  in  matrix; 

I,  = 

This  matrix  can  be  solved  by  SVD,  and  the  least  square  estimate  Cj  of  Cj  can  be  obtained: 

q=(C„Q,-,C^/ 

For  a  point  (x,  ,y^)  on  the  mirror,  there  is  a  corresponding  A^j  that  the  wave-front  is  resulted: 

g’^AljxC’j 

So,  the  fitting  error  may  be  represented  as; 

rnts  = 

The  intensity  distribution  at  focal  spot  and  the  actual  wave-fronts  by  adding  the  fitting  errors  of  the  different  DM  are  showed 
as  fig.  11,  the  computed  results  are  showed  as  table  3. 


(i)3-D  intensity  (ii)2-D  intensity  (iii)actual  wave-front 

(a).  35-actuator 


(i)3-D  intensity  (ii)2-D  intensity  (iii)actual  wave-front 

(b).  40-actuator 


(i)3-D  intensity 


ktfmsily 


(ii)2-D  intensity 
(c).  48-actuator 


(d).  60-actuator 


(iii)actual  wave-front 


Fig.  1 1 ,  Intensity  distribution  of  four  arrangements  for  fitting  wave-front 


Table  3 ,  The  capability  of  four  arrangements  for  active  nurror  to  fit  wave-front 


NO. 

actuator  number 

fitting  error:  rms  (A.) 

1 

35 

0.021293 

2 

48 

0.012763 

3 

40 

0.018575 

4 

60 

0.015813 

From  the  simulated  results  as  fig.  1 1,  MSEuv.  LEE  and  Grad  are  as  table  4. 


Table  4.  The  performance  parameters  of  four  arrangements  by  simulations 


actuator  number 

arrangement 

MSEuv 

LEE 

Grad  (pm) 

35 

rectangle 

0.071456 

91.09% 

61.63 

40 

triangle 

0.095204 

91.18% 

61.63 

48 

0.062939 

91.21% 

61.63 

60 

triangle 

0.065385 

91.23% 

61.63 

Some  conclusions  can  be  obtained  from  the  simulative  results; 

■  ®The  DM  with  48-actuator  and  rectangular  arrangement  has  the  smallest  fitting  error  for  desired  wave-front  in  four 
arrangements. 

■  ©The  wave-front  fitted  by  DM  with  triangular  actuator  arrangement  is  asymmetry  because  of  the  asymmetry  of  actuators. 

■  ©Although  the  fitting  error  from  60-actuator  with  triangular  arrangement  is  as  small  as  48-actuator,  the  many  more  actuators 
could  increase  the  complexity  of  system. 

■  5.  CONCLUSION 

I  ■  ' 

■  Computer  simulations  have  been  done  to  demonstrate  the  possibility  of  using  DM  with  adaptive  technology  for  beam  shaping 

I  ■  in  ICF  system.  To  obtain  a  high  energy  efBciency,  intensity  distribution  with  uniform  profile  and  sharp  edge  at  the  focus  spot, 
the  beam  energy  must  be  redistributed  by  inserting  the  necessary  phase  function. 

I 

I  Two  methods  of  phase  retrieval  have  been  investigated  to  compute  the  necessary  phase  distribution  for  this  application.  It  has 

j  been  found  that  GTM  based  oh  using  a  stationary  phase  approximation  to  the  far-field  integral  is  inadequate  to  get  the  desired 

I  _  intensity  distribution  because  of  the  top  modulation  and  the  inadequately  sharp  edge  of  the  focal  intensity  profile.  A  better  result 

■  has  been  obtained  by  using  IFTM. 

Calculating  equations  based  on  GTM  have  the  direct  expresses  and  help  us  understand  the  relations  of  focal  intensity  profiles 

■  with  beam  and  pupil  parameters.  IFTM  has  faster  computing  speed  by  using  Fast  Fourier  Transformation.  It  will  be  performed 

I  when  intensity  information  and  some  restraints  (pupil  and  focal  domains)  are  known.  From  the  simulation,  it  is  obvious  that 

IFTM  is  more  applicable  than  GTM. 


DM  with  different  actuator  numbers  and  arrangements  (triangle  and  rectangle)  are  investigated  to  demonstrate  its  fitting 
capabilities  for  the  given  wave-front  profiles.  Adding  the  fitting  error  produced  by  DM,  the  focal  intensity  profiles  become  worse 
than  the  computed  result  of  phase  retrieval.  DM  with  48  actuators  and  rectangular  arrangement  has  the  best  fitting  capability  for 
known  wave-front  in  four  actuator  arrangements. 
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Adaptive  Systems  with  Modal  Aperture  Probing 
for  Powerful  Technological  IR  Lasers. 


D.Yu.  Gorbenko,  Ya.I.  Malashko,  A.S.  Rumyantsev,  V.V.  Valuev 

“Sistema”  SIC,  Moscow,  Russia,  Leningradsky  prospekt,  80.  Fax:  (095)  158-5671 


The  characteristics  of  the  adaptive  systems  with  sequential  and  parallel  multifrequencies 
modal  probing  were  studies. 

Locking  of  an  adaptive  optical  system  sequential  modal  probing  with  the  CO2  gas- 
dynamic  laser  emission  of  80  kW  output  power  was  fulfilled.  The  AOS  task  was  to  maximize 
angular  on-axis  brightness  which  was  realized  by  analyzing  of  diffractionally  coupled  beam  in  the 
equivalent  far  field.  Optical  scheme  contained  a  cooled  flexible  correcting  mirror  with  19 
hexagonally  located  piezoactuators  with  the  steps  of  33  mm  and  local  correcting  amplitudes  up  to 
10  )im  (1,2].  Correcting  mirror  was  used  also  as  a  modulating  mirror. 

To  increase  AOS  time-band  perfection  up  to  1-^2  Hz  although  infringing  the  spatial  band 
the  correcting  mirror  operated  in  so-called  modal  mode  forming  the  simplest  specified  surfaces 

during  probing  and  correction:  two  tilts,  two  cylinders  Z^X^,  Z^y^  and  a  “saddle”  Z~xy.  The 
given  functions  reflected  low-frequency  spatial  aberrations  in  a  beam. 

Resonance  frequencies  of  correcting  mirror  structure  and  laser  power  fluctuation 
spectrum  were  investigated  and  their  influence  on  the  stability  of  AOS  operation. 

The  resulting  success  of  the  experiment  was  conditioned  by  the  choice  of  such  probing 
frequency  (1,1  kHz)  which  coincides  neither  with  the  own  resonances  of  mirror  nor  with  the 
spectral  bands  of  laser  power  fluctuations.  The  increase  of  average  angular  beam  brightness  was 
two  times. 

To  increase  AOS  time-band  pass  a  multifrequency  parallel  probing  was  applied  [3].' 

Although  the  described  linear  system  doesn’t  have  wide  time^and-^atial  bands,  ii  has  a 
wide  range  of  correction  whi^h  is  very  important  for  powerful  laser  system  with  low-frequency 
phase  perturbations  distortion  of  large  amplitudes. 
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Introduction 


•  The  Optical  Sciences  Company  has  developed  a  low- 
cost  wavefront  sensor  (WFS)  and  beam  profiler 

•  The  sensor  has  proved  to  be  an  important  tool  in  aiding 
optical  system  alignment  of  high  power  lasers 

•  The  sensor  is  capable  of  estimating  both  intensity  and 
phase  profiles  in  near  real  time  (currently  5  Hz)  -v 

•  The  sensor  is  capable  of  measuring  pulsed  mode 
wavefronts 


•  In  this  presentation  we  describe 


-  the  workings  of  the  sensor  hardware, 

-  features  of  the  graphical  user  interface,  and 

-  how  the  WFS  can  be  applied  to  the  measurement  of 

high  power  lasers  t 
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The  Wavefront  Sensor 


•  The  wavefront  sensor  uses  the  concept  of  the 
Hartmann  test 

•  A  two  dimensional  array  of  lenses  is  used  to 
divide  the  beam  into  subapertures 

•  Each  individual  lens  forms  an  independent 
focus  on  a  CCD  array 

•  Each  spot  location  is  an  estimate  of  the 
average  phase  gradient  (or  tilt)  over  the 
corresponding  subaperture 

•  The  sensor  works  equally  well  for  white-light 
sources 
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WFS  Description 


•  A  frame  grabber  board  records  the  CCD  focal 
plane  data  to  a  PC  Pentium  computer 

•  The  computer  calculates  spot  locations  to  the 
millipixel  level,  using  either  matched  filter  or 
centroid  techniques 

•  A  “perfect”  plane  wave  is  used  to  calibrate  the 
sensor 

•  The  spot  locations  for  the  plane  wave  are 
stored  and  used  as  reference  locations 

•  Tilts  are  introduced  to  calibrate  the  lenslet 
focal  lengths 


WFS  Optical  Block  Diagram 


•  Shown  here  is  the  simplest  configuration  of 
the  wavefront  sensor 
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I  Example  Aberrations 

I  Focus  aberrations 


LENSLET 

APPAY 


rOCUSED 

SPOTS 


LENSLET 

ARRAY 


racusEB 

SPOTS 


FOCUSED 


Diverging  F/6  Beam 
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Measuring  a  Test  Wavefront 


•  New  spot  locations  are  found  for  the  distorted 
wavefront  under  test 

•  The  new  spot  locations  are  differenced  from 
the  previously  calibrated  and  saved  “nominal” 
spot  positions 

•  These  differences  are  proportional  to  the  test 
wavefront’s  phase  differences  across  the 
subapertures 

•  Custom  software  is  used  to  reconstruct  the 
wavefront  phase  from  this  two  dimensional 
set  of  phase  differences 

u 
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Thresholded  Centroid  Subaperture  Spot  Finding 


•  Within  each  region,  K,  the  spot  is  located  using  a 
thresholded  centroid  algorithm 


II 

M 

/ 

j,iBK 

Pk 

= 

/ 

j,ieK 

^Ji>h 

•  The  sum  is  calculated  using  only  those  pixels  in 
region,  K,  for  which  the  gain-and-offset  corrected 
pixel  output,  Vjj,  exceeds  the  threshold  value,  Tjj. 
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Wavefront  Reconstruction 


•  The  subaperture  slopes  (phase  gradients),  s,  are  related  to  the 
wavefront  phase  points,  (i,  by  the  so-called  gradient  matrix,  r 

s  = 

•  An  estimate  of  the  wavefront  phase  can  therefor  be 
“reconstructed”  from  a  measurement  of  the  slopes 


•  H  is  the  wavefront  reconstructor 

•  r  (and  therefor  H)  are  be  determined  a  priori  from  a  knowledge 
of  the  subaperture/actuator  geometry,  which  is  defined  by  the 
user  during  initialization 

•  It  is  possible  to  reconstruct  Zernike  coefficients  directly  from  the 
subaperture  slopes,  without  first  reconstructing  phase.  This 
saves  time. 
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Slope  Discrepancy 


•  Geometry  lets  us  calculate  the  set  of  slopes 
consistent  with  the  calculated  phases: 

A  A 

5  =  ro 

•  The  slope  discrepancy,  \s-sf,  is  related  to  the  amount 
of  noise  in  the  original  slope  measurements. 

•  The  true  noise  is  greater  than  the  slope  discrepancy: 

-  There  are  more  slope  measurements  than  phase  points,  and 

-  The  reconstructor  manifests  its  own  noise  gain  in  converting 
slope  measurements  to  phase. 

•  After  correcting  for  these  two  effects,  we  get  a  good 
estimate  of  the  wavefront  phase  reconstruction  error. 
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Matched  Filter  Subaperture  Spot  Finding 


•  We  start  with  a  pre-defined  hypothetical  function, 

G(x,y),  which  looks  much  like  the  subaperture  spot 

•  We  minimize  the  function 

Fnn,  =  J 

jJ^K 

-  we  start  with  a  good  guess  for  the  spot  location  (perhaps 

from  the  centroid):  y„),  and  evaluate  F^.„. 

-  we  choose  four  values  for  y„.)\  left,  right,  above,  and 

below  (x^.,  y„.),  and  evaluate  for  each 

-  we  perform  parabolic  interpolation  in  each  dimension  to  find 
the  true  minimum  of  to  sub-pixel  accuracy. 
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Matched  Filter  Versus  Centroid 


•  The  advantage  of  the  centroid  is  that  it  is  fast. 

•  The  advantage  of  the  matched  filter  is  that  it  is 
optimal 

-  it  emphasizes  the  bright  pixels,  which  have  the  greatest 
signal-to-noise  ratio, 

-  so  it  is  as  immune  to  noise  as  is  theoretically  possible, 

-  especially  if  time  permits,  wherein  continued  iterations  can 
find  the  true  minimum  of  F„,n . 

•  Depending  on  the  signal  levels  and  accuracy 
requirements  of  the  system,  and  the  need  for  speed, 
one  of  these  algorithms  is  better  suited  than  the 
other. 

•  We  have  provided  a  choice 
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Features  of  the  Wavefront  Sensor 


•  The  intensity  and  phase  information  can  be  used  to  estimate 
quantities  such  as: 

-  rms  phase  error 

-  peak  to  valley  phase  error 

-  WFS  accuracy  (using  slope  discrepancy) 

-  beam  quality 

-  Zernike  polynomials 

-  diffraction  point  spread  functions 

-  MTF 

-  X  and  Y  overall  tilt 


•  Wavefront  sensor  results  can  be  presented  as; 

-  contour  plots 

-  isometric  plots 

-  vector  plots  (showing  gradients)  and 

-  simulated  interferometric  fringes  (as  with  an  interferometer) 
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Wavefront  Sensor  Hardware 


•  Either  Sony  XC-77  RS-170  camera  (~$1K) 

-  CCD  focal  plane  with  640x480  13  pm  square  (effective)  pixels 

-  6x8  mm  array  size  -  limits  measurable  beam  size  to  6  mm 

•  or  Pulnix  TM-1001  digital  camera  (~7K) 

-  CCD  focal  plane  with  1024x1024  9  pm  square  pixels 

-  9x9  mm  array  size  -  permits  beam  size  to  9  mm 

•  Lenslet  array  (~$1  K) 

-  F/25  elements 

-  0.5  mm  or  0.2  mm  diameter  subapertures 

•  MuTech  frame  grabber  (~$2K)  +  (~$1  K  for  digital  input) 

•  Pentium-class  computer  (~$2k) 
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Sensor  Noise  Sources 


•  The  main  sources  of  noise: 

-  signal  shot  noise  (when  signals  are  low) 

-  dark  current 

-  camera  read-out  noise 

*  preamp  &  digitizer 

*  pixel  jitter  (caused  by  inaccuracies  in  horizontal  line 
synchronization) 

-  fixed  pattern  noise  (when  not  compensated) 

-  camera  pixel  structure  noise  (variations  in  pixel  spacing) 

-  quantization  noise  (8  bits) 

•  In  the  case  of  strong  signals,  pixel  structure  noise 
dominates  when  the  spots  are  not  near  the  calibrated 
spot  locations 
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Sample  Wavefront  Measurements 


•  Wavefront  measurements  were  made  with  12 
subapertures  across  a  6  mm  diameter  circle 

•  The  lenslet  array  was  15.2  mm  from  the  CCD 

•  At  1.064  pm  wavelength,  the  relationship  between 
spot  displacement  and  subaperture  tilt  was  2.48 
pixels/wave 

•  The  spot  size  (full-width  at  half-maximum)  was  about 
5  pixels 

•  The  nominal  separation  between  spots  (and  therefor 
subapertures)  was  38.46  pixels. 

•  This  allows  the  measurement  of  up  to  15  waves  of 
phase  difference  across  each  subaperture 


Laser  Cavity  Wavefront 


•  Thermal  loading  of  high  energy  rod  causes  aberration 

•  Significant  trefoil  (  cos  30)  is  present 
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Adjusting  for  Collimation 


•  As  part  of  the  laser  system  alignment,  various  stages  of  the 
optics  were  brought  into  collimation 

•  Shown  above  is  a  series  of  reconstructed  wavefronts,  during 
which  a  lens  is  adjusted  to  bring  the  beam  into  best  collimation 

•  A  full  wave  of  focus  aberration  was  removed  (lower  right) 
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Aligning  Multiple  Amplifier  Rings 


•  The  wavefront  sensor  is  very  sensitive  to 
overall  beam  tilt 

•  When  tilts  where  not  removed  from  the 
wavefronts  during  reconstruction,  the  results 
were  used  to  detect  alignment  errors  between 
the  different  amplifier  rings 

•  Each  amplifier  ring  was  brought  into 
alignment  with  the  master  oscillator 
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Conclusions 


•  The  wavefront  sensor  is  able  to  provide  intensity  and  phase 
information  on  a  laser  system  in  quasi-real  time 

•  Significant  characteristics,  such  as  Zernike  modes,  Seidel 
coefficients,  WFS  accuracy,  and  wavefront  statistics,  can  be 
displayed 

•  The  rapid  update  of  optical  measurements  makes  alignment 
and  adjustment  quite  easy 

•  A  choice  of  spotfinding  and  reconstruction  algorithms  lets  the 
user  tradeoff  accuracy,  dynamic  range,  and  speed 

•  Many  other  options  (such  as  tilt  and/or  focus  removal)  aid  in 
wavefront  analysis 

•  This  wavefront  sensor  proved  to  be  an  excellent  tool  in  the 
development  and  understanding  of  laser  systems 
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COMPUTER  CONTROLLED 
ACTIVE  MIRRORS  AND  LENS  FOR 
INDUSTRIAL  APPLICATIONS 


Bimorph  adaptive  mirror 

1.  Kokorowsky  S.A.  “Analysis  of  adaptive  optical 
elements  made  from  piezoelectric  bimorph”,  JOSA,  69, 
1979 

2.  Lipson  S.G.,  Steinhaus  E.  “Bimorph  piezoelectric 
flexible  mirror”,  JOSA,  69,  1979 


Main  Requirements  to  adaptive  mirrors: 

-  deformable  range  of  mirror  surface 

-  small  number  of  controlling  channels 

-  simplicity  of  production  and  application 

-  low  cast  of  mirrors 

-  high  temperature  stability 


Glass  bimorph  mirrors 
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reflective 
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defocus  control 
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piezodisk 

controlling 

electrodes 


Main  features : 

aperture  -40nini 
surface  deformation 
number  of  actuators 
mirror  quality 
frequency  bandwidth 
reflecting  coefficient 
reflective/  protective  coatings 
control  voltage 
hysteresis 


-7p 

-8,13,17 
-0.1  (PV) 

-2  kHz 

>98%  at  A,=0.63 
A1  /  AI2O3 
±300  V 
<10% 


Water  cooled  bimorph  mirrors 

Reflective  plates 


Main  features:  aperture  - 100  mm 

deformation  of  surface  -  6  p 
number  of  actuators  - 1, 9, 18 
mirror  quality  -  0.3|a. 
frequency  bandwidth  -  >3  kHz 
material  -  copper,  molibdenum 
reflecting  coefficient  -  >99%  at  A,=10.6  jj, 
reflective/protective  coatings  -  copper/ AI2O3 
control  voltage  -  300V 


BIMORPH  ADAPTIVE  MIRROR 


cover 


mirror  frame 

reflective  plate 

piezoceramic 

tilt  control  electrodes 

connector 

glass  tube 
wires 

back  plate 

tilt  control  electrodes  holder 


General  design  of  the  adaptive  mirrors  with  tilt  and  piston  control 


Schemes  of  control  electrodes  of  the  adaptive  mirrors  for: 

a)  defocusing  control; 

b)  astigmatism  and  defocusing  control; 

c)  second-order  aberrations  control; 

d)  tilt  and  piston  control 
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Main  features  of  adaptive  mirrors 


Aperture  -  40  mm 
Stroke  (peak-to-valley)  -  4  mkm 
Generation/compensaton  of  aberrations: 
piston  - 10  mkm 
tilt  -  0.8  mrad 
defocus  -  4  mkm 
astigmatism  -  6  mkm 


coma  -  0.8  mkm 
spherical  aberration-  0.75 
Number  of  controlling  electrodes 
Histeresis  - 
Drive  voltage  - 
Capasitance  - 
Full  amplitude  bandwidth  - 
Surface  figure  - 
Residual  aberration  - 
Thermostability  of  surface  shape 
Material- 

Reflecting  coating  - 


mkm 

-  1,  9,14,18 
<10% 

±300  V 
50  nF 
IkHz 

XI4  (X<=0.633mkm) 
XllO  (X=0.633mkm) 
0.04  mkmPC 
glass  LK5,  Si 
Cr/AySi02 


Reflection  coefficient  -  >  90%  in  the  wavelength 

range  0.4  -j-  0.7  mkm 


mirror 

substrate 

material 

a ,  X 

wm 

m, 

^m/l«C 

mirror 

thickness, 

mm 

mirror 

aperture, 

mm 

quartz 

0.55 

-0.21 

2.5 

40 

Si 

2.54 

-0.035 

2.6 

40 

piezoceramic 

«  3 

0 

0.35 

40 

3.5 

0.04 

3.0 

40 

Temperature  dependences  of  defociis  changing  for  different  kinds 
of  mirror  substrate  materials  (m  -  defocus  changing  value  per  1«C) 


Surface  profile  of  influence  function  of  13  electrodes  bimorph  mirror. 
+150v  and  -150V  voltage  was  applied  to  periferical  electrode  1. 
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Type  of 
mirror 

Residual 

aberrations, 

pm 

P-Vof 

mirror 

surface, 

pm 

Sensitivity 
of  electrode 
1, 

pm/V 

Sensitivity 

of 

electrodes 

2-9, 

pm/V 

Sensitivity 

of 

electrodes 

10-13, 

pm/V 

Defocus 

0.074 

0.27 

11.2*10-3 

- 

- 

DefocusH- 

astigmatism 

0.063 

0.38 

13.86*10-3 

4.41*10-3 

- 

2nd  order 
aberrations 

0.049 

0.24 

11.5*10-3 

2.5*10-3 

1.0*10-3 

Residual  aberrations  and  sensitivity  of  the  adaptive  mirrors 


tilt  control  value,  mrad 


b 

a 


Control 
voltage,  V 


Experimental  curves  showing  tilt  c(y*.li  oI  ranges  by  using  of 
a  -  one  electrode 
b  -  two  electrodes 
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CONTROL  SYSTEM 


AMPLIFIER 

BLOCK 


Configuration  of  RTV-300  system 

Channels  20 

Max.  frequency  in  all  channels  500  Hz 

Output  voltage  ±300  V 

Max.  output  current  of  a  single  channel  80  mA 
Max.  capacitive  load  per  channel  0.2  |lF 

Velocity  of  vpltage  build-up  0.2  V/1  jls  at  0,2  jlF  load 

.  >  ■  0.7  V/l|is  at  0,02  |j,F  load 

Nonuniformity  of  amplification 
factory  in  channels  5% 

Sequential  channel  RS-485,  9600  bit/s 

Cohlsumed  power  200  W 

-  ■'*. 
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SOFTWARE 


Software  includes  the  special  subprogram-functions: 

Zero  -  to  set  0  volts  at  all  the  electrodes. 

Read  -  to  read  the  given  set  of  voltage  values  at  the  elecrodes 

from  the  file. 

Write  -  to  record  the  current  voltage  values  at  the  electrodes 
into  the  tile. 

Gain  -  to  enter  a  factor  of  amplification  for  the  set  of  voltages. 
Rotation  -  to  enter  a  rotation  angle  of  voltage  arrangements. 

Real  time  -  to  go  to  a  submenu  of  real  time  operations  . 


lens  frame 

piezoelectric  plate 
glass  plate 
glass  lens 

immersion  liquid 
back  plate 
front  plate 


General  design  of  controlled  lens  and  scheme 
of  control  electrodes 


Residual  aberrations 

|Lim 

0.12 

Sensitivity  of  electrode  1 

fim/V 

0.9*103 

Defocus 

The  same  electrical  voltage  is  applied  to  all 
the  electrodes  simultaneously. 

jurnA^ 

1.97*103 

IB8IESM9 

0.6 

Astigmatism 

The  positive  electrical  voltage  is  applied  to 
electrodes  1, 5  and  the  negative  one  to 
electrodes  3, 7  simultaneously. 

limlV 

3.35*103 

|Lim  at  300V 

1.0 

Characteristics  of  the  controlled  lens 


Applications  of  Shack-Hartmann  wavefront  sensor 

Wenhan  Jiang,  Zeping  Yang,  Hao  Xian,  Lingtao  Jiang,  Xuejun  Rao,  Bin  Xu 

Institute  of  Optics  and  Electronics,  Chinese  Academy  of  Sciences 
P.  O.  Box  350  ShuangliUf  Chengdu  610209,  China 
Fax: +86-28-5180070,  e-mail:  whjiang@cdisn2.cdisn.cd.sc.cn 

Abstract 

Shack-Hartmann  wavefront  sensor  can  measure  both  intensity  distribution  and  phase  gradients  of  optical  filed  with  high 
frame  rate.  It  is  not  only  widely  used  for  wavefront  sensing  in  adaptive  optics  system  but  also  for  inspection  of  optical 
components  and  systems,  diagnosis  of  laser  beam  quality  and  atmospheric  turbulence  measurement.  We  have  developed 
several  S-H  sensors.  The  results  of  measurements  by  using  these  sensors  in  different  areas  are  reported  in  this  paper. 

1.  Introduction 

Shack-Hartmann  (S-H)  wavefront  sensors  are  widely  used  in  adaptive  optical  systems.’  By  using  a  proper  data  acquisition 
and  processing  equipment,  the  S-H  sensor  can  measure  spatial  and  temporal  characteristics  of  intensity  and  wavefront  phase 
profile  of  optical  beam  simultaneously,  providing  detailed  dynamic  information  of  laser  beams  and  optical  systems.  From  these 
data  the  detailed  near-field  and  far-field  specifications  can  be  obtained.  It  can  be  calibrated  by  a  reference  beam  in  advance.  In 
field  measurement,  it  is  not  necessary  to  have  a  on-line  reference  and  is  not  as  sensitive  to  the  environment  as  interferometer.  It 
is  a  powerful  tool  for  laser  beam  diagnosis  and  optical  quality  control.  We  have  developed  several  types  of  S-H  sensors  for 
different  applications,  such  as  inspection  of  optical  components  and  systems,  diagnosis  of  laser  beam  quality,  and  atmospheric 
turbulence  measurement.  In  this  paper  some  results  of  these  applications  are  presented. 

2.  Construction  of  S-H  sensors 

The  S-H  sensor  for  laser  beam  diagnosis  and  atmospheric  turbulence  measurement  is  shown  in  Fig.  la.  Lenses  1  and  2  are 
the  aperture  matching  optics  for  converting  the  beam  aperture  to  the  aperture  of  lenslet  array  3.  Lenslet  array  divides  the 
aperture  into  subapertures  and  forms  a  spot  array  at  the  focus  plane.  Lens  4  matches  the  size  of  spot  array  to  the  detecting  area 
of  CCD  5.  By  A/D  converter  6  the  output  of  CCD  is  acquired  by  the  computer  7.  The  whole  system  is  calibrated  by  a  collimator 
as  a  reference  and  the  centroids  of  reference  spot  array  are  stored  in  the  computer. 

The  S-H  sensor  for  inspection  of  optical  components  and  systems  is  shown  in  Fig.  lb.  A  laser  10  is  used  as  the  light 
source  and  is  expanded  and  collimated  by  lenses  9  and  1 .  A  flat  mirror  is  used  as  the  reference  during  the  calibration  procedure. 
When  the  S-H  sensor  is  used  for  inspecting  spherical  elements,  the  beam  is  focused  by  a  lens  and  calibrated  by  a  spherical 
mirror. 


Figure  1.  Structure  of  Shack-Hartmaim  sensor 


3.  Laser  beam  diagnosis 


We  have  developed  S-H  sensors  for  UV,  visible  and  IR  wavelengths  and  used  for  laser  beam  diagnosis  and  optical  system 
inspection.  The  software  for  wavefront  reconstruction  firom  the  measured  wavefront  slopes  in  subapertures  using  different 
algorithms  has  been  developed.^  From  the  reconstructed  wavefront  and  intensity  distribution,  various  parameters  such  as  p-v 
and  rms  phase  error,  Zernike  coefficients,  PSF  and  MTF,  encircled  energy  and  can  be  calculated. 

3.1.  Beam  quality  of  a  non-stable  cavity  laseH 

In  order  to  obtain  the  dynamic  change  of  beam  quality  of  a  non-stable  cavity  laser,  wavefront  profiles  of  the  laser  beam 
were  measured  by  a  S-H  sensor.  The  sampling  rate  is  about  12  frames/sec.  One  sample  (the  5th  frame)  of  the  wavefront  phase 
error  is  shown  in  Figure  2.  The  coefficients  of  Zernike  polynomials  of  the  first  1 6  frames  are  shown  in  Figure  3  as  a  function  of 
time,  among  which  the  astigmatism  in  X  and  Y  ( terms  4  and  5  )  directions  are  the  largest  two  aberrations.  Figure  4  shows  the 
rms  of  the  wavefront  phase  of  each  samples,  the  temporal  changes  of  the  wavefront  phase  error  are  obvious.  The  upper  curve  is 

the  total  wavefront  error  Otot  and  the  lower  curve  is  the  dynamic  changes  of  wavefront  ,  oo  is  the  static 

error. 


Figure  5  is  the  PSF  pattern  for  far  field  of  16  frames  computed  from  measured  wavefront  phase  errors  and  intensity 
distribution  data  at  near  field.  The  change  of  R/Ro,  the  ratio  of  beam  radius  to  the  diffraction  limit,  is  shown  in  Figure  6. 
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Figure  3.  Temporal  change  of  Zernike  coefficients. 


Figure  4.  Temporal  change  of  rms  wavefront  errors. 
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Figure  6.  Change  of  the  ratio  of  beam  radius  to  the  diffraction  limit. 


3.2.  Beam  diagnosis  of  a  high  energy  laser  for  inertial  confinement  fusion  (ICF) 

In  the  laser  driver  for  ICF,  high  intensity  concentration  at  the  focus  and  high  effrciency  for  harmonic  transform  are 
desired.  High  power  laser  beam  with  minimum  wavefront  error  is  highly  demanded.  To  improve  the  beam  quality  and  analyze 
its  influence  factors,  a  careful  investigation  of  the  wavefront  aberration  is  necessary.  The  main  sources  of  wavefront  aberration 
in  ICF  driver  are:  (1)  static  error  of  the  optical  system,  including  figure  errors  of  optical  components,  inhomogeneity  of  optical 
materials  and  alignment  errors,  (2)  thermal  effect  of  pump  light  on  amplifiers  and  its  recovering  process,  (3)  thermal 
turbulence  in  the  optical  path  and  (4)  non-linear  effect  of  the  high  energy  laser. 

Figure  7  shows  the  optical  path  of  an  ICF  facility  in  the  Shanghai  Institute  of  Optics  and  Fine  Mechanics,  Chinese 
Academy  of  Sciences.  The  pulse  Nb  glass  laser  is  amplified  by  a  series  of  amplifiers.  In  order  to  investigate  the  wavefront 
errors  of  this  facility,  a  H-S  sensor  has  been  developed  (Fig.  la).  The  main  specifications  are: 

Size  of  lenslet  array:  26x26, 

Dynamic  range:  1.5X/subaperture  (A,=1.06nm), 

Uncertainty  of  measurement:  X/20  rms. 

Pixel  number  of  the  CCD:  5 12x512. 

Apertures  of  matching  lenses:  O50mm. 

The  wavefront  errors  measured  at  the  point  1  are  shown  in  Figure  8.  Figure  8a  is  the  static  wavefront  error  of  the 
amplifiers  al,  a2,  a3,  and  a4  measured  by  passing  a  continuous  laser  through  the  amplifiers.  The  main  components  are  defocus, 
astigmatism  and  coma.  Figure  8b,  8c,  and  8d  are  the  errors  of  a  pulse  amplified  by  al,  al  and  bl,  and  a3  respectively.  The 
static  error  is  the  main  source  in  these  stages,  the  thermal  effect  is  also  detected  but  not  significant.  In  these  stages  the  power  of 
laser  pulse  is  not  as  high  as  in  the  end  part  of  the  amplifier  chain. 

We  measured  the  recovering  process  of  the  thermal  effect  at  point  2.  The  aperture  of  the  amplifier  is  100mm.  After  a 
liigli  power  pulse  passing  through  the  system,  the  aberrations  of  in  <I>50mm  area  at  the  edge  of  the  amplifier  aperture  are 
measured  continuously.  Figure  9  shows  the  temporal  change  of  defocus  in  16  minutes  after  the  pulse.  An  attenuate  periodic 
change  of  defocus  is  obviously.  It  is  due  to  the  thermal  conduction  and  balance  process  from  the  edge  to  the  center. 


1  2  3 

Figure?.  Optical  path  of  the  ICF  facility  in  the  Shanghai  Institute  of  Optics  and  Fine  Mechanics. 


Focus 


Figure  9.  Change  of  defocus  after  laser  pulse. 

4.  Inspection  of  optical  elements  and  systems 

The  traditional  methods  for  inspection  of  large  optical  elements  and  systems  are  interferometer  and  Focaut  knife  edge. 
The  influence  of  environment  to  the  interferometric  inspection  for  large  elements  is  serious  and  the  Focaut  test  is  difficult  to 
give  quantitative  results.  The  S-H  sensor  is  less  sensitive  to  the  environment  as  the  interferometer  because  the  reference  is  off¬ 
line  and  the  dynamic  effects  can  be  eliminate  by  averaging  of  multi-firames.  A  H-S  wavefront  sensor  for  optical  inspection  is 
developed  (Figure  lb).  In  this  sensor  a  He-Ne  laser  is  projected  as  the  light  source. 

Figure  10  is  two  layouts  for  inspecting  large  spherical  and  flat  mirrors.  Figure  1 1  is  a  measurement  of  a  spherical  mirror 
of  diameter  <I)600mm.  Figure  12  shows  an  inspection  for  an  optical  system  consisting  of  two  aspherical  mirrors  (OdSOmm  and 
120mm). 


Figure  10.  Two  layouts  for  inspecting  large  spherical  and  flat  mirrors. 
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Figure  1 1.  A  measurement  of  a  spherical  mirror  of  dianieter  <t)600mm. 
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Figure  12.  A  measurement  of  an  optical  system. 

5.  EQgh  resolution  measurement  of  atmospheric  turbulence'* 

Atmospheric  turbulence  results  in  fluctuations  of  arrival  angle.  By  a  S-H  sensor  with  high  frame  rate  CCD,  the 
atmospheric  turbulence  can  be  sampled  with  high  spatial  and  temporal  resolutions.  A  large  volume  of  data  can  be  measured 
and  the  detailed  parameters  of  turbulence  can  be  obtained.  The  key  is  to  compress  the  data  in  real-time  for  storage.  We 
developed  a  S-H  sensor  with  higli  frame  rate  CCD  and  a  multi-DSP  processor  and  successfully  used  this  sensor  for  atmospheric 
turbulence  measurement. 

In  this  system,  a  37  subaperture  lenslet  array  is  arranged  in  triangle  pattern,  forming  focus  spot  in  7  rows  (Figure  13). 
Each  row  has  4  to  7  spots.  A  CCD  camera  with  830  Hz  frame  rate  and  128x128  pixels  are  used  to  detect  the  Hartmann  spots.  If 
the  spot  images  were  stored,  the  data  rate  would  be  around  14  Mbytes/sec.  The  interested  data  are  the  centroid  positions  of 
spots  which  are  the  measures  of  wavefront  slopes.  If  we  only  store  the  slope  data  of  subapertures,  the  data  rate  is  only  0.061 
Mbytes/sec  and  the  ratio  of  data  compression  is  22 1 ;  1 . 

A  multi-DSP  processor  developed  by  us  is  used  to  computing  the  wavefront  slopes  in  subapertures  in  real-time  and  the 
slope  data  are  stored  in  the  extension  memory  of  a  personal  computer.  Figure  14  is  the  diagram  of  the  system.  A  personal 
computer  is  used  as  the  host  computer.  There  are  4  processing  modules  in  the  system.  In  each  module  a  TMS  320C50  DSP  is 
used  as  processor  and  a  dual-port  memory  (2Kxl6bits)  as  the  communication  buffer  between  host  computer  and  the  DSP.  Data 
acquisition  circuit  receives  image  data  from  the  CCD  camera,  and  load  them  to  the  image  buffers  in  processing  modules 
according  the  address  from  the  Lock-up-table  (LUT).  LUT  is  a  memory  whose  input  address  is  the  coordinates  of  the  pixel  and 
output  data  is  the  code  indicating  to  which  buffer  the  image  data  should  be  send  to. 


Four  processing  modules  work  in  parallel.  Each  module  processes  1  to  2  spots  of  each  row.  In  Figure  13  the  dispatching 
scheme  of  the  processing  task  is  given.  The  first  code  in  each  square  indicating  the  active  filed  of  subaperture  is  the  number  of 
DSP  module  for  processing  this  aperture  and  the  second  code  is  the  number  of  subaperture  which  is  processed  by  this  DSP.  As 
much  as  9  to  10  subapertures  are  processed  by  each  module  in  one  frame  period.  By  using  a  16  Mbytes  expanded  RAM,  as 
much  as  200,000  frames  can  be  stored  in  one  measurement  period.  The  atmospheric  turbulence  can  be  sampled  by  37 
subapertures  with  frame  rate  830  Hz  for  about  260  seconds  continuously. 


From  these  data  the  detailed  specifications  of  optical  turbulence  can  be  deduced.  From  each  frame  the  wavefront  phase 
profile  can  be  reconstructed,  the  structure  function  D  and  coherence  length  ro  can  be  obtained.  By  large  number  of  frames  the 
statistic  parameters,  such  as  power  spectrum  of  arrival  angle  fluctuation  both  in  subapertures  and  in  the  whole  aperture  and  the 
probability  density  function  of  ro,  can  be  obtained.  Figure  15  is  the  power  spectrum  of  arrival  angle  fluctuation  measured  in  a 
500m  horizontal  path.  Figure  16  shows  the  probability  density  function  of  ro.  The  circles  are  experimental  results  and  curve  is 
fitted  by  a  theoretical  model. 


6.  Conclusion 


S-H  sensors  are  widely  used  in  Adaptive  Optical  systems  as  wavefront  sensors.  The  capability  of  measuring  both  intensity 
distribution  and  phase  gradients  of  optical  field  with  high  frame  rate  makes  it  a  powerful  tool  in  many  other  applications.  In 
this  paper,  only  a  part  of  applications  is  reported.  We  believe  that  S-H  wavefront  sensor  will  have  more  application  fields  in  the 
future. 
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This  signal  is  sampled  in  time  at  each  pixel  -  the  signal  is  examined  at  the  beat 
frequency  Aco,  to  extract  the  amplitude  and  phase  components. 

A  complex  exposure  is  formed  in  the  computer 


DOUBLE-EXPOSURE  HETERODYNE  IMAGING: 
LINE-OF-SIGHT  DEPTH  CHANGE 
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HETERODYNE  IMAGING  LABORATORY 

SCHEMATIC 


CAMERA  RESOLUTION:  640  x  480  pixels 

CAMERA  FRAME  RATE:  30  Hz 

BEAT  FREQUENCY:  8  Hz 

FRAMES  DEMODULATED:  16 


OPTICAL  LAYOUT 
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DOUBLE  PULSE  HETERODYNE 
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MEM  M IRROR  MEASUREMENT 
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/Applications  of  S-H  wavefront  sensing 

•  white  light  illumination 

-  MEMS  mirrors 

-  rear  view  mirrors 

-  (atmospheric  turbulence) 

•  laser  illumination 

-  aspheric  optics 

-  laser  gain  media 

-  atmospheric  turbulence  (laser  comm) 
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Laser  comm  field  trials 


V. 


ACfaA 


Adaptlvabptict  Auodattt 


To  determine  requirements  for  an  adaptive  optics 
system  to  be  used  in  laser  communication  over  a 
horizontal  path. 

-  determine  requirements  through  field 
measurements 

Demonstrate  laser  beam  compensation  in  a 
communications  experiment 

-  tip/tilt 

-  tip/tilt  +  higher  order 
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Laser  test  range 


^  Test  data-  Hartmann  spot  images 


Data  rates: 

-  IKHz  frame  rate 

-  1024/2048  contiguous 
images 

Intensity  statistics 

Wavefront 

reconstructions 

-  Adaptive  optics 
requirements 
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Data  reduction 


Intensity  fading  data  reduction 


Spot  image  with 
boxes 

-  Sum  of  intensities 
within  subaperture 

-  (WAVELAB  kernel 
extension) 

Average  of  spot 
images  used  as 
reference  for 
reconstruction 
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Strehl  ratio  comparisons 
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AdaptivtCptles  AuodatM 


VI 

0&45-1997 


X  =  1.3  jim 
Range  =  1  Km 

no  AO  needed  (sOl),  full  AO  needed  (s02) 


06-05-1997 


Page  11 


Shack-Hartmann  wavefront  sensing  for 

optical  metrology 


^  AO  system  recommendations 


Adaptive  optics  will  enhance  transmission  in 
85%  of  cases  at  1  Km  range. 

AO  will  enable  transmission  at  1.8  and 
2.25  Km  ranges. 


Real  time  compensation  schematic 
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Hardware  benchmark  speed  I 

Same  software  but  with  interrupt  lines 
monitored 

200  MHz  Pentium  Pro/DOS/32bit  C 


Total  period 

Frame  grab/transfer 
Computation 
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Summary 


Shack-Hartmann  wavefront  sensing  can  be 
automated  for  metrology 

Partial  atmospheric  corrections  can  enhance 
laser  beam  propagation  over  horizontal  paths 
PC  AO  implementation  feasible 
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ADAPTIVE  WAVEFRONT  SENSOR 
Valeri  A.Tartakovski,  Vladimir  P. Lukin 

Institute  of  Atmospheric  Optics  of  Russian  Academy  of  Sciences  *Tomsk  *  634055  *  Russia 
Fax:  (3822)  259086  ♦  E-mail:  tv@cnod*iao«tomsk*ru 
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The  objective  of  this  paper  is  to  direct  your  attention  on  research  of  wavefront  measurement  for  adaptive 
optics,  taken  under  strong  scintillation  conditions.  Both  the  fundamental  and  technological  problems'  concerning 
this  prionty  area  on  optics,  must  be  investigated.  We  suggest  to  solve  the  phase  measurement  problems  under 
strong  intensity  fluctuations,  in  which  case  there  is  no  continuity  of  the  wavefront,  by  means  of  selective  and 
consecutive  processing  of  the  speckles.  The  expected  practical  result  of  this  solving  is  a  gain  .in  efficiency  in 
optical  systems  for  imaging,  energy  transfer,  atmospheric  communications,  high  power  lasers,  and  astronomy. 

Lomonosov  and  Descartes  had  known  about  the  first  principles  of  the  wavefront  measurements.  Wavefront 
sensing  concerns  testing  any  kind  of  optical  system,  from  telescopes  to  microscopes,  to  industrial  imaging  and 
lasers.  Popular  sensors  are  the  Shack-Hartmaim  and  curvature  sensors.  Other  sensors  employ  lateral  shearing  or 
point  diffraction  interferometry. 

All  these  sensors  may  not  be  complicated  but  they  all  assume  a  priori  the  continuity  of  the  m^ured 
wavefront.  This  may  not  be  true  when  dealing  with  imaging  through  strong  turbulence.  This  is  because  the  best 
place  to  measure  a  wavefront  is  at  a  plane  which  is  Fourier  conjugate  (or  nearly  so)  to  the  focal  plane.  In  this 
measurement  plane,  intensity  is  assumed  to  be  rather  constant.  When  the  turbulence  level  rises,  this  assumption  is 
broken:  the  intensity  changes  (sometimes  at  a  very  fast  rate)  not  only  in  the  focal  domain  but  also  in  its  Fourier 
conjugate,  where  the  wavefront  is  measured.  These  intensity  variations,  also  called  scintillations,  develop  and 
phase  dislocations  or  wavefront  branch  cuts  spring  up. 

Under  these  conditions  the  Hartmann  sensor  fails  because  the  fixed  position,  the  small  number  of  its 
subapertures,  and  the  locations  of  the  foci  are  not  useful  any  more  to  tell  the  wavefront.  Most  of  die  information 
lies  in  the  high  intensity  scintillations  which  might  be  missed  by  the  Hartmann  subapertures.  The  curvature  sensor 
fails  here  as  it  assumes  a  wave  of  constant  intensity,  a  condition  violated  under  scintillations.  The  same  limitation 
applies  to  other  wavefront  sensors.  .  '  : 

The  wavefront  dislocations  disturb  the  processing  of  wavefront  sensors,  since  the  processing  is  based  on  the 
assumed  continuity  of  the  phase  function.  Such  disturbances  can  be  evfcn  due  to  only  one  dislocation.  In  addition, 
wavefront  dislocations  cause  disintegration  of  tlie  wavefront  into  separate  beams  which  might  be  uncorrelated. 
Moreover,  the  amount  of  energy  in  the  neighbourhood  of  dislocations  is  limited.  As  a  result  the  phase  at  these 
points  is  not  well  characterised. 

Hence  measurement  of  high-energy  sections  of  the  wavefront,  and  subsequent  separate  fitting  of  these  sections 
to  the  adaptive  mirror  should  be  the  right  way  for  measuring  and  controlling  the  wavefront  in  such  adaptive 
systems.  -  " . c. 

To  this  end,  it  is  necessary  to  formulate  properly  the  unexplained  concept  of  a  wavefront  under  strong 
scintillation  conditions.  A  proper  mathematical  description,  adequate  to  the  nature  of  this  problem,  has  to  be 
provided.  Having  done  that,  effective  algorithms  have  to  be  developed  for  data  processing  for  this  case. 

*  The  first  stage  will  be  to  find  in  what  way  to  change  the  spatiotemporal  resolution  during  measurement  and 
according  to  the  light  intensity  at  the  input  aperture.  It  should  be  borne  m  mind  that  the  ^eas  with  high  intensity 
must  be  measured  more  precisely  and  processed  with  less  delay,  as  compared  with  the  low  intensity  areas. 

*  The  second  stage  is  to  find  a  way  to  control  the  size  of  the  single  Hartmaim-Shack  subaperture  according  to 
a  level  of  intensity  and  wavefront  fluctuations  so  as  to  change  the  focal  spot  size  and  its  motion  range. 

*  In  the  third  stage,  an  algorithm  will  be  developed  to  restore  the  phase  along  closed  and  unclosed  lines  in  the 
aperture  plane,  going  through  areas  with  high  intensity,  avoiding  areas  with  very  low  intensity,  and  connecting  the 
phase  from  the  separate  areas  to  unique  phase  function. 

4  in  ,the  fourth  stage,  is  to  be  obtained  an  method  of  estimation  of  restored  phase  errors.  .  ... 

*  Finally,  new  forms  of  the  wavefront  regularisation  will  be  developed  in  order  to  realise  common  adaptive 
optics  principles.  For  example,  a  wave  calculated  from  an  er.  xct  phase  but  without  the  true  amplitude,  may  be 


1 


worse,  after  conjugation,  than  some  other  wave  which  phase  is  corrected,  using  the  true  amplitude  information,  so 
that  to  minimise  the  functionals: 

2 

dr^=mm  JJ  <ir=min. 

O  R  O 

The  unique  tracking  measurements  of  phase 

In  the  proposed  wavefront  sensor  the  measurements  of  many  Hartmann's  foci  displacements  are  replaced  by 
tracking  only  one  focal  spot,  obtained  from  the  brightest  part  of  tlie  field.  This  is  a  natural  step  forward  as  well  as 
a  solution  of  many  other  problems  in  which  the  parallel  and  sequential  processing  can  be  or  ^e  exchanged 

Each  high-energy  section  of  the  wavefront  must  be  measured  separately  and  precisely  as  possible  by:means  ,of 
dense  scanning  in  these  areas.  If  this  is  possible  and  there  exists  a  continuous  track  with  sufficiently  intensity 
between  these  areas,  a  constant  phase  shift  will  be  cjUculated  between  them.  If  the  patches  are  not  contiguous  and 
the  beams  are  uncorrelated,  the  shift  of  these  parts  of  wavefront  will  be  determined  by  successive  approximations. 
In  this  case  a  four-dimensional  analytic  signal  allows  to  circumvent  singular  points  for  .connecting 
parts  of  the  wavefront  to  create  a  unique  phase  function.  The  unique  four-dimensional  analytic  signal  is 
defined  as  follows: 

def  ■ 

W(x,y,z,t)  =  U(x,y,z,t)+iV(x,y.z,t),  V(x,y,z,t)  =  iiU(x,y,z,t)  =  HU(x,y,z,t), 

t  Z 

V[p{x,y,z),t]  =  ^U[p{x,y,z),t]  =  liU[p{x,y,z),t],  V[x(T),yix),z(x)]  =  Ht/[(D"'(a^r)] 

X  y  r 

Here  H  is  the  Hilbert  transform  operator;  W(x,  y,  z,  t)  is  the  object  wave  and  the  analytic  signal  (AS);  the 
AS  imaginary  part  is  y,  z,  t),  it  is  Hilbert  transform  of  real  part  lJ{^x,  y,  z,  t)  over  argument  x,  t,  or  z  which 
is  a  propagation  direction;  p{x,  y,  z)  is  the  turned  recording  plane,  x  is  the  parameter  of  the  curvilinear  cross- 
section  in  three-dimensional  space,  on  which  the  wave  phase  is  monotonic,  is  the  spatial  carrier  frequency  on 

this  cross-section.  0(x)  is  a  function  for  the  compression-extension  transformation  of  the  cross-section,  Some 
approximation  of  arg  W{x)  could  be  used  as  this  function,  l^aturally,  the  function  <1)(t)  could  be  known,  if  Vfx) 
had  been  defined  before.  ,  From  fourth  expression  follows  the  existence  of  the  problem  solution  and  the  solution 
itself  can  be  found  by  the  method  of  successive  approximation.  '  ' 

It  gives  a  guarantee  that  the  phase  of  the  four-dimensional  wave  function,  determined  by  the  AS,  introduce  as 
a  function  of  the  line  parameter  x,  on  the  scanning  line,  is  the  corresponding  cross  section,  of  tiie  unique  four¬ 
dimensional  phase  function  determined  by  the  AS  as  a  function  of  the  time  or  tiie  propagatiomdirection. 

The  practical  use  of  these  expressions  is  that  they  enable  one  to  design  the  algorithms  of  phase  measurement  in 
various  spatiotemporal,  one-dimensional  and  multidimensional  cross-sections  and  provide  its  coincidence  with 
the  unique  four-dimensional  phase  in  these  cross-sections.  It  is  possible  now  to  measure  the  two  or  one¬ 
dimensional  sections  of  the  four-dimensional  phase  separately  and  apart  from  each  other. 

This  phase  exists  only  as  a  consequence  of  narrow  bandwidth  of  the  temporal  and  spatial  wave  spectra. 
Determined  by  anal)^c  signal,  the  amplitude  I  W’(x)|  and  the  phase  zxgW{x)  are  invariant  relative  to  the 
replacement  of  the  Hilbert  transform  argument  and  consequently  are  unique.  ,, 

Description  of  possible  concept  of  the  sensor  is  as  follows. 

A  deflector  sliifts  laterally  the  input  beam  into  a  single  Hartmann  subaperture  which  creates  a  focal  spot.  As  in 
any  Hartmann  case,  the  focal  spot  position  (abscissa  and  ordinate)  and  intensity  can  be  measured  by  a  sensitive 
quadrant  detector  or  a  few  pixels  inside  a  sensitive  CCD  camera.  Simultaneously  the  deflector  anigular  position 
serves  as  reference  signals  (abscissa  and  ordinate).  '  '  • 

The  differences  between  the  focal  spot  co-ordinates  and  the  reference  co-ordinates  from  the  deflector  are  two 
spatial  derivatives  of  the  phase  (tip-tilt)  as  a  function  of  time,  as  the  input  beam  is  scanned  over  the  subaperture. 
These  spatial  derivatives  are  multiplied  by  the  temporal  derivatives  of  the  track  which  are  proportional  to  the 
temporal  derivatives  of  the  reference  abscissa  and  ordinate.  .  >  ■  :  . 

The  temporal  derivatives  of  the  track  are  fom.ed  m  two  feedback  loops.  The  loojps  originate  from  the  focal 
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spot  detector  (the  spot  intensity  signal),  then  in  each  loop  there  is  a  temporal  integrator,  connected  with  the 
actuator  of  the  beam  deflector.  These  feedback  loops  maximise  some  objective  (cost)  function  and  realise  the 
tracking  law.  The  products  of  two  derivatives  are  added  and  integrated  over  time.  The  result  of  this  integrator  is 
the  phase  as  a  function  of  time.  The  equation 


5cp  ^  5(p  dy^ 
dx  dt  dy  dt) 


is  used  to  integrate  the  phase  on  the  trajectory  line,  which  is  defined  by  the 


parametric  equations  x  =  x(r),  y  =  >’(?). 


The  technical  feasibility  of  the  sensor  is  not  limited  in  principle.  There  are  some  limitations  in  the  tracking 
speed,  due  to  the  deflector  actuators.  However,  these  limitations  are  less  as  time  goes  by.  Deformable  mirrors  is 
now  available  off  the  shelf  and  their  price  keeps  dropping.  They  can  reach  a  few  kHz  over  accepteble  areas. 
Thus,  applicability  of  the  sensor  is  limited  by  the  level  of  fluctuations  in  the  medium,  and  by  the  strengtli  of  the 
reference  signal. 


Fig.  1.  Scan-Hartmann  wavefront  sensor 


Object  beam  1.  Reference  beam  2.  Optical  summator  3.  Deflector  with  actuators  4.  Scale  optical  system  5. 
Only  subaperture  Hartmann  screen  6.  ,  Lens-Fourier  transformer  7.  Coordinator  to  track  total  energy  and 
displacement  of  two  focal  spots  8.  Summator  to  calculate  the  object  spot  relative  displacement  9,  10. 
Multiplicator  to  form  the  phase  derivative  by  tracking  parameter  11,  12.  Integrator  to  form  the  phase  as  a 
function  of  tracking  parameter  13.  Generator  of  derivative  of  trajectory  by  tracking  parameter  14,  15.  Integrator 
of  trajectory  to  form  the  first  and  second  actuator  control  16,  17.  Summator  for  control  forming  of  the  third 
deflector  actuator  18. 


A  posteriori  tstim&iion  of  the  restored  phase  error 

The  most  interesting  property  of  the  Scan-Hartmann  wavefront  sensor  operation  is  a'possibility  to  estimate  the 
restored  phase  error  a  posteriori  using  differing  and  orthogonal  scanning  in  the  same  areas.  •  It  is 'possible  to  use 
for  this  goal  the  double  orthogonal  evolvent  Fig.3.  It  is  possible  to  estimate  the  restored  phase  error  without 
using  the  statistical  samples,  when  only  one  two-dimensional  sample  is  available.  It  is  reached  using  some 
orthogonal  scanning  directions. 

Let  us  consider  the  following  phase  functions;  an  original  phase  %{i,  j),  phase  restored  firom  x-cross-section 
of  wave  <^{i,  j),  phase  restored  from  /-cross-section  of  wave  j),  estimate  of  the  original  phase 
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)=[%(j’jy)  + y  )]  /  2,  and  the  following  RMS  values: 

g2  _  iJ^ _  g2  _  ij^ _ I _ 

aSfJfiy)  ' 

iJ^  ij^ 

and  peak-to-valley  values  as  well: 

Poe  =  max|  (PoihJ)  -  <PeO‘.y)|.  Pxy  =  max|(p;^.(/,7)  -  (fy{i,j)\. 

Here  G  is  the  domain  of  the  wave  definition,  iVis  the  number  of  the  readings,  x  and  ynsross-sections  are 
orthogonal  each  another,  for  example,  in  a  simple  case  they  may  be  parallel  to  the  co-ordinate  axes. 

A  functional  dependence  is  discovered  between  the  real  errors  of  restored  phase  z^e » Poe  estirnates  of 

restored  phase  errors  p^y  .  On  Fig.2  is  shown  this  connection  between  the  calculated  and  unknown  in  an 
experiment  real  error.  They  fitted  well  by  the  first  order  regression  lines:  = -0.003(7)  +  1.14(3)e^  , 

0.01(1) +  0.66(6)jr,^^. 

The  rms  deviation  to  last  significant  digit  is  given  in  brackets  here.  The  correlation  coefficient  is  equal  to  0.99 
for  the  normalised  rms  error.  Fig.  2a,  and  0.94  for  the  pv  error,  Fig.26.  The  rms  deviations  from  regression  lines 
are  equal  to  0.01  and  0.04,  respectively. 


Real  errors 


Fig.2.Estimates  of  the  restored  phase  errors  versus  the  real 
errorsofrestoredphase:nonnalisedRMS.error(a)o.Peak- 

to-Valley  error  (6)  □,  parts  of  X.  Regression  lines  — 
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Adaptive  Optics  for  the  Human  Eye 


David  R.  Williams,  Junzhong  Liang,  Donald  T.  Miller  and  Austin  Roorda 
Center  for  Visual  Science,  University  of  Rochester,  Rochester,  NY  14627 


1.  Introduction 

The  first  known  attempt  to  correct  the  optics  of  the  human  eye  occurred  more  than  700  years  ago.  Spectacles  have 
been  used  to  correct  defocus  at  least  as  early  as  the  13th  century  ^>2  and  the  eye's  astigmatism  has  been  corrected  since  the  19th 
century^.  Based  on  subjective  observations  of  a  point  source  of  light,  Helmholtz  argued  that  the  eye  suffered  from  a  host  of 
aberrations  that  are  not  found  in  conventional,  man-made  optical  systems^.  There  have  been  a  number  of  methods  developed 
to  quantify  these  aberrations^"^.  Liang  et  al^  developed  a  technique  based  on  the  Hartmann-Shack  principle  to  measure  the 
eye's  wave  aberration.  This  method  provides  a  rapid,  automated,  and  objective  measure  of  the  wave  aberration  simultaneously 
at  a  large  number  of  sample  points  across  the  eye's  pupil.  By  increasing  the  density  of  samples  taken  of  the  wavefront  slope 
in  the  pupil  and  increasing  the  number  of  Zernike  modes  used  in  reconstruction,  Liang  and  Williams  ^  ®  provided  a  more 
complete  description  of  the  aberrations  of  the  eye.  These  measurements  quantified  the  irregular  aberrations  predicted  by 
Helmholtz's  and  subsequent  investigators^"^.  While  any  of  a  number  of  methods  could  be  used  to  recover  the  wave  aberration 
in  the  eye,  we  focus  here  on  the  Hartmann-Shack  method  due  to  its  recent  success  in  an  adaptive  optics  system  for  the  eye. 

The  effect  of  these  aberrations  on  retinal  image  quality  is  illustrated  in  Fig.  1,  which  shows  the  point  spread  function 
of  the  eye  at  various  pupil  sizes.  These  point  spread  functions  were  calculated  from  wave  aberration  data  obtained  with  a 
Hartmann-Shack  wavefront  sensor^®  for  a  single  subject.  The  effects  of  defocus  and  astigmatism,  which  are  correctable  with 
conventional  spectacles,  have  been  removed.  At  pupil  diameters  of  2-3  mm,  which  correspond  to  the  diameter  of  the  pupil  in 
very  bright  viewing  conditions,  diffraction  is  the  single  largest  source  of  retinal  image  blur.  As  pupil  size  increases, 
aberrations  begin  to  dominate,  reducing  optical  quality.  A  small  pupil  roughly  3  mm  in  diameter  generally  provides  the  eye 
with  the  best  optical  performance  in  the  range  of  spatial  frequencies  (<20-30  cycles/deg)  that  are  important  for  normal  vision. 
If  the  eye's  aberrations  could  be  corrected  across  the  fully-dilated  8  mm  pupil,  resolution  could  be  increased  more  than  2.7-fold 
over  that  obtainable  with  a  3  mm  pupil.  This  increase  in  resolution  allows  single  cells  in  the  retina,  which  have  dimensions 
on  the  order  of  a  few  microns,  to  be  imaged  in  the  living  eye. 

Spectacles  and  contact  lenses  do  not  correct  any  aberrations  except  defocus  and  astigmatism.  There  have  been  various 
attempts  to  correct  these  additional  aberrations  in  the  laboratory.  For  viewing  sinusoidal  gratings,  all  aberrations  of  the  eye 
can  be  avoided  with  interference  fringes  imaged  on  the  retina^  However,  this  technique  is  impractical  for  viewing  stimuli 
other  than  gratings  and  is  of  no  use  for  improving  normal  vision  and  the  quality  of  fundus  images.  Another  approach  is  to  use 
a  contact  lens  to  null  the  refraction  at  the  first  surface  of  the  cornea.  This  approach  has  the  advantages  of  simplicity  and  low 
cost.  However,  contact  must  be  made  with  the  optics  of  the  eye  and  the  effectiveness  of  the  method  depends  on  the  extent  to 
which  the  corneal  surface  dominates  in  the  combined  effects  of  refractive  index  variations  throughout  the  eye's  optics.  A  recent 
attempt  to  use  a  special  contact  lens  to  increase  the  axial  resolution  of  a  confocal  laser  scanning  ophthalmoscope^^  showed 
only  modest  improvement. 

Adaptive  optics  is  well-suited  to  cope  with  the  large  variation  in  the  pattern  of  aberrations  from  eye  to  eye  and  can 
correct  the  wave  aberration  of  the  eye  without  making  contact  with  it.  Dreher  et  al.  first  used  a  deformable  mirror  in 
conjunction  with  the  human  eye^^.  They  succeeded  in  correcting  the  astigmatism  in  one  subject's  eye  based  on  the 
prescription  provided  by  a  conventional  refraction.  The  recent  development  of  a  Hartmann-Shack  wavefront  sensor  for  the  eye^ 
has  improved  the  measurement  of  the  eye's  aberrations^^  Liang,  Williams,  and  Miller^^  have  combined  a  Hartmann-Shack 
wavefront  sensor  with  a  deformable  mirror  to  correct  these  aberrations,  and  have  applied  this  device  to  study  both  the  visual 
performance  and  retinal  images  of  eyes  corrected  to  provide  supernormal  image  quality.  In  normal  eyes  in  which  light  scatter 
is  not  an  important  limitation  of  image  quality,  adaptive  optics  offers  the  possibility  of  obtaining  retinal  image  quality 
limited  only  by  diffraction  at  the  dilated  pupil. 

2.  Measuring  the  Wave  Aberration  of  the  Eye 

2. 1  Hartmann-Shack  Wavefront  Sensor  for  the  Eye 

The  position  of  the  subject's  eye  is  fixed  with  respect  to  the  optical  system  by  having  him  clench  his  teeth  on  a 
dental  impression.  The  optical  axis  of  the  wavefront  sensor  passes  through  a  chosen  reference  axis  in  the  eye,  such  as  the  line 
connecting  the  center  of  the  dilated  pupil  to  the  foveal  center. 

A  collimated,  linearly  polarized  beam  of  a  low  power  laser  incident  on  the  pupil  is' focused  to  a  small  spot  on  the 
retina,  which  the  subject  fixates.  This  small  spot  is  equivalent  to  the  guide  star  or  laser  beacon  in  astronomical  applications 
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of  adaptive  optics.  Its  angular  subtense  is  kept  small  enough  that  its  retinal  image  corresponds  to  the  point  spread  function  of 
the  eye.  The  use  of  a  small  diameter  for  the  beam  at  the  pupil,  approximately  1.5  mm,  ensures  that  the  image  on  the  retina 
is  diffraction-limited  and  relatively  unaffected  by  changes  in  the  refractive  state  of  the  eye.  Liang  and  Williams  used  light  from 
a  5  mW,  632.8nm  He-Ne  laser  and  we  currently  use  750  nm  light  from  a  diode  laser.  Semiconductor  diode  lasers  in  the  near 
infrared  are  particularly  convenient  because  they  are  compact  and  the  point  source  appears  less  bright  for  the  subject.  In  a 
typical  measurement  at  750  nm,  the  retina  is  illuminated  for  a  2  seconds  with  a  total  power  entering  the  eye  of  less  than  15 
micro-watts,  which  is  more  than  an  order  of  magnitude  below  the  ANSI  maximum  permissible  exposure^^. 

Depolarized  light  from  the  retina  is  transmitted  by  a  polarizing  beamsplitter  that  rejects  the  unwanted  reflections 
from  the  relay  lenses  and  the  first  surface  of  the  cornea,  which  retain  their  polarization.  The  relay  lenses  image  the  subject's 
pupil  onto  an  array  of  217  lenslets.  We  represented  the  wave  aberration  with  the  sum  of  65  Zemike  polynomials, 
corresponding  to  aberrations  up  to  and  including  10th  order 

2.2  The  Wave  Aberration  of  the  Normal  Human  Eye 

Despite  the  large  differences  in  the  wave  aberration  between  eyes  of  different  observers,  the  wave  aberrations  of  left 
and  right  eyes  of  the  same  observer  are  often  rather  similar,  when  the  left  eye  imaged  is  flipped  about  a  vertical  axis  for 
comparison  with  the  right  eye.  This  indicates  that  the  aberrations  are  not  random  defects. 

These  aberrations  generally  increase  toward  the  margin  of  the  pupil  and  often  do  not  correspond  to  the  classical 
aberrations  of  man-made  optical  systems.  Fig.  2  shows  the  RMS  wavefront  error  contributed  by  each  Zemike  order  for  a  3.4 
mm  pupil  (triangular  symbols),  a  7.3  mm  pupil  (square  symbols),  and  an  artificial  eye  with  a  6.7  mm  pupil  (round  symbols). 
The  average  RMS  wavefront  error  decreases  monotonically  as  the  Zemike  order  increases  for  both  pupil  sizes  in  the  human 
eye,  though  the  pattern  varies  somewhat  among  individual  observers.  The  RMS  error  for  the  small  pupil  lies  3-4  times  lower 
than  that  for  the  large  pupil  of  real  eyes.  This  illustrates  the  well-known  fact  that  aberrations  grow  with  increasing  pupil  size. 

Fig.  2  also  shows  the  RMS  wavefront  error  for  an  artificial  eye  (round  symbols)  measured  with  the  same  insfrument. 
The  artificial  eye  consisted  of  a  achromatic  doublet  (f =16  mm)  and  a  diffuser  to  mimic  the  retina.  Though  our  measurements 
reveal  some  expected  spherical  aberration  in  the  artificial  eye,  the  mean  RMS  value  of  each  order  averages  about  one  order  of 
magnitude  lower  than  in  human  eyes  for  the  large  7.3  mm  pupil.  Evolution  has  not  yet  created  optics  in  the  human  eye,  at 
least,  that  can  rival  those  fabricated  by  man. 

3.  Adaptive  Compensation 

3.1  Technique 

The  experimental  setup  is  shown  on  figure  3.  A  deformable  mirror  (Xinetics,  Inc.)  compensated  for  the  eye’s  wave 
aberration.  The  mirror  consisted  of  an  aluminized  glass  faceplate  with  37  lead  zirconate-lead  titanate  (PZT)  actuators  mounted 
in  a  square  array  on  the  back  surface.  The  stroke  of  the  mirror  beneath  each  actuator  was  plus  and  minus  2  mm,  allowing  a 
wavefront  shift  of  8  mm  in  the  reflected  beam.  The  mirror  lay  in  a  plane  conjugate  with  both  the  eye's  pupil  plane  and  the 
lenslet  array  of  the  wavefront  sensor.  The  spacing  of  adjacent  actuators  on  the  mirror  was  7  mm.  The  eye's  pupil  plane  was 
magnified  6.25  times  at  the  deformable  mirror,  so  that  the  actuator  spacing  in  the  pupil  plane  was  1.12  mm. 

Compensation  was  achieved  with  closed-loop  feedback  control.  Loops  were  repeated  until  the  rms  wavefront  efror 
could  be  reduced  no  further,  which  usually  required  10-20  loops.  The  improvement  in  the  optical  quality  of  the  eye  provided 
by  the  deformable  mirror  was  evaluated  by  measuring  the  observer's  contrast  sensitivity  for  gratings  viewed  through  the 
deformable  mirror,  and  by  examining  the  quality  of  the  images  of  the  living  retina  taken  through  the  deformable  mirror.  All 
experimental  measurements  used  a  6  mm  artificial  pupil  conjugate  with  the  entrance  pupil  of  the  eye.  The  psychophysical  and 
retinal  imaging  experiments  were  performed  by  sliding  a  mirror  into  the  path  between  the  wavefront  sensor  and  the  deformable 
mirror. 

3.2  Results 

To  determine  which  aberrations  were  corrected  by  the  deformable  mirror,  the  wave  aberration  of  the  eye  was 
decomposed  into  65  Zemike  modes  up  to  10th  order.  A  comparison  of  the  measured  rms  wavefront  error  of  each  Zemike  order 
for  the  eyes  without  adaptive  compensation  with  the  error  after  adaptive  compensation  showed  that  lower  Zemike  orders  up  to 
and  including  4th  order  were  significantly  reduced  while  aberrations  beyond  6th  order  remain  almost  unchanged.  Our  results 
show  that  with  adaptive  compensation  we  not  only  can  correct  the  eye's  defocus  and  astigmatism  but  also  coma,  spherical 
aberration,  and  other  irregular  aberrations  in  the  eye. 

Vision  through  Adaptive  Optics.  Stimuli  such  as  edges  viewed  through  the  compensating  deformable  mirror 
have  a  strikingly  crisp  appearance  consistent  with  the  supernormal  quality  of  the  retinal  image.  Liang,  Williams  and  Miller 
showed  that  observers  with  normal  vision  had  enhanced  contrast  sensitivity  when  looking  through  adaptive  optics.  Contrast 
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sensitivity  measurements  for  two  observers  before  and  after  adaptive  compensation  are  shown  in  Fig.  4.  Filled  symbols  show 
contrast  sensitivity  without  adaptive  compensation,  i.e.  with  the  deformable  mirror  flat.  Unfilled  symbols  show  the  contrast 
sensitivity  when  the  mirror  corrects  the  eye's  aberrations.  At  55  cycle/deg,  neither  observer  could  detect  the  grating  when  the 
deformable  mirror  was  flat  even  at  100%  contrast.  This  was  true  no  matter  what  pupil  size  they  tried  and  despite  attempts  to 
correct  defocus  and  astigmatism  with  trial  lenses.  With  adaptive  compensation,  the  observers  required  about  40%  contrast  on 
average  to  detect  the  grating.  At  another  spatial  frequency,  27.5  cycles/deg,  contrast  sensitivity  was  improved  by  about  a 
factor  of  6  by  adaptive  compensation. 

Imaging  the  Retina.  Estimates  of  the  dimensions  of  microscopic  structures  in  living  human  retinas  have  been 
obtained  using  coherent^O  and  incoherent  light^^  under  conditions  in  which  only  defocus  and  astigmatism  in  the  eye  were 
corrected.  Cone  spacing  of  living  retinas  can  be  estimated  based  on  the  power  spectra  of  retinal  images  in  coherent  light,  but 
laser  speckle  associated  with  coherent  imaging  and  imaging  blur  due  to  the  eye's  aberrations  prevent  the  resolution  of  the 
retinal  structures  in  single  exposures.  Miller  et.  al.  showed,  using  incoherent  light,  that  when  defocus  and  astigmatism  are 
carefully  coirected,  the  cone  mosaic  can  sometimes  be  resolved  in  individual  retinal  images.  Liang,  Williams,  and  Miller^  ^ 
showed  that  adaptive  optics  offers  substantial  improvement  in  the  quality  of  these  images. 

In  their  experiments,  the  observer  fixated  a  crosshair  which  determined  the  retinal  location  to  be  imaged.  When  the 
subject  was  ready,  a  krypton  flash  lamp  delivered  a  4  msec  flash,  illuminating  a  retinal  disc  1  deg  in  diameter.  The  short  4 
msec  exposure  helped  to  prevent  motion  blur  due  to  the  movement  of  retina  during  each  exposure.  The  lamp  output,  which 
was  broad  band  white  light,  was  filtered  with  narrow  band  interferenee  filter  to  the  desired  imaging  wavelength,  "^e  retinal 
irradiance  of  the  flash  was  5.7  mW/mm^  which  is  about  90  times  less  that  the  ANSI  maximum  permissible  exposure  for  this 
retinal  image  size^®.  To  acquire  images  of  the  retina,  a  scientific  grade  CCD  was  positioned  conjugate  with  the  retina,  in  the 
plane  that  previously  contained  the  grating  stimuli  used  in  the  psychophysical  experiments.  Each  CCD  pixel  subtended  8  arc 
sec,  corresponding  to  a  sampling  frequency  at  the  retina  of  450  cycles/deg.  Images  were  acquired  when  the  mirror  was  flat 
(uncompensated  imaging)  and  when  the  mirror  was  warped  to  correct  the  eye's  wave  aberration  (compensated  imaging).  For 
both  the  compensated  and  uncompensated  imaging,  astigmatism  was  corrected  with  trial  lenses  using  the  amount  of 
astigmatism  measured  with  the  wavefront  sensor  when  the  mirror  was  flat.  For  uncompensated  imaging,  the  optimum  focus 
was  determined  by  translating  in  tandem  the  eye  and  the  nearest  lens  to  the  eye  to  achieve  the  highest  contrast  retinal  images. 
For  compensated  imaging,  the  optimum  focus  was  provided  directly  by  adjusting  the  focus  to  compensate  for  the  chromatic 
aberration  between  the  wavefront  sensor  wavelength  and  the  imaging  wavelength. 

Fig.  5  shows  how  image  quality  is  improved  by  using  adaptive  compensation.  The  two  images  are  of  the  same 
retinal  location  (1  degree  eccentricity)  for  subject  JW,  without  and  with  adaptive  compensation,  respectively.  The  image 
without  compensation  was  obtained  with  the  best  focus  and  astigmatism  correction.  After  adaptive  compensation,  the 
photoreceptor  array  is  better  resolved  and  has  a  higher  contrast. 

The  depth  of  field  is  short  compared  to  the  retinal  thickness  and  by  adjusting  the  focus  different  structures  can  be 
observed.  Figure  6  shows  two  images  of  the  same  retinal  location  for  a  single  subject  IP.  The  capillaries  lie  on  top  of  the 
retina  and  appear  as  shadows,  back-illuminated  by  the  out  of  focus  photoreceptor  array.  By  focusing  approximately  100  (im 
deeper  into  the  retina  the  photoreceptor  array  is  clearly  seen  with  the  blood  vessels  appearing  as  dim  shadows  across  the 
image.  In  these  images,  there  is  no  evidence  of  variation  in  image  quality  across  the  image.  This  indicates  that  the 
isoplanatic  patch,  which  limits  the  compensated  field  of  view  to  only  a  few  arcseconds  for  ground-based  telescopes^^,  is  at 
least  as  large  as  the  one  degree  field  of  the  present  instrument. 

Wavefront  sensing  in  the  eye  may  have  important  applications  in  refining  laser  surgical  procedures  such  as  PRK  or 
in  improving  the  design  of  contact  lenses.  Adaptive  optics  allows  microscopic  features  to  be  resolved  in  the  living  human 
retina,  which  may  have  interesting  application  in  basic  and  clinical  research  on  the  eye,  and  ultimately  in  the  diagnosis  and 
treatment  of  eye  disease. 
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Figure  Captions 

Fig.  I.  PSFs  verses  pupil  size  for  a  single  subject.  The  PSFs  are  calculated  based  on  the  subject's  aberrations  measured  with 
the  Hartmann-Shack  wavefront  sensor. 

Fig.  2.  Zernike  description  of  the  eye’s  aberrations.  The  upper  curve  (square  symbols)  shows  the  RMS  wavefront  error  of 
each  Zernike  order  for  a  7.3  mm  pupil  averaged  across  14  human  eyes.  Error  bars  indicate  the  standard  deviation  among  eyes. 
For  the  2nd  order  Zernike  modes,  only  astigmatism  is  shown.  The  average  amount  of  astigmatism  in  these  observers  was  0.6 
diopters,  corresponding  to  a  mean  RMS  value  of  0.77  mm.  The  middle  curve  (triangular  symbols)  shows  the  data  for  a  3.4 
mm  pupil  averaged  across  12  eyes.  The  lower  curve  (round  symbols)  is  for  an  artificial  eye.  The  error  bars  shows  the  standard 
deviation  of  10  repeated  measurements. 

Fig.  3.  System  for  vision  testing  or  retinal  imaging  with  adaptive  optics. 

Fig.‘4.  Contrast  sensitivity  plots  for  two  subjects  before  and  after  adaptive  compensation. 

Fig.  5.  Images  before  and  after  adaptive  compensation  for  aberrations  for  JWs  right  eye.  Both  images  are  of  the  same  retina 
area  at  1  deg  from  the  central  fovea.  Images  were  taken  with  550  nm  light  (25  nm  bandwidth)  through  a  6  mm  pupil.  Each 
image  has  been  normalized  to  its  mean  intensity  and  both  are  displayed  within  the  same  gray  scale  range.  The  dark  vertical 
band  across  each  image  is  an  out-of-focus  shadow  of  a  blood  vessel. 

Fig.  6.  1  deg  images  of  IP's  right  eye  of  the  same  retinal  location  at  two  focal  planes.  In  the  left  image,  capillaries  on  the 
order  of  10  pm  are  seen.  By  focusing  deeper  into  the  retina,  the  underlying  photoreceptor  mosaic  can  be  seen  and  the 
capillaries  appear  as  faint  shadows. 
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Abstract:  Inexpensive  micromachined  varifocal 
mirrors,  providing  very  wide  range  of  focal  distance 
control  are  described.  The  focal  distance  of  a  mirror 
with  a  clear  aperture  of  1cm  can  be  controlled  in  the 
range  oo  to  +0.25m  (four  diopters)  in  a  frequency 
band  of  up  to  75Hz  by  applying  control  voltages  in 
the  range  of  0  to  300V.  Applications  for  focus  con¬ 
trol  in  a  CCD  camera  and  for  direct  correction  of  the 
accommodation  depth  of  a  human  eye  are  reported. 


Al -coated 
nitride  membrane 


Figure  1:  Scheme  of  micromachined  adaptive  mirror. 

Fast  focus  adjustment  is  desirable  in  many  imag¬ 
ing  optical  systems.  Adaptive  mirrors  [1,  2],  used 
for  correction  of  aberrations  with  low  amplitude  are 
in  general  much  faster  than  existing  mechanical  aut¬ 
ofocus  systems  used  for  deep  focus  adjustment.  Re¬ 
cently  developed  silicon  micromachined  adaptive  mir¬ 
rors  [3,  4]  provide  very  wide  range  of  surface  deforma¬ 
tions,  making  possible  direct  focus  control  in  imaging 
optical  systems.  Incorporation  of  an  adaptive  mirror 
into  imaging  lens  system  introduces  a  certain  design 
problem,  because  an  additional  space  should  be  pro¬ 
vided  to  separate  the  incident  and  the  reflected  light 
paths  and  because  the  oblique  incidence  onto  the  de¬ 
formed  adaptive  mirror  results  in  additional  aberra¬ 
tion.  Here  we  consider  two  configmations  which  use 
varifocal  adaptive  mirrors  for  quick  adjustment  of  the 
distance  to  the  sharply  imaged  plane  in  the  range 
oo . . .  0.3m  with  acceptable  optical  quality. 

The  schematic  of  micromachined  membrane  adap- 


Figure  2;  Micromachined  defocus  corrector  mounted 
onto  PCB. 


tive  mirror  [5]  with  a  wide  range  of  focal  distance  con¬ 
trol  is  shown  in  Fig.  1.  The  mirror  is  fabricated  using 
bulk  silicon  micromachining.  The  reflective  surface 
is  formed  by  a  circular  0.5  .. .  1/im-thick  silicon-reach 
nitride  membrane,  which  is  coated  by  a  reflective  alu¬ 
minum  layer.  The  chip  with  circular  stretched  mem¬ 
brane  is  mounted  75pm  over  metalized  printed  circuit 
board  (PCB).  The  radius  of  the  membrane  curvature 
is  controlled  in  the  range  oo . . .  0.5m  by  applying  volt¬ 
age  in  the  range  0 . . .  160V  between  the  membrane 
and  the  metalization  layer  of  the  PCB.  The  com¬ 
plete  device  measures  20x30x5mm  and  weights  ap¬ 
proximately  5g  -see  Fig.  2.  Interferometric  patterns 
correspondent  to  the  initial  and  deformed  surface  of 
the  mirror  are  shown  in  Fig.  3. 

The  dynamics  of  the  mirror  response  depends  onto 
the  resonant  frequency  of  the  membrane  and  the  air 
damping.  Under  normal  atmospheric  pressure  the 
mirror  is  essentially  over-damped,  demonstrating  lin¬ 
ear  frequency  response  in  the  range  0 . . .  75Hz  (the 
particular  mirror  used  in  experiments).  Cutoff  fre¬ 
quency  can  be  extended  into  kHz  range  by  increasing 
the  membrane  tension  during  the  membrane  fabrica¬ 
tion,  which  will  also  rise  the  necessary  control  volt¬ 
age.  To  measure  the  dynamics  of  mirror  response,  the 
variation  of  local  near-field  intensity  of  a  collimated 
light  beam  was  registered  after  reflection  from  the  dy- 
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Figure  3;  Interferogram  of  initial  mirror  surface  with 
active  diameter  of  10mm,  the  same  under  control 
voltage  of  90V,  corresponding  to  the  focal  distance 
of  75cm. 


Figure  4:  Response  of  the  mirror  to  25Hz  meander 
input  signal. 


namically  deformed  mirror,  using  method  similar  to 
one  described  in  [3].  The  oscillogram,  correspondent 
to  the  mirror  response  to  the  rectangular  input  sig¬ 
nal  with  amplitude  150V  and  frequency  25Hz,  corre¬ 
sponding  to  switching  the  focus  distance  in  the  range 
oo . . .  0.3m  50  times  per  second,  is  shown  in  Fig.  4. 
The  device  is  sufficiently  fast  to  provide  independent 
focus  adjustment  for  each  frame  with  a  standard  rate 
of  25  frames  per  second.  Micromachined  adaptive 
mirrors  are  reliable  -  similar  devices  have  been  contin¬ 
uously  driven  during  six  months  with  a  frequency  of 
lOOHz  and  the  range  of  deflection  of  10/zm  (~  1.5-10® 
deformation  cycles)  without  any  observable  change  of 
optical,  mechanical  or  electrical  characteristics. 

Optical  setups  used  for  correction  of  image  fo¬ 
cus  are  shown  in  Fig.  5.  In  the  first  configura¬ 
tion,  the  adaptive  mirror  replaces  the  reflective  coat¬ 
ing  on  one  of  pentaprism  surfaces.  The  angle  of 
axial  incidence  onto  the  mirror  surface  equals  to 
if  =  arcsin(n  ■  sin  22.  5°)~  34°,  where  n  is  the  refrac¬ 
tion  coefficient  of  the  prism  material,  so  the  deformed 
mirror  introduces  considerable  amount  of  astigma¬ 
tism  into  the  wavefront  transmitted  by  the  system. 
The  angular  field  is  limited  to  approximately  24°  by 
the  prism  system  [6]  and  the  image  is  inverted.  De¬ 
scribed  system  was  used  in  combination  with  a  CCD 
camera,  equipped  with  50mm  Fl/5  lens  focused  to 
infinity.  The  system  was  quickly  and  efficiently  fo¬ 


cused  using  the  adaptive  mirror  in  the  range  from  in¬ 
finity  (zero  control  voltage)  to  30cm  (control  voltage 
of  150V).  Example  of  fast  focus  switching  (switch¬ 
ing  time  approximately  10ms,  ‘still  exposition  30ms) 
between  two  objects  placed  at  distances  75cm  and 
50cm  respectively  is  shown  in  Fig.  6.  The  observed 
image  quality  at  distances  shorter  Im  was  slightly 
lower  than  that  obtained  by  direct  manual  focusing 
of  the  camera  lens,  but  still  acceptable  as  all  impor¬ 
tant  details  were  resolved.  Focusing  in  the  range  Im 
to  00  produced  images  with  a  quality  indistinguish¬ 
able  from  direct  lens  focusing. 


Plane 


Figure  5:  Two  configurations  used  for  focus  adjust¬ 
ment. 

The  second  setup,  used  for  direct  correction  of  the 
accommodation  depth  of  the  human  eye,  consists  of 
two  mirrors  in  periscope-like  configuration  -see  Fig.  5. 
Due  to  the  large  angle  of  incidence  (~  45°),  the 
deformed  adaptive  mirror  introduces  a  considerable 
amount  of  astigmatism  when  the  eye  is  accommo¬ 
dated  to  infinity  and  the  system  is  focused  to  a  short 
distance.  Ray-tracing  analysis  showed  that  the  ax¬ 
ial  resolution  of  the  system  (eye  lens  was  considered 
ideal  with  a  focal  distance  of  17mm  and  diameter 
5mm)  drops  down  to  15  lines  per  millimeter  in  the 
case  of  focusing  to  a  distance  of  0.3m.  Compensation 
of  astigmatism  introduced  to  the  model  improved  the 
resolution  to  40  lines/mm.  Visual  experiments,  per¬ 
formed  by  the  author  and  a  group  of  colleagues,  aged 
21  to  65,  demonstrated  practical  efficiency  of  the  re¬ 
ported  scheme.  In  the  experiment,  two  text  patterns 
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Figure  6:  Two  consequent  frames  (time  separation 
40ms),  obtained  by  focusing  the  system  to  the  dis¬ 
tance  of  75cm  (top)  and  50cm  (bottom).  Both  images 
are  mirrored  to  improve  readability. 


were  placed  in  front  of  the  system:  one  at  the  distance 
of  four  meters  (poster,  56pt  Helvetica)  and  the  sec¬ 
ond  (textbook  lOpt  Times-Roman)  at  the  distance  of 
approximately  15cm.  When  no  voltage  was  applied 
to  the  adaptive  mirror,  the  poster  was  clearly  read¬ 
able  by  all  participants  of  the  experiments,  the  text  in 
the  near  zone  was  readable  by  approximately  40%  of 
all  participants  with  a  considerable  eye  effort.  After 
applying  the  voltage  of  150V  to  the  adaptive  mirror, 
all  participants  could  read  the  textbook  without  any 
considerable  eye-strain,  but  the  poster  was  absolutely 
unreadable  for  all  participants.  All  participants  char¬ 
acterized  the  quality  of  the  corrected  image  as  accept¬ 
able,  with  slightly  distorted  field  of  view.  The  cor¬ 
rection  setup  for  one  eye  provides  field  of  view  of  ap¬ 
proximately  30°,  measures  15x25x35mm  and  weights 
less  than  40g.  The  deformable  mirror  is  connected  to 
the  voltage  source  via  resistor  of  5  MOhm  for  safety 
reasons. 

The  optical  quality  of  both  systems  can  be  further 
improved  by  using  astigmatic  varifocal  mirror  to  re¬ 
duce  the  aberration  introduced  by  the  oblique  inci¬ 
dence  onto  the  surface  of  the  deformed  mirror.  Ellip¬ 
tical  membrane  will  provide  necessary  alteration  of 


the  optical  power  without  introducing  on-axis  astig¬ 
matism  at  oblique  angles  in  the  whole  range  of  fo¬ 
cus  control.  More  complex  aberrations  can  be  cor¬ 
rected  using  multi-electrode  control  structures.  Both 
reported  setups  were  not  optimized  with  respect  to 
resolution,  aberrations  and  field  of  view.  Special  de¬ 
sign  of  a  system  including  the  varifocal  mirror  will 
certainly  provide  much  higher  optical  performance. 

Reported  experiments  show  that  micromachined 
varifocal  mirrors  can  be  used  for  fast  focus  adjust¬ 
ment  in  imaging  systems  like  surveillance  cameras, 
industrial  vision  systems,  binoculars,  photo-cameras, 
camcoders,  eye-testing  equipment  and  devices  for  di¬ 
rect  correction  of  vision  acuity. 

This  work  has  been  supported  by  the  long-term 
EC  ESPRIT  project  No.  21063,  “Micro-optical  sili¬ 
con  systems”.  Author  is  grateful  to  P.M.  Saxro  from 
DIMES  for  support  of  this  work. 
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Introduction 

In  vivo  tissue  optical  properties  control  is  very  important  for  many  applications.  A  number  of  laser  surgery,  therapy  and 
diagnostic  technologies  include  tissue  compression  and  stretching  used  for  better  transportation  of  laser  beam  to  underlying 
la\  ers  of  tissue.  The  human  eye  compression  technique  allows  to  perform  transscleral  laser  coagulation  of  the  ciliaiy  body  and 
retina/choroid.  The  possibility  of  selective  clearance  of  the  upper  tissue  layers  should  be  veiy  useful  for  developing  of  the  eye 
ball  diaphanoscopy  techniques  and  for  detecting  of  local  inhomogeneities  hidden  by  a  highly  scattering  medium  in  functional 
tomography.  In  this  report  we  present  results  on  the  human  sclera  optical  and  temporal  properties  controlled  by  employing 
administration  of  osmotically  active  chemicals.  These  results  should  be  very'  useful  for  developing  of  computer-aided  eye 
tissues  photocoagulation  systems  as  well  as  eye  ball  tomographic  systems. 


Biotissues  optical  characteristics  control 
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Scleral  tissue  model 


The  human  sclera  is  a  turbid,  nontransparent  medium  covers  about  80%  of  the  eye  ball  and  serves  as  a  protective 
membrance.  It  is  a  strong  fibrous  tissue  which  is  mainly  consists  of  cojunctive  collagen  fibers  packed  in  lamellar  bundles 
which  are  immersed  within  amoiphous  ground  (interstitial)  substance.  For  normal  sclera  tissue  tyj  «  p,.  For  example,  for 
porcine  sclera  at  633  nm  pa  =  0-2  mm'',  Ps  s  80  mm  ',  g  s  0.9  (M. Hammer  et  al.,  Phys. Med. Biol.  Vol.40,  pp.963-978,  1995). 


Refrective  indices  matching  -  optical  enlightening  effect 

The  ratio  of  the  scattering  coefficients  for  different  degree  of  collagen  fibrils  and  interstitial  substance  refractive  indices 
matching  has  a  view 


n|(t) 
n,(0) 

where  ni  are  the  refractive  indices  of  collagen  and  interstitial  substance,  respectively; 

n|(t)  =  n,(0)C|(t)  +  n.;C^(t), 

C|(t)  +  Cjt)  =  l, 

C|  and  c^  are  relative  volume  concentrations  of  initial  interstitial  medium  and  administrated  chemical  agent  (osmoKle). 
respectively;  ni(0)  and  n^  are  the  corresponding  values  of  refractive  indices. 

For  a  highly  matching  what  is  easely  attainable  in  experiments: 

n^/n|(t)  =  1.001.  When  n(,/n,(0)  s  1.1. 

T<.(t)~exp|-[p5(t)  +  p3]i}, 

where  p-j «  Ps(0)- 1  is  the  tissue  thickness. 


The  human  sclera  spectral  characteristics 


Ri,  % 


T>  pical  experimental  spectra  of  the  human  sclera  samples  measured  for  different  time  intervals  of  60%  Trazograph  solution 
administration.  Spectra  numbered  1  were  measured  in  1  min  after  sample  immersing  into  solution,  spectra  2-8  were  measured 
subsequently  at  2  min  intervals.  An  individual  spectrum  measuring  time  scanning  fi-om  higher  to  lower  wavelengths  was  about 
8.^  sec.  a  -  the  collimated  transminance,  Tc.  sample  thickness  -  0.6  mm;  b  -  the  total  transmittance,  T,  sample  thickness  -  0.7 
mm;  c  -  the  diffusion  reflection  R^.  sample  thickness  -  0.7  mm  (heavy  pigmented  tissue). 

Optical  monitoring  of  maner  diffusion  through  the  scleral  tissue 

Basing  on  the  first  Pick’s  law  for  the  large  external  volume  of  a  chemical  agent  we  can  find  that  its  concentration  within  a 
tissue  sample  has  an  exponential  temporal  dependence 

c,m  =c,,,{l-exp(-t/T)}. 

where 

T  =  1-/D . 

I  is  the  tissue  layer  thickness,  c^xi  =  const,  D  is  the  chemical  agent  difiusion  coefficient 


The  time-dependent  collimated  transmisttance  (dots)  of  the  sclera  sample  of  thickness  0.5  mm  measured  at  633  nm 
concuirently  with  60%  Trazograph  administration.  The  solid  line  and  dashed  lines  represent  the  mean  value  and  the  upper  and 
lower  limits  of  T^  obtained  by  a  calculation  of  Tc  using  experimental  data:  Dt=1.4610’’  ±  I.^IO"*  cnr/s  (see  Table,  sample 
#1). 

Mean  values  and  rms  of  diffusion  coefficient  for  Trazograph  Dj  in  the  human  sclera  samples  defined  on  the  base  of 

experimental  data  for  Tc(t) 


Sample 

Thickness  of  sample, 
mm 

Coefficient  of  diffusion, 
cm'/s 

RMS  of  diffusion  coefficient, 
cm'/s 

#1 

0.5 

1.4610- 

1.90- 10-* 

#2 

0.5 

4.001  O'*  ... 

I.LSIO-* 

#3 

0.4 

5.291  O'* 

7.58- 10'^ 

#4 

0.4 

1.20- 10'' 

l.OMO'* 

Efficiency  of  the  bovine  scleral  samples  optical  enlightening 
(Tcmax .  %)  due  to  various  solutions  administration  at  room  temperature 


nm 

Solution 

450 

500 

600 

700 

Trazograph,  60% 

13 

18.7 

70.2 

76.3 

Trazograph.  76% 

39.2 

50.3 

84.1 

99.2 

Glucose.  45% 

21.6 

29.5 

69.4 

73.8 

PEG.  80%  (20000) 

46.4 

58.3 

81.7 

99.2 

Efficiency  of  the  scleral  samples  optical  enlightening  (Tcmax )  and  time  of  enlightening  (t„,i^) 
via  temperature  of  PEG  (20000)  80%  solution  at  A.=700  nm 


Solution 
temperamre.  "C 

18.5 

20 

m 

30 

33 

38 

40 

43 

46.5 

Tcmix-  % 

81.0 

82.2 

83.5 

85.4 

87.2 

89.6 

90.1 

94.8 

98.5 

tmax-  min 

27 

25 

23 

21 

18 

12 

10 

.  8 
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Estimated  coefficient  of  diffusion  for  difffent  chemical  agents  in  the  bovine  sclera  samples 

Trazogiaph  (76%) 


Sample 

Coefficient  of  diffusion. 
cmVs 

RMS  of  diffusion  coefficient,  cm‘/s 

#3 

1.1810"' 

1.20- 10-' 

#4 

6.28- 1  O'-' 

1.50- 10'' 

#5 

9.08- 1  O'-' 

1.85- 10'' 

#6 

1.2610"' 

2.67- 10-' 

#7 

5.88- 10'' 

7.07- 10-* 

Mean  value 

9. 13' 10-' 

2.75-1  O'-' 

Glucose  (45%) 


Sample 

Coefficient  of  diffusion, 
cm"/s 

RMS  of  diffusion  coefficient,  cm"/s 

m 

6.23-1  O'-' 

7.48-10'* 

#9 

8.29-10-' 

2.45-10-' 

#10 

1.46-10"' 

1.20-10'' 

#11 

4.37-10"' 

2.29-10"' 

#12 

3.38-10"' 

2.32-10'' 

Mean  value 

2.13-10"' 

1.48-10"' 

PEG  (80%) 


Sample 

Coefficient  of  diffusion. 
cm‘/s 

RMS  of  diffusion  coefficient,  cm‘/s 

#13 

1.64-10"' 

2.71-10'' 

#14 

3.92-10"' 

1.10-10"' 

#15 

1.52-10"' 

7.56-10' 

#16 

5.54-10'' 

1.23-10' 

Mean  value 

1.91-10"' 

1.24-10"' 

4 
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Optical  scheme  of  the  scanning  polarization  sensitive  spatial  speckle  correlometer.  1  -  monomode  He-Ne  laser  (633  nm.  2  mW) 
w  ith  linear  polarization  of  its  radiation;  2  -  telescope;  3  -  focusing  micro-objective;  4  -  tissue  under  study;  5  -  computer 
controlled  2D-scanning  device;  6  -  manually  rotating  polarizer;  7  -  photomultiplier  with  a  pinhole. 

Because  in  the  correlation  experiments  only  the  fluemating  component  ^(t)  =  Is(t)  -  <Is(t)>  of  the  scattered  light  intensity  is 
usualK’  determined,  the  corresponding  autocorrelation  function  can  be  obtained  as; 

g.(T)  =<  '^(t)  ^(t  +  !)>/<  X'(t)  > . 

In  the  region  of  the  small  temporal  scales  the  anah  sis  of  the  intensity  fluemations  can  be  more  precisly  carried  out  by  the  smdy 
of  the  asymptotic  behavior  of  corresponding  structure  functions.Difx) 

D,(T)  =  <[I(t+x)-I(t)]->. 

Generalized  parameter  which  can  be  used  for  characterizing  of  this  behavior  is  the  exponential  factor 

^  _ln(D,(T,)/D,(T,)) 

'  Infti/T,) 

where  values  of  T|  and  determine  the  boundaries  of  the  region  of  temporal  scales  for  which  V|  is  estimated. 


CCD  images  of  the  far-zone  speckle  patterns  for  early  (~  2.5  min)  (a)  and  later  (~  10  min)  (b)  stages  of  the  sclera 
enlightenment  using  Trazograph  (60%)  administration. 


0  100  200  300  400  x.fis 

The  human  sclera  samples  autocorrelation  functions  of  the  far-zone  speckle  intensity  flucmations  measured  using  scanning 
spatial  speckle  correlometer.  a  -  the  large  time  scale,  sample  partly  enlightened  by  Trazograph  (60%);  b  -  the  time  evolution  of 
the  fonn  of  autocorrelation  peak  (the  small  time  scale)  with  the  sclera  turbidity  decay,  sample  thickness  0.55  mm, 
measurements  were  done  for  different  time  intervals  of  Trazograph  (60%)  administration:  •  -  120  s;  +  -  220  s;  *  -  320  s;  □  - 
420  s;  X  -  520  s;  0  -  620  s;  A  -  720  s;  S  -  820  s. 


0.0  6.0  12.0  18.0  24.0  At,mm 


lA  pical  dependencies  of  exponential  factor  V|  (a)  and  normalized  value  of  the  average  intensity  <Is>/Io  (b)  measured  in  the 
observ'ation  point  on  time  period  of  Trazograph  administration. 


Ol- 

1.0 


9.4 


17.S 


26.2 


t,  iiiin 


Transmission,  % 


Transmission,  “t 


TTie  time-dependent  averaged  over  the  scanning  trace  (1.5  mm)  values  of  the  mean  speckle  intensity  <!,>  (curves  3)  and  its 
polarization  components  <Ij.>  (curves  1)  and  <I||>  (curves  2)  measured  in  the  paraxial  region  for  the  human  sclera  sample  of 
thickness  0.4  mm:  a.  b  and  c  -  are  the  subsequent  measurements  for  sample  kept  at  first  in  Trazograph  (60%)  solution  (a),  then 
in  physiological  solution  (0.9%  NaCl)  (b).  and  finally  again  in  Trazograph  (60%)  solution  (c);  X.=633  nm. 

•  Conclusion 

The  results  of  this  study  show  that  administration  of  osmolytes  to  the  sclera  effects  the  refractive  index  matching  of  the 
cnllagen  fibrils  and  interstitial  fluid  what  leads  to  dramatic  changing  (reducing)  of  the  scattering  properties  of  sclera. 

Dynamics  of  the  scleral  tissue  optical  enlightening  using  osmolytes  is  defined  by  the  characteristic  time  response  about  3- 
1 0  min. 

Presented  results  are  general  and  can  be  applicable  for  description  of  many  other  fibrous  tissues. 

The  wide  range  of  the  sclera  scattering  parameters  changing  allows  to  understand  and  to  describe  the  transition  from 
diffuse  mode  of  light  scattering  to  a  coherent  light  propagation  through  tissue. 

The  intensity  polarization  measurements  should  be  very  helpful  in  study  of  tissue  structure  changing.  The  optical 
enlightening  of  fibrous  tissue  technology  using  the  suitable  chemical  agent  should  be  very  important  for  developing  of 
computer-aided  adaptive  laser  systems  for  coagulation  of  tissue  disorders  hidden  under  the  high  scattering  tissue  layers. 
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Overview 


•  We  have  used  a  127-element  liquid  crystal  Spatial 
Light  Modulator  (SLM)  to  produce  phase  screens  that 
simulate  atmospheric  turbulence. 

•  We  have  Introduced  the  turbulence  as  a  disturbance 
to  a  Shack-Hartmann  wavefront  sensor  (WFS). 

•  The  error  signal  from  the  WFS  has  been  used  to 
drive  a  deformable  mirror  so  as  to  correct  the 
disturbance. 


The  performance  of  the  adaptive  optical  system  is 
demonstrated  by  measuring  open  and  closed  loop 


Strehl  ratios  for  various  turbulence  strengths. 


iix 


Background 


•  Liquid  Crystal  (LC)  phase  modulators  or  retarders  can 
induce  controlled,  repeatable,  dynamic  aberrations  into 
optical  systems  at  low  cost,  low  complexity,  and  high 
flexibility. 

•  Used  as  wave  front  correctors  they  have  several 
advantages 

-  low  power  density  (1 - 1 0'^  W/cm^) 

-  wide  range  of  phase  delay  (up  to  6  mm) 

-  high  transparency  (>  70%)  for  0.5  <  >.  <  6  pm 

•  LC’s  show  promise  for  atmospheric  wave  front  correction, 
but  as  yet  are  limited  to  low  temporal  and  spatial 
frequencies.  Using  the  device  for  simulating  atmospheric 
turbulence  is  an  important  demonstration  of  the  limitations 
and  potentialities  of  such  modulators. 


BH888888B888a 
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127  hexagonal  elements 
13x13  array 
D  =  1.1  cm 
1  full  wave  of  phase 
retardance  at  ^=670  nm 

0  -  1 0V  AC  operation 

Meadowlark  controller 
operates  from  PC 

Hollow  dots  represent 
element  center  locations 
for  modeling 


Liquid  Crystai  Phase  Moduiator 

(LCPM) 


Meadowlark  Optics  Hexl27 

•Element  Area  -  2.81  mm^ 

•1  HeNe  Wave  (0.63  |im)  Phase  Delay 
•Transmission  >95%  @  0.67  )xm 
•Diffraction  Losses  ~  4% 

•Separation  Between  Elements  -  25  pm 
•LC  Thickness  -  2.5  pm 
•Threshold  Voltage  - 1.2  V 
•Control  Voltage  Range  0-10  V 
•Control  Voltage  Frequency  -  2  kHz 


Present  Experimental 
Configuration 


SLM  Disturbance 

Present  Laboratory  Configuration:  DM  Correction 

Shack-Hartmann  WFS 


Moving  Turbulence 


•  We  represent  non-slewing  turbulence  at  any  instant 
as  a  single  phase  screen  equivalent  to  a  stack  of 
“frozen”  virtual  screens,  each  moving  with  wind 
velocity,  v,  in  Its  own  random  direction  perpendicular 
to  the  line  of  sight  (LOS). 

•  For  each  frozen  virtual  screen,  the  phase  fluctuation 
at  position,  r,  and  time,  t,  Is  the  same  as  at  position 
r+vAt  at  time  t+At. 

•  Isotropy  dictates  that  the  temporal  statistics  of 
turbulence  are  a  function  of  the  magnitude  of  v,  and 
not  its  direction.  For  a  non-slewing  LOS, 

V  =  2.34fQro 
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Moving  Turbulence  (Continued) 


•  We  use  an  FFT  approach  to  generate  the  Zernike 
coefficients,  aj(t) : 

-  initialize  a  64-element  complex  array  of  temporal  frequencies 
with  uniformly  distributed  random  variables  of  unity  standard 
deviation, 

-  force  the  negative  frequency  elements  to  be  the  complex 
conjugate  of  the  positive  frequencies,  and  the  DC  term  to  be 
zero, 

-  multiply  by  a  filter  function  proportional  to  (co/v)-''^/® , 

-  take  the  Fourier  transform,  and 

-  multiply  the  resulting  time  sequence  of  numbers  by 

[cADIrJ'T 


BC-905 

-19- 


ptica! 
kclanccs 
» ompany 


Experimental  Setup  for  Atmos¬ 
pheric  Phase  Screen  Simulation 


Diode  laser  wavelength  =  670  nm,  lens  focal  length  =  300  mm,  pinhole  diameter  =  25  pm 
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Far-Field  Intensity  Patterns  for  Time- 
Correlated  Moving  Turbulence 


Sequence  is  horizontal  raster  scan,  top  left  to  bottom  right 
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Dynamic  Characteristics  of  SLM 


•  Transition  time  depends  on  initial  voltage,  step 
amplitude,  and  step  direction 

•  LCs  are  normally  slower  switching  from  high  to  low 
voltage  than  reverse 

•  LCs  are  normally  slower  switching  over  a  small 
positive  step  than  a  large  one 

•  For  our  SLM,  this  means  a  long  (-200  msec)  settling 
time  from  one  arbitrary  phase  screen  to  the  next 

•  Major  improvements  In  switching  times  will  have  to 
be  made  on  these  devices  before  they  are  useful  as 
real-time  atmospheric  turbulence  simulators 
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Liquid  Crystai  Phase  Moduiator/ 
Adaptive-Optics  Experiment 


/ 


127  Element 
>>^LCPM 


Focal  Plane  Camera 
128X128  Si  CCD 
250  FPS 
Low  Noise 

Pupil  Plane  Camera 
64X64  Si  CCD 
1000  FPS 
Low  Noise 


Mirror  I 

Stroke  6  |xm 
163  Active  Actuators 
Hexagonal  Pattern 
Digital  Storage  of  all  Frames 
Image  Stabilizing  Tracker 


Meadowlark  Optics  Hexl27 
LC  Thickness  -  2.5  |im 
Threshold  Voltage  - 1.2  V 
Control  Voltage  Range  0-10  V 
*  Control  Voltage  Frequency  -  2  kHz 
Element  Area  -  2.81  mm^ 


Dynamic  Response  of  SLM 

— iiMiiii  n  il  1 1  I  I  ir  11"'  — "  .-r-v'-.i';;  •..vvxv'-a'mmu-aBmcnaaMMaaBHMaMMBB 

•  We  switched  the  SLM  from  a  flat  phase  screen  to  1 .0  waves  of 
focus  error 

•  1 00  msec  is  required  to  switch 
from  flat  to  1  wave  defocus 

•  400  msec  is  required  to  switch 
from  defocus  back  to  flat 

•  The  long  switching  time  in  the 
“fast”  direction  is  due  to  those 
SLM  elements  which  experience 
a  small  phase  change 


•  These  two  switching  times  suggest  a  slow,  2  Hz  disturbance 
bandwidth 

BC-905 
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Disturbance,  Sensor,  and  Corrector 
Geometric  Reiationships 


•  Hexagonal  SLM,  used  for  disturbance,  introduces 
piston  in  each  element 

•  Hexagonal  lenslet  array  used  for  Hartmann  wavefront 
sensor 

-  magnified  1:1  to  match  SLM  elements 

-  shifted  in  both  x  and  y,  since  Hartmann  WFS  can’t  detect 
piston 

-  each  SLM  actuator  piston  appears  as  slopes  in  three  WFS 
subapertures 

•  Hexagonal  deformable  mirror  with  169  elements  used 
as  corrector 

-  magnified  1:1  to  match  WFS  subapertures 

-  shifted  in  x  and  y  for  a  triangular  geometry 


Simulation  Results 

Cross-Section  of  Shack-Hartmann  Spot 
Two  Segments  Per  SH  Sub-Aperture 


Hartmann  Spot 
Becomes  Bi-Modal  As 
SLM  Segment 
Displacement 
Approaches  7J2 

Centroid 

Measurement  Error 
Goes  to  100  Percent 


Wave-Front  Sensor  Pixels 


Percent  Error  In  Centroid 


Two  Segments 
Across  Shack-Hartmann  / 
Sub-Aperture 


/  Three  Segments 
Across  Shack-Hartmann 
Sub- Aperture 


Two  Segments  Three  Segments 

Across  Shack-Hartmann  Across  Shack-Hartmann 
Sub-Aperture  Sub-Aperture 


SLM  Total  Displacement  in  Waves 


'W 


SLM  Element  and  WFS  Subaperture 
Geometric  Relationship 
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Active  Area  of  SLM 


127  SLM  Elements 
—  148  Subapertures 


Active  Area  of  Mirror 


W  Correcting  the  SLM  Disturbance 


25X\D 


2X\D 


Open 

Loop 


Strehl=0.22 


Strehl=0.30 


Strehl=0.46 


Strehl=0.34 


Strehl=0.65 


Strehl=0.63 


D/ro  =  25 


D/ro  =  12 


D/rn  =  8 


Open  and  Closed  Loop  Comparison 
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Open 


Closed 


Conclusions 


•  SLM  -  Low  Cost  Adaptive  Optics  for  Instrumentation 

•  When  Used  With  Shack-Hartmann  Sensor 

-  High  Order  Shack-Hartmann  Spot  Aberrations 

-  Inter-Segment  Discontinuity  Must  Be  Constrained  <  Ul 

•  Phase  Shifting  Interferometer(PSI) 

-  Measures  Phase  Directly 

-  More  Complicated 

-  May  be  More  Appropriate  0TUr[^ 

•  Future  Work 

-  Laboratory  Experiments  With  PSI 

-  Closed  Loop  Operation  of  SLM  with  DM  as  a  Disturbance  Source 


SLM  Wave-Front  Corrector 
Proposed  Configuration 


Laser 

Source 


Beam 

Splitter 


SLM 

Disturbance 


Mirror 


A  Fast  SLM  for  the  Future 


•  Goal 

-  X/1 0  performance 

-  1  msec  settling  time  to  provide  500  Hz  operation 

•  Solution 


-  Dual-frequency  LC  material 

-  TNE  timing  calibration 

*  200-element  pre-calibrated  look-up  table 

•  1 0  voltage  steps 

•  10  starting  voltages 

•  2  directions 

*  electronics 

•  1 0  |isec  resolution 

•  “smart”  DACs  that  each  access  the  look-up  table  and  switch 
voltages  based  on  the  initial  and  final  phases  between  two  states 

•  30  VAC  driving  capability  with  two  frequencies  on  all  channels 

•  temperature  stabilization  of  SLM 
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The  electrical  and  optical  properties  of 
liquid  devices  for  adaptive  optics. 


Gordon  D.  Love*^  David  F.  Buscher^  Alan  Purvis’ 
David  J.  Robertson^  and  Ray  M.  Sharpies^. 


School  of  Engineering  and  ^Department  of  Physics 
University  of  Durham,  South  Road,  Durham,  DHl  3LE,  U.K. 


Tel:  +44-191-374-7496 
Fax:  +44-191-374-3709 


Email:  g.d.love@durham.ac.uk 


The  use  of  adaptive  optics  (AO)  for  wavefront  compensation  other  than  for  atmospheric 
turbulence  is  an  emerging,  exciting  field.  The  required  system  specifications  will  obviously  be 
dependent  on  the  particular  application,  and  indeed  may  not  be  well-characterized  for  many  of 
the  potential  uses.  Liquid  crystal  spatial  light  modulators  (LC-SLMs)  are  also  an  emerging 
technology  within  the  field  of  AO.  With  this  in  mind  we  present  here  a  summary  of  the 
specifications  of  both  currently  available  and  predicted  LC-SLMs  which  may  be  used  as 
wavefront  correctors,  and  discuss  how  they  may  be  particularly  suited  for  use  in  non- 
astronomical  AO  systems. 

Liquid  crystal  devices  have  generally  been  used  as  amplitude  modulators  in  the  display 
industry,  but  they  may  equally  well  be  used  as  phase  modulators.  In  general  there  are  two 
classes  of  liquid  crystals  of  interest,  namely  nematics  and  ferroelectrics.  Here  we  concentrate 
on  describing  nematic  devices,  because  currently  they  are  more  suitable  for  analog  phase 
modulation  (this  state  of  affairs  may  change  in  the  future  as  fast  analog  ferroelectric  devices 
become  available).  The  available  devices  can  also  be  sub-divided  into  electrically  addressed  and 
optically  addressed.  We  describe  the  electrically  addressed  devices,  simply  because  they  can  be 
thought  of  as  analogous  to  a  deformable  mirror  in  a  conventional  AO  system.  Other  authors 
have  described  work  on  all-optical  AO  systems  (e.g.  the  review  by  Pepper*). 

Table  1.  shows  a  summary  of  the  specifications  of  nematic  LC-SLMs.  We  describe  the 
properties  using  usual  AO  terminology,  rather  than  typical  LC  parameters  which  are  generally 
aimed  towards  display  applications.  These  specifications  are  stroke  or  throw,  number  of 
actuators,  actuator  impulse  function,  fill  factor,  optical  quality,  transmittance,  dispersion, 
temperature  effects,  hysteresis,  polarization  effects,  response  time,  electrical  drive  requirements 
and  cost.  We  invite  correspondence  from  workers  who  are  interested  in  parameters  which  we 
have  not  mentioned  here.  We  also  note  that  table  1 .  represents  the  current  state-of-the-art  as 
we  see  it,  of  an  evolving  area. 


Parameter 

Typical  Value 

Notes 

Stroke 
(Phase  Shift) 

Typically  several  visible  waves 
(~2pm) 

Stroke  =  birefiingence  x  cell  thickness 
Typical  birefiingence  =  0.2 

Typical  cell  thickness  =  10pm 

Larger  strokes  are  available  simply  by 
increasing  cell  thickness 

Number  of 
actuators  (pixels) 

127  total  in  currently  available 
devices  specifically  for  AO^. 
LCTVs  have  also  been  used 
for  AO  with  very  large  no’s  of 
pixels. 

Much  larger  arrays  are  available. 

E.g.  128x128  or  256x256  in  currently 
available  VLSI  designs  which  potentially 
are  suitable  for  AO. 

Actuator  impulse 
function 

Piston-only 

Tip-tilt  actuators  have  also  been  made 
which  could  be  adapted  into  arrays^. 

Fill  factor 

>90% 

Optical  Quality 

<X/5  P-V 

Typically  LC  displays  do  not  have  good 
optical  quality.  It  is  possible  to  build  very 
high  quality  devices,  e.g.  X/26  P-V,  X/142 
rms.  (at  0.633  pm)  has  been  measured'*. 

Transmittance 

>90% 

Devices  need  to  be  made  fi'om  index 
matched  materials  with  AR  coatings. 

Note  also  that  LCs  can  be  made  to 
operate  in  reflection  or  transmittance. 

Wavelength 

Bandwidth 

Full  visible  and  most  of  IR 

LCs  transmit  from  the  blue  to  >10pm 
except  for  particular  molecular  absorption 
bands  in  the  IR^. 

Dispersion 

Typically  5%  over  150nm 
bandwidth  in  the  visible 

In  the  IR  the  average  dispersion  is  small, 
except  near  absorption  bands. 

Temperature 

Effects 

Change  in  birefiingence  with 
temperature  is  small 

LCs  are  only  in  the  liquid  crystalline 
state  at  certain  temperatures,  typically 

0  to  100°C. 

Hysteresis 

None 

Polarization 

Effects 

Generally,  LCs  can  only  phase 
modulate  polarized  light 

Can  use  device  with  2  LC  layers^,  or  a 
quarter-wave  plate  in  reflection®  to 
modulate  unpolarized  light. 

Response  time 

Typically  45ms  for  a  IX,  round 
trip  time,  if  total  excursion  is 

Faster  materials  are  being  assessed,  in 
particular  using  the  dual  frequency 

restricted  to  be  1  X.  For  several  effect’’*, 
waves  time  -lOO’s  ms. 

~0  to  lOV  AC  (~10KHz)  LCs  draw  negligible  current  and  hence 

do  not  require  large  power  supplies 

~  5K  to  30K  US$.  Cost  varies  greatly  depending  on  device. 

They  have  the  ability  to  become  low-cost 
if  large  numbers  are  required. 

Table  1.  Summary  of  the  specifications  of  liquid  crystal  devices  for  adaptive  optics. 

The  above  specifications  show  that  LC-SLMs  have  a  number  of  features  which  mean  they  may 
be  particularly  useful  for  industrial  and  medical  applications  of  adaptive  optics.  They  are  low 
cost  in  comparison  to  typical  deformable  mirrors  (although  not  necessarily  in  comparison  to 
new  micro-machined  mirrors).  They  have  been  proven  as  precise  wavefront  controllers  whose 
pixels  can  be  simply  controlled.  They  operate  (in  general)  in  transmission  and  so  can  be 
incorporated  into  optical  systems  simply.  For  astronomical  AO  systems  the  largest  drawback 
of  LCs  has  been  their  relatively  slow  response  time  (although  this  point  is  now  being 
addressed),  which  may  not  be  an  issue  for  many  non-astronomical  applications.  Devices  are 
currently  being  planned  which  have  a  very  large  number  of  actuators  (>  64x64)  pixels,  and  so 
if  one  is  interested  in  a  very  high  degree  of  correction  LCs  may  be  the  only  option  available. 

In  contrast  LCs  cannot  handle  very  large  optical  power  densities  and  so  may  be  unsuitable  for 
certain  laser  applications.  They  may  also  not  be  suitable  where  extreme  temperature 
fluctuations  are  to  be  expected  (>100°C). 


Driving 

requirements 

Cost 
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INTRODUCTION 

The  commercially  available  liquid  crystal  (LC)  light  valves  were  developed  in  the  middle 
of  1970s  and  had  been  extensively  used  for  display  applications  [1],  later  for  laser  optics 
[2,3],  for  radar  applications  [4]  and  even  for  generation  of  femtosecond  optical  pulses  [5]. 
However  in  adaptive  optics  LCs  became  to  use  recently.  Probably,  not  all  resources  of 
deformable  mirrors  were  exhausted.  Now  the  situation  has  changed.  Already  there  are  a 
high  technology  of  manufacturing  of  commercially  available  liquid  crystal  wavefront 
correctors  for  adaptive  optics  specially  [6]  and  hardware  for  their  control.  I  shell  not  repeat 
all  merits  and  demerits  such  devices.  You  can  find  its  in  Ref.  7  and  references  therein.  But, 
as  rule,  no  one  take  attention  features  of  LCs  except  for  a  necessity  of  alternating  control 
voltage  using  to  prevent  the  electrolytic  dissociation  of  LC  molecules.  In  my  report  I  would 
like  to  attract  attention  to  features  of  adaptive  optical  elements  based  on  nematic  LCs  with 
initial  planar  orientation  and  to  regard  use  of  these  features  for  adaptive  optic  applications. 

0.  FEATURES  LC  ADAPTIVE  OPTIC  ELEMENTS 

1.  The  necessity  to  use  of  an  alternating  control  voltage  to  prevent  the  electrolytic 
dissociation  and  polarization  of  LC  molecules. 

2.  Phase  delay  in  the  LC  layer  is  determined  by  the  effective  value  of  the  voltage  modulus. 

3. '  Nonlinear  character  of  the  response  time  of  LC  correctors,  depending  on  the  control 

voltage  and  the  layer  thickness. 

4.  It  is  necessary  to  take  into  consideration  the  spectral  composition  of  the  control  voltage. 

5.  LC  correctors  have  parametric  capacitance. 

6.  Scattering  near  the  threshold  and  at  the  moment  of  inversion  of  dielectric  anisotropy  for 
LCs  with  low-frequency  sign  inversion  of  dielectric  anisotropy. 

7.  Optoelectronic  characteristics  have  strong  temperature  dependency. 

8.  Existence  of  the  threshold  voltage  and  possibility  variation  of  it  by  initial  orientation 
selection . 

1.  ALTERNARING  CONTROL  VOLTAGE 

The  control  birefringence  of  LC  layer  can  be  induced  both  by  static  and  alternating 
electrical  field.  In  the  case  of  static  fields  the  current  must  flow  through  electrodes,  so 
electrodes’  processes  become  important  The  electric  field  in  the  liquid  crystal  layer 
becomes  inhomogeneous  and  is  determined  by  anisotropy  conditions.  Double  layers  are 
formed  at  the  electrodes,  which  decrease  the  field  in  the  cell,  so  a  higher  applied  voltage  is 
needed  to  reorient  the  layer  as  in  the  case  of  alternating  fields.  Development  of  the  double 
layers  takes  several  seconds  [8],  so  they  can  be  completely  neglected  for  alternating  fields. 
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The  use  of  unipolar  voltage  is  no  desirable  too.  The  unipolar  pulse  voltage  causes 
polarization  of  liquid  crystal  and  this  circumstance  increases  the  relaxation  time  after 
switching  off  of  control  voltage  [9]. 


Low  frequency  dual  polarity  control  voltage  can  stimulate  electro-hydrodynamical 
instability  [10]  when  the  frequency  fixl/te  [11]  where  Te»e|/47iCT||  is  the  dielectric 
relaxation  time,  en  is  dielectric  permittivity  along  initial  orientation,  and  ct||  is  specific 
electric  conductivity  measured  in  the  same  direction.  Besides,  conductive  currents  can 
simulate  hysteresis  [12].  This  circumstance  can  be  avoided  by  changing  LCs  with  high  a||. 

In  practice,  in  order  to  avoid  ion  displacements  we  shall  use  an  alternating  bipolar  voltage 
of  frequency  greater  then  the  carrier  relaxation  frequency.  In  the  case  of  alternating  field  we 
assume  the  frequency  to  be  high  enough  that  the  alignment  of  the  molecules  cannot  follow 
the  alternation  of  the  field  (typically>100Hz)  but  low  enough,  that  the  polarization  can 
follow  it  (typically<10^Hz). 


2.  MODULUS  OF  CONTROL  VOLTAGE 

The  orientation  of  nematic  LC  molecules  (the  phase  delay  value)  depends  on  the  effective 
value  of  the  voltage  modulus.  This  allows  to  develop  a  new  type  of  wavefront  correctors 
[13].  I  have  named  theirs  modal  wavefront  correctors  because  theyhave  control  distributed 
in  aperture  like  modal  deformable  mirrors.  The  simplest  modal  LC  corrector  is  a 
cylindrical  adaptive  lens  (Fig.  la).  The  cylindrical  lens  has  the  transparent  iow-Ohmic 
electrode  and  the  transparent  high-Ohmic  electrode  coated  by  two  equidistant  line 
conducting  contacts  over  aperture  ceil  edge.  Electrically  this  element  can  be  represented  by 
means  of  distributed  capacity  and  resistance  (Fig.  lb).  Electrical  voltage  distribution  is 
described  by  the  Laplace  equation: 


dU 

RCiU)-. 


(I) 


In  approximation  RC  =  const,  the  stationary 
distribution  of  voltage  can  be  obtained  for  harmonic 
control  voltages  Ux^U^smim)  and  C/2=C/osin(co/+(p) 
analjdically.  But  now  we  are  developing  the  more 
common  model  which  includes  the  dynamics  of 
modal  LC  correctors.  The  results  of  numerical 
observations  are  shown  in  Fig.  2a,b. 


Fig.  1  Adaptive  cylindrical  lens  -  a)  set-up,  b)  equivalent  electrical  scheme:  (1)  glass  substrates,  (2)  oriented 
covers,  (3)  transparent  low-Ohmic  electrode,  (4)  transparent  high-Ohmic  electrode,  (5)  equidistant  conductive 
contacts,  (6)  LC  layer,  (7)  spacers.  (Ui  and  Ui)  control  voltages. 
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Fig.  2.  Dynamics  of  voltage  variation  in  cylindrical  adaptive  lens. 


Fig.  3.  Phase  replay  of  LC  layer. 

I  have  cut  off  transition  process  from  switching  on  moment  of  control  voltage  till  S-RC 
where  5  is  square  of  the  LC  cell.  I  have  kept  one  period  of  voltage  variation  and  transition 
process  of  voltage  switching  off  This  will  be  important  in  our  subsequent  consideration. 
LC  phase  replays  at  shown  in  Fig.  2a,b  voltage  distributions  are  represented  in  Fig.  3a,b, 
accordingly.  It  is  necessary  to  say  that  these  figures  were  calculated  without  taking  account 
for  electrical  transition  process.  This  is  very  rough  approximation  but  temporal  scales  in 


Fig.  2  and  Fig.  3  differ  in  one  thousand  times.  The  processes  of  switching  on  and  switching 
off  are  shown  separately.  The  phase  delays  have  cylindrical  shapes  and  LC  cell  works  like 
a  cylindrical  lens.  The  shape  of  phase  delay  distribution  depends  on  RC  and  feeding 
voltage  parameters. 


Modal  LC  correctors  can  be  control  by  variation  of  feeding  voltage  spectrum.  For  example, 
the  rectangular  and  serrated  dual  polarity  feeding  voltage  give  different  phase  distributions 
(Fig.  4a, b).  Represented  results  have  been  obtained  for  first  ten  harmonics  and 
experimental  dependency  AO(L).  Searching  of  the  controlling  voltage  shape  at  the 
desirable  profile  of  phase  delay  distribution  is  enough  challenge.  In  the  simplest  case  we 
can  attempt  to  find  the  feeding  voltage  shape  to  minimize  the  mean  square  deviation  from 
parabola  phase  delay  distribution  [14].  The  phase  delay  distribution,  voltage  modulus  one 
and  the  intensity  distribution  in  focal  plane  are  shown  in  Fig.  5  for  the  ideal  LC  cylindrical 
lens  (Fig.  5a),  for  real  lens  with  simple  sine  feeding  voltage  (Fig.  5b)  and  with  feeding 
voltage  of  complex  shape  (Fig.  5c).  This  shape  is  represented  in  (Fig.  5d). 


Focus  distance  F=5tn,  rms  from  parabola 
5=0.89471,  Strehl=11.8%. 


b) 


Fig.  4.  Control  voltage  shape,  modulus  of  them  distribution  (dotted  line)  and  LC  phase  replay. 

The  adaptive  spherical  lens  has  the  same  design  but  the  conducting  contacts  should  be 
made  in  the  circle  form. 

The  another  more  complex  modal  correctors  is  shown  in  Fig.  6.  It  has  16x16  electrodes 
with  0.3mm  gaps  between  them.  The  resistance  between  neighbours  electrodes  is  20kOm, 
and  the  LC  thickness  is  12p.  From  Eq.l  at  boundary  condition  lj\^=UQsm(Git)  and  L4.„=0 
we  can  obtain  the  voltage  distribution  for  any  separate  electrode: 

U(x,t)  =  Uf,  exp(-ax)  cos(©f  -  or ) ,  (2) 

where  a  =  ^^aRC/2  .  By  increasing  feeding  voltage  frequency  the  mutual  influence 
between  electrodes  is  decreasing.  This  is  illustrated  by  means  interferograms  in  Fig.  7. 
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Thus,  to  suppress  wavefront  aberrations  by  means  such  corrector  in  the  first  place  it  is 
necessary  to  introduce  low  frequency  control  voltages  for  low-order  aberrations,  and  then 
to  add  high  frequency  control  voltages. 


u-l - ! - 1 - 1 - 1  T' 

0.0  0.2  0.4  0.6  X  cm  1.0 

F=4.95m.  5=07t.  Strehl=100%. 


C) 


F=4.95m,  5=0.089Tt,  Strehl=83.8%. 


b) 


V”! - 1 - 1 - 1 - 1 - T’ 

0.0  •  0.2  0.4  0.6  X,  cm  1.0 

F=4.95m.  5=1.0047:,  Strehl=6.47%. 


;  U.  in  per-unit  t 


Fig.  5.  Selection  of  optimal  control  voltage:  a)- 
c)  phase  delay  and  voltage  modulus  distributions 
(dotted  line)  and  intensity  distribution  in  focal 
plane  (upper  part  of  drawing);  d)  shape  of 
optimal  control  voltage.  The  real  interferogram 
is  placed,  in  the  right  lower  comer  of  b) 


Fig.  6.  Many  electrodes  modal  corrector:  (1)  LC  layer;  (2)  glass  substrate;  (3) 
transparent  low-Ohmic  layer;  (4)  high-Ohmic  layer;  (5)  control  electrodes;  (6) 
spacers;  (7)  oriented  covers;  (8)  dielectric  mirror. 


3.  RESPONCE  TIME  OF  LC  CORRECTORS 

The  rapidity  of  LC  correctors  is  determined  by  the  liquid  crystal  layer  thickness  as  ~l/i^ 
therefore  thin  cells  are  preferable  for  fast  phase  modulation.  However,  there  are  two 
problems.  The  first  is  a  small  depth  of  modulation  which  is  more  limited  then  given  by 
theory  [11].  The  second  is  high  probability  of  breakdown  or  short  circuit.  There  is  a  critical 
value  60V/p  [9]  of  the  strong  electrical  field.  In  the  real  substrates  the  liquid  crystal 
molecules  surface  alignment  is  seldom  uniaxial  due  to  the  surface  roughness,  defects  or 
irregular  alignment  treatments.  The  "virtually  inactive"  thickness  is  defined  by  the  quality 
and  preparation  of  the  substrates.  Under  the  conditions  of  strong  anchoring  and  good 
surface  quality,  the  thickness  of  the  partially  disordered  surface  layers  is  reduced.  On  the 
other  hand,  the  large  "virtually  inactive"  thickness  is  expected  for  weak  anchoring  and 
rough  surface  alignment.  In  Ref  15  "virtually  inactive"  thickness  is  estimated  as  ~0.3)i  for 
each  surface,  and  in  Ref  16  an  estimate  of  0.5p  for  the  thickness  of  the  inactive  volume  per 
cell  is  thus  obtained. 

The  another  important  parameter  for  rapidity  is  amplitude  of  control  voltage.  Well  known 
that  the  rise  time  decreased  on  increase  in  the  control  voltage  as  and  the  decay  time  do 
not  depend  from  control  voltage. 


Adaptive  systems  with  simple  amplitude  control 
have  not  high  rapidity.  The  typical  characteristic  of 
such  system  is  shown  in  .  " 

adaptive  system  represented  stabilized  interferometer 
operating  on  the  basis  of  the  aperture  probing 
principle  [17] 

The  dynamic  control  can  be  realized  by  varying 
switching  on  and  switching  off  time  at  fixed  voltage 
amplitude.  Such  control  can  be  used  for  modal 
correctors.  For  example,  the  desirable  focus  distance 
can  be  achieved  by  this  method.  Fig.  9  shows  focus 
variation  of  cylindrical  lens  after  switching  on  and 
off  of  control  voltage. 

1000^  F,  m 

:  OFF  /'■' 
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Fig.  8.  External  perturbations  suppression 
coefficient  Q  in  system  with  aperture 
probing,  f  is  frequency  of  external 
perturbations. 
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Fig.  9.  Focus  variation  according  to 
Fig.  3a. 


This  method  was  modified  for  LCs  with  a  low-firequency 
sign  inversion  of  dielectric  anisotropy.  In  this  case  the  fotf 
value  is  comparable  with  the  ton  and  the  first  is 
determined  by  applied  voltage  as  well.  To  decrease  the 
phase  delay  by  the  quantity  AO  it  is  necessary  to  apply  a 
low-frequency  voltage  at  the  time  rs=AO/Os  where  Os  is 
the  S-effect  rate.  In  order  to  decrease  the  phase  delay  by 
the  same  quantity  it  is  necessary  to  apply  a  high- 
frequency  voltage  at  the  time  /b=AO/Ob  where  Ob  is  the 
B-effect  rate.  The  ratio  between  ts  and  ts  determines  the 


birefnngent  state  of  the  liquid  crystal  layer  in  dynamics.  The  static  value  of  AO  depends  on 
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harmonics  of  control  voltage.  Dual  frequency  control  allows  to  increase  rapidity 
significantly  (Fig.  10)  [18]. 

The  idea  of  transient  nematic  liquid  crystal  effect 

[19]  is  to  utilize  the  fast  decay  time  due  to  the  small 
relaxation  angle  because  of  highly  deformed  liquid 
crystal  directors.  A  relatively  high  ac  voltage  is 
initially  applied  to  the  liquid  crystal  cell.  As  a  result, 
almost  all  the  molecules  are  aligned  by  the  electric 
field  almost  orthogonal  to  the  substrate  surfaces 
except  in  the  boundary  layers.  When  these  highly 
deformed  directors  start  to  relax,  i.e.,  the  voltage  is 
removed  completely,  the  directors  undergo  free 
relaxation.  However,  when  the  phase  delay  reaches 
the  desirable  value  (this  is  possibly  because  of  2n 
phase  degeneration),  a  voltage  is  applied  to  the  cell 
to  stop  the  directors'  motion. 

.■\uthors  of  [19]  make  difference  between  this  mode  of  operation  and  bias  voltage  effect 

[20]  (sometimes  called  surface  mode),  because  bias  voltage  is  constantly  present  during  the 
relaxation  process.  The  presence  of  this  holding  voltage  provides  an  intermediate  stopping 
force  to  the  LC  directors  and.  thus,  decreases  the  relaxation  time.  In  reality  these  effects 
nave  the  same  cause  -  bias  voltage.  In  case  of  transient  effect  the  control  voltage  represents 
the  periodic  or  quasiperiodic  pulse  sequence.  The  periodic  puise  voltage  has  low-frequency 
harmonic  compositions  which  play  a  role  of  ''bias"  voltage.  This  will  be  ciear  from  next 
section. 


Fig.  10.  External  perturbations 
suppression  coefficient  in  phase 
conjugation  in  stabilized  interferometer. 


4.  SPECTRAL  COMPOSITION 


When  we  use  bipolar  rectangular  control  voltage  we  can  drive  a  modulator  by  variation  of 
pulse  period-to-pulse  duration  ratio  q  =  T't;  T  is  period  of  pulse  sequence,  x  is  pulse 
duration.  Amplitude  of  m  harmonics  is  given  by 


.  f  x' 

- sm  mn- 

Txm  \  r. 


3-  A<h, 


2- 


1- 


°  0.1  0.3  0.5  1/q  1.0 

Fig.  12.  Dependence  of  the 
electrooptical  response  of  LC  on 
Mq-,  ?7o=10V. 


(2) 


1.0 


Fig.  11.  Variation  of  harmonics’  modulus  as 
function  of  l/q\  (I)  first  harmonics,  (II)  second 
harmonics,  (III)  third  harmonics. 
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The  given  voltage  represents  a  number  of  harmonics  and  each  harmonics  makes  its 
contribution  in  the  birefringence  magnitude.  Variation  of  parameter  q  leads  to  harmonics’ 
amplitude  change  and,  consequently,  to  variation  of  phase  modulation  depth  in  liquid 
crystal  layer.  For  example,  the  total  effect  from  harmonics  in  cases  "a"  and  "b"  (Fig.  1 1)  is 
different.  The  experimental  dependence  AC){l/q')  is  shown  in  Fig.  12 

5.  PARAMETRIC  CAPACITANCE 

The  dependence  Ad)(l/^)  is  non-monotone.  It  repeats  the  sum  of  first  harmonics’  modulus. 
However,  higher  harmonics  are  shunted  by  a  modulator  capacitance. 

Capacitance  variation  has  the  same  nature  as  birefringence  variation.  This  is  orientation 
change  of  LC  molecules.  But  in  really  this  connection  is  more  complex  because  of 
dielectric  losses  tangent. 

Results  of  measurements  of  capacitance  and  electrooptic  replay  I O I  =A0^3j(-A0  are 
presented  in  Fig.  13.  Here  AO„j^=2£/AnA,  in  tz  unit,  is  LC  layer  thickness,  An  is  optic 
anisotropy,  and  X  is  wavelength.  I  have  used  the  sine  control  voltage.  In  the  case  of  voltage 
increasing  curves  |<b|  (C/)  and  C{U)  are  similar  (Fig.  13a).  For  LCs  with  a  low-frequency 
sign  inversion  of  dielectric  anisotropy  frequency  dependence  is  more  complex  (Fig.  13b). 
Since  capacitance  dependence  correlates  with  the  phase  delay  dependence,  this  complicates 
the  feedback  design  significantly. 

C,  nF  ---  -iOl,- 


£  —  <1 


“  :  :  :  u.  V  3'  “ :  2  -  :  y,  V  '2' 

Fig.  13.  Capacitance  dependence  of  (solid  line)on  applied  voltage  for  LC  cell  with  LC654  (a) 
(/^IkHz)  and  frequency  for  LClOOl  (b)  Electrooptic  relay  is  shown  by  dotted  line. 

6.  SCATTERING 

Two-frequency  addressing  first  have  been  proposed  for  dynamic-scattering-type  liquid 
crystal  displays.  The  cutoff  frequency  marks  the  transition  from  conductivity  anisotropy 
regime  to  a  dielectric  anisotropy  induced  alignment  regime  [23]. 

In  dual-frequency  control  experiment  we  also  observed  the  light  scattering  at  the  instant  of 
control  voltage  frequencies  switching.  However,  the  false  error  signal  was  not  registered 
because  it  had  the  frequency  beyond  the  transmission  bandwidth  of  the  feedback  loop.  The 
scattering  can  limit  the  rapidity  of  adaptive  systems. 
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The  another  scattering  types  are  connected  with  threshold  existence  and  temperature 
fluctuations.  They  are  wrote  in  literature  in  detail  [24, 25]. 

7.  TEMPERATURE  DEPENDENCY 

The  rapidity  and  phase  delay  value  have  strong  dependency  from  temperature.  The  first  is 
connected  with  dynamic  viscosity  yi  which  can  be  decreased  by  increasing  temperature: 
Yi~exp(/4/^7^  where  ^  is  an  activation  energy  [1 1],  and  k  is  the  Boltzmaim  constant.  But  in 
this  case  elastic  constant  Ku  is  decreasing  also:  where  S  is  the  order  parameter  of 

the  liquid  crystal.  The  linear  decrease  of  yi  /  ATn  (this  relation  define  the  rapidity  [22])  is 
observed  when  temperature  is  increased  to  some  optimal  temperature  Topt .  Above  Topt ,  the 
elastic  constant  drops  more  significantly  than  the  effective  viscosity  does,  resulting  in 
slower  decay  time  [26]. 


8.  THRESHOLD  VOLTAGE 

In  wavefront  sensors  based  spatial  filtration  there  are  some  difficulties  associated  with 
formation  and  alignment  of  spatial  filters.  Space  light  modulators  (SLM)  which  operate  the 
light  passing  through  them  can  be  able  to  overcome  these  problem.  When  SLM  operating 
in  self-consistent  filtration  regime  are  used  as  the  filters  then  the  mentioned  difficulties  can 
be  partially  eliminated.  Here  the  filter  is  produced  automatically  in  accordance  with  the 
optical  signal  and  no  alignment  is  needed  when  the  light  beam  incidence  angle  at  the  input 
of  the  optical  system  is  changed.  Existence  of  the  threshold  voltage  can  be  used  to  select 
zero  spatial  component.  We  used  this  possibility  for  visualization  phase  distonions  by  the 
Foucault  knife-edge  method,  by  the  autoreference  interferometer  and  by  the  Zemike  valve 
method  [27],  The  last  was  used  for  compensation  of  weak  phase  distortions  as  well  [28]. 

The  threshold  voltage  can  be  changed  by  initial  orientation  of  LC  because  of  elastic 
constants’  values  depend  on  the  method  of  orientation.  In  the  case  of  planar  orientation 
(i=l,  2,  3)  are  greater  than  in  the  tilted  one,  and  the  decay  time  is  shorter  [9]. 
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Motivation 


Open  loop  adaptive  optics  correction  based  on  turbulence 
prediction  rather  than  closed-loop  measurement  and 
correction  at  all  apermre  points 

-  Measurement  at  a  few  points  predicts  aberration  over  tbe  apermre 

-  Beacon  measurements  are  unnecessary 

Outline 

•  Introduction  to  prediction-based  adaptive  optics 

•  Spatially  and  temporally  structured  turbulence 

•  Optical  correction  experiment  geometry:  beam  jitter  signal 

•  Prediction  method:  LSE  description,  jitter  prediction 

•  Simulation  and  experiment  results 

•  Future  directions 


I  want  to  present  our  experiments  in  adaptive  optics  correction  that  are  based 
on  the  prediction  of  turbulence  rather  than  measurement  of  the  aberration  at 
every  point  in  the  aperture.  In  this  way  only  a  few  measurements  at  the  time  of 
correction  are  necessary  so  that  the  system  could  potentially  operate  at  very 
high  speed.  And  there  is  no  need  for  a  beacon  measurement,  that  is  a 
measurement  of  the  uncorrected  signal  or  of  the  corrected  signal.  The  system 
operates  open  loop. 


This  talk  will  first  introduce  predicton  based  adaptive  optics.  I’ll  show  the 
kinds  of  spatially  and  temporally  structured  turbulence  we  want  to  correct  as 
well  as  the  system  used  to  correct  it  The  method  of  prediction  we  use  is  called 
Linear  Stochastic  estimation.  I’ll  present  the  results  of  our  simulations  and 
experiments  showing  that  the  technique  actually  works  and  then  I’ll  discuss 
some  implications  of  the  work. 


Introduction 


Prediction  -  Based  Optical  Wavefront  Correction 


Linear  Stochastic  Estimation 

•  Correction  of  structured  turbulence  otherwise  not  easily  described  by  statistical  models 

•  Measurement  of  flow  at  a  few  points  predicts  aberradon  at  other  locations  and  time 

•  Open  loop  correcdon  system:  no  beacon  measurement  or  feedback  loops 

•  Lower  operadonal  bandwidth  required 


Highly  structured  turbulence  such  as  aircraft  boundary  and  shear  layers 
requires  very  high  bandwidth  measurement  and  correction  systems  especially 
when  high  spatial  resolution  is  required.  The  highly  organized  features  of  such 
turbulence  and  their  intermittancy  are  difficult  to  describe  using  statistical 
turbulence  models 

We  propose  a  system  by  which  stochastic  estimates  of  flow  properties  are 
developed  from  previous  flow  measurements.  This  information  is  combined 
with  a  few  point  measurements  of  the  flow  at  the  time  of  correction  to  predict 
the  necessary  correction  that  takes  place  at  another  location  in  the  flow. 

So  rather  than  making  measurements  at  aU  points  in  the  aperture  and  driving 
all  apertures  of  an  adaptive  optics  system  based  on  those  measurements  in  a 
closed  loop  system,  we  are  making  open  loop  corrections  based  on  an 
stochastic  estimate  of  the  flow,  imposing  a  condition  on  that  flow  from  a  few 
point  measurements  and  then  inferring  the  optical  aberration  in  the  area  of 
interest.  Shown  here  is  a  point  temperature  measurement  taken  at  the  front  of 
an  aircraft  turret,  used  to  predict  the  tilt  along  a  beam  path  at  the  back  of  the 
turret. 

Key  advantages  of  predictive  systems  are  that  only  a  few  measurements  and 
calculations  need  to  occur  at  the  time  of  correction  lowering  the  operational  bw 
requirements  of  the  system.  And  no  beacon  measurements  of  the  turbulence 
are  required.  We  don’t  measure  the  beam  we  simply  correct  it. 

Potentially  we  may  be  able  to  use  a  few  point  measurements  to  infer  larger 
spatial  areas  usable  for  image  correction.  In  this  case  we  have  pictured  a  point 
measurement  inferring  a  path  integral  for  simple  jitter  correction. 


Here,  I  want  to  illustrate  the  differences  between  our  system  and  a  traditional 
closed  loop  system 

In  the  closed  loop  case,  the  mirror  position  is  measured,  fed  back,  subtracted 
from  the  disturbance  signal  and  residual  signal  is  sensed.  This  is  the  beacon 
signal  I  mentioned,  a  signal  that  is  constantly  monitoring  the  performance  of 
the  adaptive  optics  and  feeding  it  back. 

In  the  predictive  system  a  guess  is  made  of  the  form  of  the  disturbance  based 
on  past  information.  In  this  case,  LSE  is  the  means  by  which  the  past 
information  is  used  to  make  the  prediction.  The  tilt  sensor  is  only  a 
performance  monitor  for  experimental  convenience  in  this  system.  It  does  not 
contribute  to  the  correction  signal. 


Notes: 

PD  Controller=  Proportional  derivative  controller,  produces  weighted  value  of 
the  current  mirror  position  and  the  difference  between  the  current  position  and 
the  last  position. 

Predicts  mirror  movement  ahead  of  time  so  that  mirror  response  dynamics 
(phase  lag  times)  can  be  accounted  for. 

Suggest  the  possibility  that  these  types  of  systems  might  work  together  in  a 
single  adaptive  optics  system.  An  example  will  be  shown  at  the  end. 


Typical  probe  locations  and  how  those  locations  can  be  chosen.  Probe  location 
is  in  general  related  to  the  areas  with  the  highest  frequency  content  and  the 
strongest  signal  fluctuations.  The  RMS  temperature  profile  shows  that  the 
strongest  fluctuations  occur  near  the  geometrical  edge  of  the  flow  defined  by 
the  nozzle  exit  edge.  Typically,  one  probe  is  located  inside  the  edge  and  the 
other  probe  is  located  outside. 


Probe  location  is  not  always  critical  as  long  as  enough  frequency  components 
are  present  to  truly  represent  the  flow  events  that  cause  optical  aberration. 


Steering  Mirror  Control  System 


Uncorrected  beam  jitter  signal 


10  uR  tilt  imparted  by  IJ  *  temperature  variations 


System  diagram 


The  system  is  composed  of  a  He  Ne  laser  and  a  piezoelectric  steering  mirror. 
The  2  mm  wide  beam  propagates  through  a  disturbance,  a  heated  air  flow  in 
this  case,  and  the  position  of  the  beam  is  monitored  using  a  position  sensitive 
detector.  Dynamic  temperature  variations  in  the  air  flow  cause  the  beams 
position  to  wander.  The  piezo  mirror  is  used  to  compensate  for  this  beam  jitter. 
The  mirror  is  driven  by  a  tilt  prediction  produced  by  a  combination  of  a 
statistical  knowledge  of  the  flow  combined  with  2  point  temperature 
measurements  taken  just  before  correction  is  to  occur. 


Piezo  mirror  is  driven  in  1  D  not  by  direct  flow  measurements  but  by  the  jitter 
signal  estimated  from  2  anemometer  probe  temperature  measurements 


Heated  Round  Jet  flow  structure 

•  Kelvin  -  Helmholtz  instabilities  form  highly 
structured  vortices  that  create  optical 
distortion 

•  flow  visualizations 

•  flow  can  be  controlled  using  acoustic 
forcing  from  ordinary  speaker 


Rather  than  use  an  aircraft  geometry  we  have  chosen  a  laboratory  air  flow  that 
provides  interesting  features  and  coherent  strucutures  while  maintaining  some 
level  of  controllability 


I  will  show  a  visualization  of  the  flow  showing  the  kinds  of  structure  that  can 
cause  optical  distortions  in  the  beam. 


In  another  kind  of  visualization,  where  quantitative  temperature  measurements 
are  made  using  a  1 -dimensional  Hartmann  sensor  we  can  see  some  of  the 
different  states  of  the  flow  that  we  can  induce  by  acoustically  forcing  the  flow. 


More  details  about  Hartmann  sensor  measurements  of  the  flow  are  given  in: 

L.  McMackin,  B.  Masson,  N.  Clark,  K  Bishop,  R.  Pierson  and  E.  Chen, 
"Hartmann  wavefront  sensor  smdies  of  dynamic  organized  structure  in  flow 
fields,"  AIAA  Journal  33  (11),  2158-2164,  1995. 


Linear  Stochastic  Estimation 


•Predict  path-averaged  beam  jitter  trom  several  point  temperature 
measurements  near  path 

•Estimate  formed  Iroin  unconditional  correlations  of  the  temperature 
and  the  Jitter. 


-  Take  linear  term  from  power  series 
expansion  of  expected  value  of  jitter,  J^{t) 
given  temperanire  at  positions  ,tj,  T,(t) 

-  minimize  mean  square  error 
between  estimate  and  jitter 

-  results  in  correlation  expression 
solved  for  coefficients,  Aj  “Weights" 


4(f)=(y,(r)|r..(f-T,)) 


•LSE  for  2  reference  points 


These  calculations  are  performed 
before  correction 


■(r,r.)  (r,A)‘ 

■A,' 

.(^a4). 

This  slide  explains  we  really  doing  with  the  temperature  probe  signals: 

What  is  LSE? 

We  are  trying  to  predict  the  jitter  caused  by  temperature  variations  that  occur 
along  the  entire  path  of  the  beam  through  the  flow  from  just  a  couple  of  point 
temperature  measurements.  What  makes  this  kind  of  estimate  very  useful  is 
that  it  is  produced  from  unconditional  correlations  of  temperature  and  jitter  so 
that  we  do  not  have  to  wait  for  a  particular  “  event”  to  trigger  the  prediction 
from  a  conditional  average. 

The  usual  expression  for  the  expected  value  of  the  jitter  given  a  temperature 
measurement  T,  can  be  expanded  in  a  power  series  in  T.  It  has  been  shown  that 
the  linear  term  in  the  expansion  is  sufficient  to  describe  coherent  organized 
structure  in  flow.  So  our  estimate  then  takes  the  form  of  the  linear  combination 
of  weighted  temperature  signals.  We  need  to  calculate  the  weights,  that  is  the 
“A”  coefficients.  By  minimizing  the  mean-square  error  between  the  estimate 
and  the  jitter  we  obtain  an  expression  that  depends  upon  the  correlations 
between  the  measured  temperature  and  jitter  signal  and  also  depends  upon 
correlations  between  the  temperature  signals.  We  solve  for  the  A’s  by  solving  a 
set  of  linear  equations  shown  here  in  matrix  form  for  two  temperature 
measurements  T1  and  T2. 1  want  to  emphasize  that  these  calculations  are 
accomplished  before  correction  occurs.  This  is  a  kind  of  ’’training”  done  ahead 
of  time  and  is  based  on  actual  measurements  of  the  jitter. 


LSE  Cont’d 
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Examine  the  temporal  delays,  x,  that  maximize  the 
probability  that  the  estimated  jitter  is  best  represented  by 
the  measured  temperature 

•  Choose  “lags"  ti  to  maximize  the  Temperature-Jitter  correlations 

r,:5£/P(r,4) 

•  Evaluate  Temperature-Jitter  correlations  at  t; 
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•  Evaluate  Temperature-Temperature  correlations  at  v, 


In  order  to  find  the  best  values  of  the  correlations  to  put  into  the  matrix 
expression  we  have  to  find  the  particular  time  tau  at  which  to  evaluate  the 
correlations.  The  best  tau  is  the  time  lag  at  which  the  temperature-jitter 
correlation  is  maximized.  We  evaluate  both  the  temperature-jitter  correlation 
and  the  temperature-temperature  correlations  at  this  time  lag  and  input  the 
results  into  the  matrix  expression. 


SUP=  mathematical  expression  for  supremum  (maximum  value) 


LSE  References: 

T.  Luna,  "Linear  Stochastic  Estimation  of  Optical  Beam  Deflection 
through  a  Heated  Round  Turbulent  Jet,"  MS  Thesis,  Dept  of  Mechanical 
Engineering,  University  of  New  Mexico,  1996. 

T.  Luna,  C.  R.  Truman,”Linear  stochastic  estimation  of  optical  beam  deflection 
through  a  heated  jet,”  AIAA  paper  97-0072,  35th  Aerospace  Sciences  Meeting, 
Reno,  NV  Jan  6-10,  1997. 


This  is  the  LSE  procedure  shown  in  block  diagram  form  where  two 
temperature  signals  and  a  measured  optical  jitter  signal  are  input.  The  “lags” 
are  found  from  the  time  at  which  the  T-J  correlations  are  maximum.  The  T-T 
correlations  are  then  evaluated  at  this  lag.  The  set  of  linear  equations  defined 
by  the  matrix  expression  given  previously  are  solved  for  the  coefficients  Aj 
“weights”.  This  is  the  way  LSE  is  “trained”  on  the  flow. 

These  weights  and  lags  are  then  used  in  the  actual  correction  stage. 

To  make  the  actual  prediction  the  weights  and  lags  are  applied  to  two  new 
temperature  measurements  made  at  the  time  of  correction.  The  linear 
combination  of  the  weighted  signals  at  the  appropriate  time  delays  creates  the 
prediction. 


In  the  coirection  stage  the  weights  and  lags  are  given.  These  weights  and  lags 
are  applied  to  new  temperature  measurements  that  are  taken  at  the  time  of 
correction  to  produce  the  estimated  jitter  signal  J'',  which  is  a  linear  summation 
of  the  ATi’s,  i.e.  J''=  AiTi  +A2T2 

This  diagram  is  to  illustrate  the  difference  between  the  simulations  that  we 
have  used  to  examine  how  well  LSE  can  do  to  correct  our  jitter  signal,  and  the 
experiments  we  have  actually  performed. 

In  the  simulations  the  2  temperature  signals  and  the  jitter  signal  are  used  to 
create  the  weights  and  lags,  the  weights  and  lags  are  then  applied  to  the 
prediction  algorithm  in  which  the  same  temperature  signals  are  used.  The 
resulting  estimated  jitter  signal  is  then  compared  to  the  measured  jitter  signal 
used  in  the  LSE  to  determine  the  residual  jitter. 

In  the  experiment,  the  weights  and  lags  are  applied  to  the  new  temperature 
signals  measured  during  the  correction.  The  jitter  estimate  is  used  to  control  a 
steering  mirror  and  the  residual  jitter  is  measured.  To  determine  the  amount  of 
jitter  correction  that  has  taken  place,  the  measured  residual  is  compared  to  the 
jitter  signal  used  to  create  the  weights  and  lags. 


The  next  five  slides  show  the  simulation  and  experimental  results  illustrating 
that  the  LSE  prediction  can  be  effective  in  reducing  beam  jitter. 


Simulation  Results — Forced  Flow 
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In  simulation,  2  temperature  signals  and  a  measured  jitter  signal  are  measured 
and  used  to  produce  the  weights  and  lags  describing  the  flow  at  a  downstream 
position  of  2  nozzle  diameters  (~  2.5  cm).  The  weights  and  lags  are  then 
applied  to  the  same  temperature  measurements  to  predict  the  jitter.  When 
compared  to  the  measured  tilt,  the  predicted  jitter  signal  would  remove  59%  of 
the  RMS  jitter  signal. 


Comparison  to  experimental  results.  In  this  case  the  predicted  jitter  is  used  to 
move  a  piezo  mirror  to  compensate  for  the  jitter.  The  results  show  a  41%  (on 
average)  correction  compared  to  nearly  60%  in  simulation. 


Simulation  Results — Unforced  Flow 
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We  also  looked  at  an  unforced  flow  case  where  the  flow  structures  are 
intermittant  rather  than  periodic  as  in  the  forced  case.  As  expected  the 
prediction  is  not  as  accurate  in  the  case  averaging  about  17%  as  compared  to 
almost  60%  for  the  forced  case.  Although  the  reduction  in  prediction  accuracy 
is  not  completely  understood  and  these  results  are  preliminary,  it  may  in  part 
be  due  to  smaller  correlations  between  the  jitter  and  the  temperature 
measurements  as  well  as  a  lack  of  symmetry  in  the  unforced  flow  along  the 
beam  path.  In  this  case  only  2  measurements  do  not  predict  events  along  the 
entire  beam  path  very  well. 
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Comparison  to  experimental  results.  In  the  unforced  case  experimental 
correction  is  about  16%  compared  to  the  simulated  value  of  17%. 


Simulation  and  Experimental  Results  for  2-pt  LSE  vs.  4-pt  LSE 
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Comparison  of  2-point  LSE  to  4-Point  LSE  .  As  expected  prediction  accuracy 
increases  as  the  number  of  point  measurement  inputs  increases. 


Implications  for  Adaptive  Optics 


•  Extension  of  thin  beam  to  spatially  resolved 
correction 

-  Tomography  system  allows  multiple  non- 
intrusive  point  measurments 

-  MEMS  devices,  concept  of  multiple  aperture 
adaptive  optics 

•  Combination  of  predictive  and  closed  loop 
correction 


High  Speed  Optical  Tomography 

2-D  temperature  measurements  at  5  KHz 


•spatial  and  temporal  measurement  can  allow  higher  resolution 
predictionover  larger  spatial  areas 

•suitable  for  MEMS  or  other  multiple  aperture  adaptive  optics  systems 
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In  other  work  at  the  Phillips  Lab  we  have  achieved  high  spatial  and  temporal 
resolution  optical  tomography  for  measuring  the  temperature  of  dynamic 
flows.  Although  there  is  no  time  to  discuss  the  details  of  this  system,  what  we 
get  from  the  system  are  2  dimensional  measurements  of  the  temperature.  In 
these  tomographic  reconstructions  we  represent  the  temperature  as  a  grey  level 
so  that  it  looks  like  a  picture.  But  we  really  have  a  64  x  64  spatially  resolved 
temperature  measurement  taken  with  a  speed  of  5  KHz. 

Using  the  spatial  and  temporal  fidelity  of  this  measurment  could  allow  the 
predicitive  algorithms  to  be  applied  to  many  spatial  points  at  once. 

Tomography  references: 
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measurement  of  optical  turbulence  by  tomographic  imaging  of  a  heated  air 
flow,"  SPIE  Procedings  vol.  2827,  p.  130-141,  Digital  Image  Recovery  and 
Synthesis  HI,  Denver,  CO,  Aug  5-6,  1996. 

L.  McMackin,  B.  Masson,  N.  Clark,  K  Bishop,  R.  Pierson  and  E.  Chen, 
"Hartmann  wavefront  sensor  studies  of  dynamic  organized  structure  in  flow 
fields,"  AIAA  Journal  33  (11),  2158-2164,  1995. 

C.  R.  Truman,  T.  Luna,  L.  McMackin,  B.  Masson.  K.  P.  Bishop,  E.  Y.  Chen, 
"Optical  tomographic  study  of  the  effect  of  excitation  of  a  heated  round  jet," 
AIAA  Paper  96-2323,  Plasma  Dynamics  and  Lasers  Conference,  New  Orleans, 
June  1996. 


Improve  Current  Control  Systems 
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Example  of  how  predictive  procedures  (such  as  LSE)  can  be  used  to  improve 
the  correction  results  from  closed  loop  adaptive  optics  systems. 

In  this  case  rather  than  simply  feeding  back  the  residual  disturbance  to  the  tilt 
sensor,  the  error  is  also  used  to  drive  an  adaptive  prediction  system  in  which 
the  temperature  signals  as  well  as  the  feedback  signal  are  used  to  compute  the 
weights.  In  this  way  both  feedback  and  prediction  are  used  to  adaptively 
correct  for  dynamic  aberrations. 


Fiber  Phase  Conjugator 
as  Wave  Front  Corrector 
for  Industrial  Laser  Systems 
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Fibers  with  DLD^erent  Core  Diameters 


Choice  of^active  medium 
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PHASE  CONJUGATE  MIRROR  BASED  ON  NDiYAG 
CRYSTAL  FOR  HIGH  POWER  LASER  SYSTEMS 
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Reflection  coefficient 
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Experimental  depend¬ 
ences  of  the  reflection 
coefficient  of  the  SPEC 
system  at  the  peak  power 
on  the  total  gain  of 
"nonlinear"  and  FL 
amplifiers  with  the  peak 
power  of  the  incident 
beam  0.075  W  (a),  and 
on  the  peak  power  of  the 
incident  beam  with  the 
total  gain  7.2  at  the 
single  pulse  regime  (b). 
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Phase  conjugation  of  the  phase  -  distorted  beam 


Transverse  ctructure  of  recovered  image  "net" 
for  the  same  pumping  beam  diameters  (a) 
and  for  the  different  diameters  (b). 
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Three-loop  scheme 
gain  a/  =  15 
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Different  schemes  of  laser  with  distributed  feedback  on  the 
population  inversion  gratings. 


ANISOPLANATIC  EFFECT  IN 


Residual  errors  for  two  separated  points  in 

the  object  plane  Extended  object  correction 
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correlation  of  Zemike  coefficients 
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•  Residual  error  as  a  function  of  the  distance  from  the  referent 

I  p  point 


1,2,3  -  corrector  in  the  centre  of  the  layer 

1  -  Zc=Zs=0.8 

2  -  Zc=Zs=0.5 

3  -  Zc=Zs=0.2 

4,5,6  -  corrector  in  the  input  aperture  plane 

4  -  2s=0.8 

5  -  Zs=0.5 

6  -  Zs-0.2 


